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FOREWORD 
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to  1  May  1968.  The  work  was  conducted  under  Project  1472, 
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"Optimization  of  Dynamics  Data  Analysis.  "  Technical  administration 
by  AFFDL  was  provided  by  Mr.  Charles  E.  Thomas  and  Mrs.  Phyllis 
Bolds  (FDDS). 

The  work  was  performed  by  MAC  personnel  under  the  direction 
of  Mr.  Allan  G.  Piersol  with  the  assistance  of  Mr.  W.  Frank  van  der 
Laan.  Principal  contributors  included  Dr.  Richard  L.  Barnoski 
(structural  dynamics).  Dr.  Sung  C.  Choi  (statistics),  Mr.  Loren  D. 
Enochson  (digital  computing),  Mr.  Ronald  D.  Kelly  (analog  data  anal¬ 
ysis),  and  Dr.  Prichard  H.  White  (vibroacoustics).  The  contractor's 
report  number  is  MAC  Report  No.  616-05. 
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ABSTRACT 


Techniques  for  predicting  the  vibration  environments  of  future 
aircraft  based  upon  statistical  extrapolations  from  data  measured  on 
past  aircraft  are  investigated.  As  a  first  step,  principal  sources  of 
aircraft  vibration  are  identified,  and  analytical  relationships  for  the 
resulting  vibration  environment  are  approximated.  Available  AFFDL 
data  are  then  summarized  and  evaluated.  Extensive  statistical  studies 
are  performed  on  the  available  data  to  investigate  the  average  properties 
of  aircraft  vibration  among  the  three  orthogonal  directions,  various 
structural  zones,  various  aircraft  models,  and  various  aircraft  groups. 
The  available  data  are  also  used  to  study  the  spatial  distribution  of 
vibration  levels  within  a  given  structural  zone.  Specific  prediction 
models  are  then  suggested  and  regression  analysis  procedures  to  arrive 
at  conservative  prediction  levels  are  detailed.  The  suggested  techniques 
are  illustrated  using  available  AFFDL  data.  Procedures  for  deriving 
vibration  test  specifications  based  upon  environmental  vibration  predic¬ 
tions  are  suggested.  Finally,  possible  extension  of  the  techniques  to  the 
prediction  of  internal  acoustic  noise  are  discussed. 

This  abstract  is  subject  to  special  export  controls  and  each  trans¬ 
mittal  to  foreign  governments  or  foreign  nationals  may  be  made  only 
with  prior  approval  of  the  Air  Force  Flight  Dynamics  Laboratory, 
Wright- Patterson  AFB,  Ohio. 
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1.  INTRODUCTION 


The  most  critical  step  in  the  derivation  of  vibration  and  acoustic 
design  criteria  and/or  test  specifications  for  modern  aircralt  is  the  pre¬ 
diction  of  the  anticipated  flight  vibration  and  acoustic  environ: 7 lents. 

This  usually  involves  the  prediction  of  power  spectra  (or  som-  other 
measure  of  frequency  composition)  for  the  environments.  Cor  .lider 
first  the  problem  of  predicting  flight  vibration  environmi  nts. 

There  are  two  fundamental  ways  to  approach  the  vibration  predic¬ 
tion  problem.  For  the  purpose  here,  these  two  ways  will  be  referred 
to  as  (1)  analytical  prediction  procedures  and  (2)  empirical  extrapolation 
prediction  procedures.  Analytical  prediction  procedures  inch  de  all 
those  techniques  for  predicting  flight  vibration  environments  which  are 
based  upon  calculating  or  measuring  the  response  of  a  derived  structural 
model  to  an  assumed  excitation  function.  Such  prediction  proc  idures 
have  been  classified  by  Gray  and  Piersol  [1]  into  four  categ  ,rif.s,  as 
follows : 

1.  Classical  (mathematical  model)  approach 

2.  Multiple  input  model  approach 

3.  Physical  model  approach 

4.  Statistical  energy  approach 

Analytical  prediction  procedures  have  been  used  with  varying  degrees  of 
success  in  the  past.  They  all  have  the  common  characteristic,  however, 
of  requiring  detailed  information  about  the  structural  design  of  the  air¬ 
craft  in  question,  and  the  excitation  functions  to  which  it  will  be  sub¬ 
jected.  Unfortunately,  such  information  is  rarely  available  in  the 
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necessary  detail  at  the  time  >)vhen  vibration  predictions  are  most  de¬ 
sired  (early  in  the  design  phase). 

Empirical  extrapolation  prediction  procedures  include  those  tech¬ 
niques  which  are  based  upon  studies  of  data  collected  during  flight  tests 
of  previous  aircraft.  There  are  two  basic  types  of  extrapolation  tech¬ 
niques.  The  predictions  may  be  based  upon  data  from  a  single  vehicle 
of  similar  design  (specific  extrapolation),  or  they  may  be  based  upon 
pooled  data  from  one  or  more  general  vehicles  (general  extrapolation). 

The  most  commonly  used  specific  extrapolation  procedures  are  those 
suggested  by  Condos  and  Butler  [2],  barrett  [3],  and  Winter  [4] .  All 
three  procedures  employ  a  scaling  formula  to  extrapolate  vibration  data 
measured  on  a  previous  vehicle  to  predict  the  vibration  of  a  new  vehicle 
of  similar  design  based  upon  differences  in  the  excitation  pressure  levels 
and  structural  surface  weight  densities. 

General  extrapolation  procedures  are  based  upon  empirically 
derived  correlations  between  the  average  vibration  response  character¬ 
istics  for  a  general  class  of  structure  and  one  or  more  parameters  re¬ 
lated  to  the  excitation  forces,  structural  properties,  and/or  flight  con¬ 
ditions  for  the  aircraft  of  interest.  A  number  of  such  procedures  have 
been  proposed  in  recent  years  for  missiles  and  spacecraft,  a^  well  as 
aircraft.  One  of  the  earliest  to  appear  in  the  literature  was  the  method 
suggested  by  Mahaffey  and  Smith  [5],  which  is  based  upon  observed  corre¬ 
lations  between  structural  vibration  and  jet  engine  acoustic  noise  for  data 
collected  from  the  B-58  aircraft.  These  same  data  were  modified  by 
Brust  and  Himelblau  [6]  to  develop  a  prediction  rule  for  the  SKYBOLT 
missile.  Similar  studies  were  performed  by  Eldred,  Roberts,  and  White 
[7]  using  data  collected  from  the  SNARK  missile.  Curtis  [8]  developed 
an  empirical  prediction  rule  by  observing  correlations  between  aircraft 


structural  vibration  and  free  stream  dynamic  pressure  using  data 
collected  from  the  F-AU,  B-59,  F-101,  and  F-106  aircraft.  From 
JUPITER  and  TITAN  I  missile  vibration  and  acoustic  data  collected 
during  static  firings,  Franken  [9]  developed  a  procedure  which  predicts 
radial  skin  vibration  for  missiles  as  a  function  of  rocket  engine  acoustic 
noise  and  skin  surface  weight  density.  A  similar  technique  was  de¬ 
veloped  by  Winter  [4]  for  predicting  ring  frame  and  stringer  vibration 
for  missiles  and  spacecraft  based  upon  JUPITER,  TITAN,  MINUTEMAN, 
SKYBOL.T,  and  GEMINI  vibration  and  acoustic  measurements. 

All  of  the  above  general  extrapolation  procedures  are  similar  in 
that  they  permit  the  prediction  of  structural  vibration  levels  in  a  future 
flight  vehicle  without  a  detailed  knowledge  of  the  specific  structural  de¬ 
sign,  or  the  need  for  detailed  data  from  a  previous  vehicle  of  similar 
design.  The  advantage  of  such  procedures  is  clear.  They  can  be  readily 
applied  to  anticipated  flight  vehicles  even  at  the  preliminary  design  stage 
before  the  detailed  structural  design  has  been  established.  The  funda¬ 
mental  disadvantage  is  equally  clear.  Since  the  procedures  do  not  use 
detailed  structural  information  for  the  flight  vehicle  design  in  question, 
or  specific  data  for  a  similar  vehicle,  they  do  not  provide  the  potential 
accuracy  which  an  analytical  procedure  or  a  specific  extrapolation  pro¬ 
cedure  could  theoretically  produce  under  ideal  conditions.  In  the  opinion 
of  many  contemporary  environmental  engineers,  however,  the  fundamental 
advantage  of  general  extrapolation  prediction  procedures  far  outweighs 
their  disadvantage.  This  fact  is  clearly  established  by  the  current  wide¬ 
spread  use  of  such  procedures. 

Noting  the  current  availability  of  previously  developed  general 
extrapolation  procedures,  an  obvious  question  arises  at  this  point; 
namely,  is  there  a  need  for  an  improved  general  extrapolation  procedure 
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for  predicting  aircraft  vibration  environments?  It  is  believed  that  the 
answer  to  this  question  is  yes.  A  recent  survey  of  the  better  known 
flight  vehicle  vibration  prediction  techniques  [4]  indicates  that  all  of  the 
general  extrapolation  procedures  reviewed  have  limitations  and  de¬ 
ficiencies  of  one  form  or  another.  The  principal  deficiencies  are  as 
follows : 

1.  All  of  the  procedures  (with  the  notable  exception 
of  the  Curtis  procedure  [5])  require  predictions 
for  the  excitation  environment.  The  resulting 
vibration  predictions  can  be  no  better  than  the  ex¬ 
citation  predictions. 

2.  All  of  the  procedures  (again  excluding  the  Curtis 
procedure)  were  developed  using  acoustic  induced 
vibration  data.  Such  procedures  are  fully  useful 
only  for  the  prediction  of  takeoff  (or  liftoff)  vibration 
environments . 

3.  Most  of  the  procedures  were  developed  from  data  for 
only  a  few  (in  some  cases,  only  one)  flight  vehicles. 

This  clearly  limits  the  generality  of  the  procedures. 

4.  None  of  the  procedures  include  provisions  for  pre¬ 
dicting  periodic  contributions  to  the  vibration  environ¬ 
ment  from  the  rotation  of  jet  engine  shafts  and/or 
auxiliary  equipment. 

5.  The  statistical  techniques  used  in  the  development  of 
the  procedures  were  not  always  as  thorough  as  would 
appear  to  be  warranted  by  the  importance  of  the  prob¬ 
lem. 


The  above  deficiencies,  along  with  others,  have  tended  to  limit  the 
effectiveness  of  the  various  procedures  in  practice.  Available  compari¬ 
sons  between  predictions  and  actual  measured  data  [4]  indicate  that 
errors  in  the  predicted  power  spectra  of  20  dB  or  more  are  common. 
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In  light  of  the  above  discussions,  it  is  believed  that  the  develop¬ 
ment  of  an  improved  and  more  efficient  general  extrapolation  prediction 
procedure  is  a  worthwhile  task.  For  the  case  of  aircraft  vibration  pre¬ 
dictions,  AFFDLi  is  in  a  unique  position  to  pursue  such  a  task  because 
of  their  large  and  ever  growing  library  of  aircraft  vibration  data. 

Now  consider  the  problem  of  predicting  flight  acoustic  environ¬ 
ments  (inside  the  aircraft).  Considerable  attention  has  been  given  in 
recent  years  to  the  prediction  of  the  external  acoustic  noise  generated 
by  jet  and  rocket  engines,  as  well  as  the  aerodynamic  noise  generated 
by  boundary  layer  turbulence  during  flight  [10].  As  previously  noted, 
such  information  is  needed  as  an  intermediate  step  in  the  prediction  of 
flight  vibration  environments.  Except  for  the  case  of  commercial  air¬ 
craft,  however,  somewhat  less  attention  has  been  given  to  the 
problem  of  predicting  internal  acoustic  noise  environments.  This 
relative  lack  of  interest  is  unquestionably  due  in  large  part  to  the  wide¬ 
spread  (and  in  most  cases,  accurate)  belief  that  the  internal  acoustic 
noise  environment  in  aircraft  is  far  less  damaging  to  equipment  than  the 
vibration  environment.  Furthermore,  since  the  crew  of  noncommercial 
aircraft  generally  wear  protective  head  gear,  relatively  high  acoustic 
noise  levels  can  be  permitted  inside  the  aircraft  without  posing  a  serious 
hazard  to  the  crew  members. 

For  the  case  of  commercial  aircraft,  internal  acoustic  predictions 
are  usually  based  upon  calculations  for  the  attenuation  of  aeroacoustic 
loads  by  the  aircraft  structure  and  sound  proofing,  combined  with  esti¬ 
mates  for  the  noise  generated  by  the  airconditioning  system  and  other 
noise  producing  equipment.  Because  of  the  stringent  competitive  re¬ 
quirements  for  passenger  comfort,  the  predictions  must  be  relatively 
accurate.  For  the  case  of  noncommercial  aircraft,  however,  it  appears 
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that  acoustic  predictions  with  the  accuracy  provided  by  a  general  ex¬ 
trapolation  procedure  might  be  adequate.  The  statistical  techniques 
required  to  develop  an  extrapolation  procedure  for  acoustic  predictions 
would  be  basically  the  same  as  those  required  to  develop  a  procedure 
for  vibration  predictions.  Hence,  it  is  believed  that  this  constitutes  a 
worthwhile  task  to  pursue  along  with  the  development  of  a  general  ex¬ 
trapolation  procedure  for  aircraft  vibration  predictions. 

The  primary  purpose  of  the  studies  reported  herein  is  to  formu¬ 
late  a  well  defined  program  for  the  development  of  an  improved  extrapo¬ 
lation  procedure  for  the  prediction  of  aircraft  flight  vibration  environ¬ 
ments.  The  intent  is  that  AFFDL  will  implement  the  program  using 
AFFDL  collected  flight  vibration  data  and  computer  facilities. 

Secondary  objectives  of  the  studies  include  (a)  a  general  evaluation 
of  currently  available  AFFDL  aircraft  vibration  data  and  (b)  the  in¬ 
vestigation  of  improved  techniques  for  converting  aircraft  flight  vibra¬ 
tion  predictions  into  appropriate  test  levels  and  durations  for  aircraft 
component  vibration  test  specifications. 

In  order  to  establish  proper  statistical  procedures  and  illustrate 
their  use,  it  was  necessary  during  the  course  of  these  studies  to  perform 
a  considerable  amount  of  data  analysis.  For  convenience  and  clarity, 
only  summaries  of  pertinent  results  are  included  in  this  report.  The  de¬ 
tailed  results  of  the  analysis  are  presented  under  separate  cover  . in 
Appendix  A.  Also  as  part  of  the  work  reported  herein,  a  digital  com¬ 
puter  program  for  the  efficient  computation  of  power  spectra  using  fast 
Fourier  transform  techniques  was  developed  and  delivered  to  AFFDL. 

The  documentation  for  this  fast  Fourier  transform  power  spectrum  pro¬ 
gram  is  also  presented  under  separate  cover  in  Appendix  B. 


2.  PROBLEM  FORMULATION  AND  APPROACH 


The  basic  approach  to  be  pursued  herein  centers  around  the 
development  of  a  linear  model  for  the  power  spectrum  of  aircraft  flight 
vibration.  In  general  terms,  the  model  will  be  of  the  form 


G(f)  =  A^(f)  +  A^(f)  Xj  +  A^({)  +  .  .  .  +  Aj^(f) 

(1) 

N 

=  A  (f)  +  ^  A  (f)  X 
i=l 


where 


G(f)  =  the  power  spectral  density  for  the  vibration 
response 

Xj  =  the  ith  independent  variable 

A.(f)  =  the  weighting  factor  for  the  ith  independent 
variable 

AQ(f)  =  the  power  spectral  density  for  the  residual  vibration 
(if  any)  when  all  independent  variables  equal  zero 

The  variables,  X^  (i  =  1,  2,  3,  .  .  .  ,  N)  ,  ideally  would  cover  all 
factors  which  influence  the  flight  vibration  environment.  Included  would 
be  pertinent  descriptive  parameters  of  the  aircraft  structure,  the  engine 
operating  conditions,  and/or  the  aircraft  flight  conditions.  For  example, 
the  X.  variables  might  be  X  =  structural  weight  density  (w)  ,  X  =  en- 
gine  exhaust  gas  velocity  (V  ) ,  and  X,  =  dynamic  pressure  (q)  .  Note 
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that  a  single  variable  might  be  a  power,  product,  and/or  quotient  of 

2  8 

several  parameters.  For  example,  =  (q/w)  and  .  Hence, 

although  Eq.  (1)  is  a  linear  model,  it  can  be  applied  to  nonlinear  re¬ 
lationships  as  long  as  they  are  anticipated  and  properly  incorporated 
into  the  independent  variables. 

The  weights,  A^(f)  (i  =  1,  2,  3 . N) ,  establish  the  relative 

contribution  of  each  variable  to  the  vibration  power  spectrum.  For  ex¬ 
ample,  if  Xj  =  w  and  Aj(f)  =  0,  this  would  mean  that  the  structural 
weight  density  has  no  influence  on  the  vibration  power  spectrum  at  fre¬ 
quency  f .  The  AQ(f)  weight  accounts  for  any  contributions  to  the 
vibration  power  spectrum  which  occur  when  X^  =  0  (i  =  1,  2,  3,  .  .  ,  ,  N) 

The  development  of  a  model  of  the  form  given  in  Eq.  (1)  involves 
two  general  problems.  The  first  is  to  select  the  variables  to  be  used 
in  the  model,  and  the  second  is  to  determine  appropriate  values  for 
the  weights.  These  two  problems  along  with  statistical  considerations 
and  the  specific  approach  which  hopefully  will  solve  them  are  now  dis¬ 
cussed. 

2.  1  SELECTION  OF  VARIABLES 

The  selection  of  appropriate  independent  variables  could  be 
approached  by  purely  empirical  procedures.  This  would  be  done  by 
simply  guessing  at  all  possible  factors  which  might  influence  structural 
vibration.  Those  factors  which  are  actually  related  to  structural  vibra¬ 
tion  in  a  statistically  significant  way  would  then  be  determined  through 
a  multiple  correlation  study  of  all  available  vibration  data.  Such  an 
approach,  however,  is  considered  unsuitable  for  the  problem  of  concern 


It  is  believed  that  a  superior  approach  is  to  select  the  variables 
based  principally  upon  theoretical  knowledge  and  well-established  ex¬ 
perience.  This  second  approach  provides  important  advantages.  First, 
assuming  sufficient  theoretical  and  experimental  information  is  avail¬ 
able,  it  increases  the  likelihood  of  selecting  only  those  variables  which 
are  relevant  to  the  vibration  environment.  Second,  it  permits  nonlinear 
relationships,  as  well  as  product  and  quotient  relationships,  for  various 
factors  to  be  anticipated  and  properly  included  into  £q.  (1)  as  a  single 
variable  in  the  linear  model.  The  elimination  of  irrelevant  variables 
and  nonlinear  relationships  in  Eq.  (1)  will  greatly  increase  the  statistical 
significance  of  the  regression  analyses  required  to  calculate  the  A^(f) 
weights.  Of  course,  there  is  the  possibility  that  relevant  variables  may 
be  omitted  or  that  nonlinear  relationships  for  various  factors  may  be 
improperly  interpreted.  These  risks,  however,  should  be  minimal  for 
two  reasons.  First,  theoretical  and  experimental  definitions  for  air¬ 
craft  vibration  sources  and,  hence,  the  factors  which  influence  flight 
vibration,  are  readily  available  and  quite  thorough.  Second,  special 
tests  for  significant  correlation  and  linearity  of  all  selected  variables 
can  and  will  be  performed  as  part  of  the  regression  analysis  to  establish 
the  A,(f)  weights. 

2.  2  DETERMINATION  OF  WEIGHTS 

The  A^(f)  weights  in  Eq.  (1)  will  be  determined  using  a  conventional 
regression  analysis  procedure.  The  general  procedure  is  as  follows. 
Assume  that  M  different  aircraft  flight  vibration  measurements  are  avail¬ 
able  for  various  different  combinations  of  structural  locations,  engine 
operating  conditions,  and/or  flight  conditions.  Hence,  M  power  spectra 
can  be  computed,  and  each  can  be  considered  an  estimate  of  the  power 
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spectrum,  G(f)  ,  given  by  Eq.  (1).  That  is,  the  kth  measured  spectrum, 
/s 

G,  (f)  (k  =  1,  2,  3,  .  .  .  ,  M)  can  be  considered  an  estimate  (denoted  by  the 

iC 

hat  ^  )  for 


G^(f)  =  A^(f)  +  A^(f)  +  A2(f)  +  .  .  .  ^  Aj^(f) 

(2) 

=  A„(£)  .  Z  A.(£)  X.^  ;  k=l,2,  3 . M 


where  X.,  is  the  value  of  the  ith  variable  for  the  kth  measurement, 
ik 

The  problem  is  to  solve  for  values  of  A.(f)  which  will  minimize  the 

variance  of  the  measured  values  G,  (f)  about  the  model  values  G,  (f) . 

k  ^  k 

Assume  that  the  spectral  density  measurements, 
normally  distributed  with  a  mean  value  of  Gj^(f) ,  as  given  in  Eq.  (2) . 


Under  this  assumption,  the  best  minimum  variance  unbiased  estimates 
for  the  weights  A^(f)  are  determined  as  follows.  Let  Q(f)  denote  the 
sum  of  the  squared  deviations  at  a  given  frequency,  as  follows. 


Q(f)  = 


M  »• 

=  S  [ 

k=l  L 


12 

G, _(f)  -  G,  (f) 


(3) 


Using  Eq.  (2)  and  dropping  the  (f)  notation  for  clarity,  it  follows  that 


M 
k=  1 


^Nk 


(4) 
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HH 


By  taking  the  derivatives  of  Q  with  respect  to  (i  -  0,  1,  2,  .  .  .  ,  N) 
and  equating  them  to  zero,  the  following  set  of  sini iltaneous  equations 
is  generated. 


“^0  ^1  ^  ^Ik 

k 


k  k 


A  A  __  2 

X,,  +  A,  ^  X,,  + 
0-^  Ik  Ik 

k  k 


^Ik  ^Nk  ■  ?  ^Ik  S 


■^0  ^  ^2k  “^1  2  ^2k  ^Ik 

k  k 


2  ^2k  ^Nk  ^  ^2k  °k 
k  k 


(5) 


k  k 


Nk 


k 


The  relationships  in  Eq.  (5)  are  known  as  the  "normal  equations"  in 

A 

multiple  regression  theory.  The  terms,  (i  =  0,  1,  2,  ....  N)  ,  are 

estimates  for  the  A,  weights.  The  problem  now  is  to  solve  Eq.  (5)  for 
A  ^ 

the  A.  terms.  To  accomplish  this,  let 


B 

r,  s 


M 

S  X 

k=l 


rk  ^sk 


r,  s  =  1,  2,  3,  .  .  .  ,  N 


(6) 


11 


Excluding  the  first  equation  and  the  first  colunnn  from  each  equation, 
Eq.  (5)  may  now  be  written  in  matrix  form,  as  follows. 


(7) 


Equation  (7)  solves  for  the  A.  weight  for  i  =  1,  2,  3,  ....  N.  The 
weight  can  be  solved  for  using  the  first  equation  in  Eq.  (6)  after  the 
other  weights  have  been  determined.  Computer  programs  for  perform¬ 
ing  multiple  regression  analyses  are  widely  available. 


2.  3  STATISTICAL  CONSIDERATIONS 

Having  estimated  the  weights  for  the  N  selected  variables,  a  pre¬ 
diction  model  for  the  power  spectrum  of  aircraft  flight  vibration  can  now 
be  written,  as  follows. 


N 


G(f)  =  A  (f)  =  2  A,(f)  X 
i=l 


(8) 
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At  this  point,  it  should  be  made  clear  that  the  development  of  a  "perfect" 
prediction  model  for  aircraft  vibration  is  beyond  practicality.  There  are 
certain  factors  (particularly  those  related  to  structural  characteristics) 
which  significantly  influence  '/ibration,  but  which  cannot  be  conveniently 
defined  in  terms  of  simple  variables.  In  other  words,  there  will  be 
some  pertinent  terms  missing  in  Eq.  (8).  This  will  appear  as  scatter 
in  the  actual  power  spectra  about  the  predicted  power  spectra  determined 
using  Eq.  (8). 

To  illustrate  this  point,  assume  the  model  consists  of  only  one 
variable,  such  that 


G(f)  =  A(f)  X 


where  the  regression  equation  used  to  estimate  the  weight  A(f)  is 


(10) 


Now  assume  the  M  power  spectra  mea  ired  for  the  various  values  of  X 
are  compared  to  the  prediction  model  given  by  Eq.  (9).  The  results  at 
a  specific  frequency  might  be  as  shown  in  Figure  1.  The  scatter  in 
Figure  1  may  be  thought  of  as  an  indication  of  the  "efficiency"  of  the 
prediction  model.  If  X  were  the  only  variable  relevant  to  the  vibration 


environment,  there  would  be  no  scatter;  that  is,  all  measurements 
would  fall  on  the  prediction  line  (assuming  no  measurement  error  is 
present).  This  would  constitute  a  perfectly  efficient  model. 

A  convenient  measure  for  the  efficiency  of  the  prediction  model 
is  given  by  the  "correlation  coefficient,  "  y  »  ior  the  regression  analysis. 
The  correlation  coefficient  is  a  number  bounded  by  zero  and  unity, 

2 

(0  <  y  <  1),  where  unity  indicates  a  perfect  model.  The  quantity  1  -  y 
is  a  measure  of  the  "power"  contributed  to  the  observed  vibration  by 
variables  which  have  not  been  included  in  the  model.  For  correlation 
coefficients  of  less  than  unity,  the  scatter  of  the  measured  data  about 

A 

the  model  can  be  described  by  a  standard  deviation  for  the  measured  G 
values  at  various  values  of  X.  Assuming  some  specific  distribution 
function  for  this  scatter,  prediction  limits  for  future  measurements  of 
the  vibration  environment  at  any  value  of  X  can  be  determined.  For  the 
case  of  a  flight  vehicle  vibration  model,  conservative  estimates  are  re¬ 
quired.  Hence,  some  upper  limit  with  a  relatively  large  probability  of 
exceeding  future  vibration  levels  will  be  required  in  practice. 

Other  analyses  which  shovild  be  performed  during  the  regression 
study  include  a  test  for  significance  of  the  calculated  correlation  coeffi¬ 
cients  and  a  test  for  linearity  of  the  model.  The  former  can  be  per¬ 
formed  using  a  Fisher  "Z"  transformation  and  the  latter  using  an  analysis 
of  variance  test. 

These  varr'ous  statistical  analysis  procedures  will  be  outlined  and 
illustrated  in  Section  8. 

2.  4  OUTLINE  OF  SUGGESTED  APPROACH 

The  specific  steps  which  will  be  pursued  to  develop  a  prediction 
model  for  aircraft  vibration  using  AFFDL  data  are  as  follows. 
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1.  Development  of  Vibration  Relationships 


This  initial 


step  involves  three  requirements.  The  first  is  to  identify  the  various 
sources  of  aircraft  vibration  and  assess  their  relative  importance. 

The  second  is  to  provide  a  preliminary  definition  of  the  structural  ex¬ 
citation  parameters  which  theoretically  influence  aircraft  vibration 
environments.  The  third  is  to  analytically  relate  the  primary  sources 
of  vibration  to  basic  flight  and/or  engine  performance  parameters. 
Section  3  covers  this  subject. 


2.  Summary  of  Available  Data  —  The  flight  vibration  and 
acoustic  data  currently  available  for  study  must  now  be  defined  and 
collected.  Fortunately,  AFFDL  has  carefully  documented  and  pub¬ 
lished  the  results  of  all  flight  vibration  surveys  performed  by  their 
activity.  Hence,  this  first  step  requires  only  a  summary  of  this 
documentation,  which  is  presented  in  Section  4. 


3.  Preliminary  Evaluation  of  Available  Data  —  The  final  pro¬ 
cedure  to  be  recommended  in  this  report  will  require  no  small  amount 
of  digital  computer  computations  by  AFFDL.  Hence,  a  careful  editing 
and  evaluation  of  all  available  data  to  eliminate  unrepresentative  data 
and/or  data  of  questionable  quality  is  warranted.  This  evaluation  is 
summarized  in  Section  5. 


4.  Homogeneity  Studies  —  This  step  is  needed  to  reduce  the 
redundancy  of  the  available  data  and  the  resulting  prediction  model. 

For  example,  there  is  no  point  in  deriving  separate  prediction  models 
for  the  vibration  environment  of  different  types  of  aircraft  if,  in  fact, 
the  vibration  environments  are  not  significantly  different.  It  follows 
that  careful  investigations  of  the  available  data  for  homogeneity 
(equivalence)  should  be  performed.  Of  particular  interest  are  possible 
differences  in  the  vibration  environment  among  the  three  orthogonal 
axes,  the  various  structural  zones  of  the  aircraft,  the  various  aircraft 
in  a  given  aircraft  group,  and  the  various  aircraft  groups.  These 
homogeneity  tests  are  developed  and  illustrated  in  Section  6. 

5.  Spatial  Distribution  Studies  —  In  some  cases,  it  may  not 
be  feasible  to  include  variables  in  the  prediction  model  which  will 
account  for  differences  in  the  vibration  from  one  location  to  another 
on  the  aircraft  structure.  For  such  cases,  the  predictions  must  be 
based  upon  some  conservative  upper  limit  determined  using  the  distribu¬ 
tion  function  for  the  variation  of  the  vibration  over  the  structure.  The 
determination  of  this  spatial  distribution  function  is  the  subject  of 
Section  7. 


6.  Prediction  Model  Studies  —  Based  upon  the  material 
developed  in  Sections  3  through  7,  the  next  step  is  to  select  the  vari¬ 
ables  for  the  regression  analysis,  and  to  illustrate  the  analysis  pro¬ 
cedures  on  properly  edited  data.  This  is  done  in  Section  8. 
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7.  Test  Level  and  Duration  Selection  Studies  —  Having 
established  an  empirical  procedure  for  predicting  structural  vibration, 
the  next  logical  step  is  to  formulate  specific  techniques  for  converting 
the  resulting  vibration  predictions  into  appropriate  test  specifications. 
This  problem  is  not  straightforward  for  the  case  of  aircraft  vibration 
environments  for  two  reasons.  First,  the  vibration  service  life  of 
many  aircraft  is  measured  in  terms  of  thousands  of  hours.  However, 
laboratory  vibration  tests  must  be  performed  in  much  shorter  periods 
of  time,  preferably  in  terms  of  minutes.  This  means  that  techniques 
must  be  employed  to  compress  the  long  service  life  into  short  test 
durations.  The  second  problem  relates  to  the  basic  nonstationary  chai 
acter  of  the  aircraft  flight  vibration  environment.  The  vibration  leveln 
experienced  by  the  aircraft  during  takeoff  are  considerably  different 
from  those  which  occur  during  cruise.  It  is  desirable,  however,  that 
the  vibration  tests  be  specified  in  terms  of  only  one  vibration  level  (fit 
least  for  a  given  structural  zone).  Hence,  it  is  necessary  to  convert 
the  various  different  vibration  levels  which  occur  during  the  service 
life  to  a  single  vibration  level  for  the  vibration  tests.  These  matters 
are  pursued  in  Section  9. 

8.  Applications  to  Acoustic  Predictions  —  A  final  step  is  to 
investigate  how  the  techniques  for  arriving  at  a  vibration  prediction  pro 
cedure  might  be  modified  to  develop  a  procedure  for  predicting  internal 
acoustic  noise.  This  subject  is  considered  in  Section  10. 


3.  DEVELOPMENT  OF  VIBRATION  RELATIONSHIPS 


3.  1  IDENTIFICATION  OF  VIBRATION  SOURCES 

The  sources  of  aircraft  vibration  may  be  broadly  divi.ded  into  two 
classifications,  random  and  deterministic.  From  various  references 
including  [11-28],  the  principal  random  and  deterministic  sources  of 
vibration  which  might  occur  for  various  types  of  aircraft  are  as  sum¬ 
marized  in  Tables  1  through  4.  All  vibration  sources  taken  together 
are  summarized  in  Table  5.  Note  that  the  aircraft  are  divided  into  five 
types,  where  each  type  is  identified  by  a  "group"  number  consistent 
with  designations  used  by  AFFDL.  Also  note  that  fixed  wing  propeller 
aircraft  (Groups  1  and  2)  are  included  in  the  summary  only  for  com¬ 
pleteness.  The  development  of  prediction  models  for  these  aircraft 
groups  is  probably  not  warranted  since  the  future  design  of  such  air¬ 
craft  is  unlikely. 

From  Table  5,  there  are  eight  principal  sources  of  random  vibra¬ 
tion  and  twelve  principal  sources  of  deterministic  vibration  in  aircraft. 
Each  of  these  sources  will  now  be  briefly  reviewed  in  terms  of  their 
relative  importance  to  the  overall  vibration  environment,  and  their 
suitability  for  consideration  in  the  development  of  a  general  extrapola¬ 
tion  procedure  for  the  prediction  of  aircraft  vibration. 

R-1.  Jet  Acoustic  Noise  —  The  operation  of  turbo-jet  engines 
is  accompanied  by  intense  random  acoustic  noise  caused  by  the  mixing 
of  the  high  velocity  exhaust  gases  with  the  ambient  air.  This  acoustic 
noise  is  a  major  source  of  vibration  in  jet  powered  aircraft  during  take¬ 
off  and  early  climb,  particularly  for  those  aircraft  with  wing  mounted 
engines  (Group  3).  It  becomes  less  significant  as  the  aircraft  speed  in¬ 
creases. 
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•  Air  conditioning  noise 


•  Unstable  aerodynamic  con¬ 
ditions 


Table  2.  Environmental  Sources  for  Turboprop  Transports  (Group  2) 


Source  of 
Random  Vibration 


•  Jet  acoustic  noi 


Source  of 

Deterministic  Vibration 


se 


Propeller  blade  passage 


•  Aerodynamic  boundary 
layer  noise 

•  Propeller  blade  vortex  noise 


•  Propeller  rotation 

•  Turbine  shaft  rotation  and 
blade  passage 


•  Atmospheric  turbulence 
(gust  loads) 

•  Runway  roughness 

•  Wake  turbulence  caused  by 
auxiliary  lift  devices, 
drag  devices,  and  spoilers 

•  Air  conditioning  noise 


•  Engine  accessory  equipment 

•  Auxiliary  rotating  equipment 

•  Gear  box  noise 

•  Acoustic  cavity  resonances 

•  Unstable  aerodynamic  con¬ 
ditions 


Table  3.  Environmental  Sources  for  Jet  Bombers  (Group  3)  and 
Century  Jet  Fighters  (Group  5) 


Sources  of 

Deterministic  Vibration 


•  Jet  acoustic  noise  •  Turbine  shaft  rotation  and 

blade  passage 

•  Aerodynamic  boundary 

layer  noise  •  Engine  accessory  equipment 

•  Atmospheric  turbulence  •  Auxiliary  rotating  equipment 

(gust  loads) 

•  Gunfire 

•  Runway  roughness 

•  Acoustic  cavity  resonances 

•  Wake  turbulence  caused  by 

avixiliary  lift  devices,  •  Unstable  aerodynamic  con- 

drag  devices,  and  spoilers  ditions 

•  Transonic  shock  wave - 
boundary  layer  interaction 


Sources  of 
Random  Vibration 


•  Airconditioning  noise 


n 

n 

n 


Table  4.  Environmental  Sources  for  Helicopters  (Group  10) 


n 

n 

n 

n 

n 

n 

n 

n 

0 

D 

D 

0 

0 

0 

I 


Source  of 
Random  Vibration 


Source  of 

Deterministic  Vibration 


Jet  acoustic  noise 
Rotor  blade  vortex  noise 
Atmospheric  turbulence 
Airconditioning  noise 


Rotor  blade  passage  (main 
and  tail) 

Rotor  rotation 


Jet  engine  turbine  shaft  rota¬ 
tion  and  blade  passage’^ 


•  Reciprocating  engine  exhaust 


•  Reciprocating  engine  shaft 
rotation** 


Engine  accessory  equipment 
Auxiliary  rotating  equipment 


•  Gear  box  noise 

•  Rotor  blade  slope 

•  Gunfire 


* 

Jet  engine  powered  helicopters  only 

jjcjjt 

Reciprocating  engine  powered  helicopters  only 
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Table  5.  Summary  of  Environmental  Sources  of  Vibration 


R-2.  Aerodynamic  Boundary  Layer  Noise  —  The  flight  of  air¬ 
craft  at  high  speeds  produces  a  turbulent  boundary  layer  at  the  inteixace 
between  the  aircraft  and  the  atmosphere.  This  turbulence  produces 
vibration  which  increases  in  severity  as  dynamic  pressure  increases. 
For  high  speed  aircraft  (Groups  3  and  5),  it  is  a  major  source  of  vibra¬ 
tion  during  cruise  and  other  high  speed  flight  conditions. 

R-3.  Propeller  (or  Rotor)  Blade  Vortex  Noise  —  Fixed  wing 
aircraft  propellers  and  helicopter  rotors  are  responsible  for  two  types 
of  acoustic  excitation.  The  first  is  a  periodic  excitation  related  to 
propeller  or  rotor  blade  passage,  which  is  discussed  later.  The  second 
is  a  random  excitation  produced  by  vortex  shedding  from  the  propeller 
or  rotor  blades.  Vortex  noise  may  be  quite  pronounced  when  propellers 
or  rotors  are  operated  at  high  rpm  with  low  pitch  (high  rpm  but  low 
power).  For  the  more  interesting  case  where  flight  power  is  delivered 
to  the  propellers  and  rotors,  the  vortex  noise  is  generally  not  significant 
compared  to  the  vibration  induced  by  blade  passage  and  other  sources. 

R-4.  Atmospheric  Turbulence  —  Atmospheric  turbulence  is  a 
principal  source  of  low  frequency  vibration  in  aircraft.  However,  it  is 
not  a  continuous  source  of  vibration.  Its  occurrence  is  a  function  only 
of  atmospheric  conditions  and  not  of  aircraft  flight  parameters.  Hence, 
there  is  no  way  to  generally  predict  the  vibration  induced  by  atmospheric 
turbulence  based  upon  parameters  of  the  aircraft  structure,  engine  oper¬ 
ating  conditions,  and/or  aircraft  flight  conditions.  Of  course,  given  the 
knowledge  that  atmospheric  turbulence  will  occur,  one  can  calculate  the 
aircraft  response  to  the  turbulence,  as  outlined  in  [12]  .  Predictions 
of  this  type  should  be  handled  separately  and  should  not  be  included  in  a 
general  prediction  procedure  of  the  type  being  studied  here. 


R-5.  Runway  Roughneis  —  Runway  roughness  is  a  source  of 
low  frequency  vibration  during  takeoff,  landing  and  taxi.  The  resulting 
vibration  is  a  function  of  the  specific  runway  condition  and  aircraft. 

Hence,  like  atmospheric  turbulence,  the  runway  roughness  problem 
should  be  handled  separately  and  should  not  be  included  in  a  general 
prediction  procedure. 

R-6.  Wake  Turbulence  —  Wake  turbulence  caused  by  auxiliary 
lift  devices,  drag  devices  and  spoilers  can  be  a  significant  source  of 
aircraft  vibration,  particularly  during  landing  approaches.  It  should  be 
considered  in  a  vibration  prediction  model  for  any  type  of  aircraft  for 
those  flight  co  ■  Utions  where  such  devices  are  used. 

R-7.  Transonic  Vibration  —  For  missiles  and  spacecraft,  the 
interaction  of  shock  waves  and  boundary  layer  pressures  during  transonic 
flight  is  often  the  predominate  source  of  flight  vibration.  The  relatively 
clean  geometry  of  modern  aircraft,  however,  greatly  restricts  the  influ¬ 
ence  of  transonic  shock  wave-boundary  layer  interaction.  Any  sig¬ 
nificant  vibration  induced  by  transonic  flight  would  probably  be  very 
localized  and  a  function  of  the  detailed  geometry  of  the  aircraft. 

Hence,  such  vibration  should  not  be  included  in  a  general  prediction 
procedure. 

R-8.  Airconditioning  Noise  —  Airconditioning  noise  is  usually 
more  of  an  acoustical  comfoi  roblem  than  a  structural  vibration  prob¬ 
lem.  The  vibration  produced  by  airconditioning  noise  is  generally  local 
and  high  frequency  in  nature.  It  is  rarely  a  significant  contributor  to  the 
overall  vibration  environment  of  aircraft,  except  perhaps  in  very  re¬ 
stricted  structural  areas.  It  is,  however,  a  major  source  of  cabin  acovistic 


D-1.  Propeller  (or  Rotor)  Blade  Passage  —  As  discussed 
under  R-3,  fixed  wing  aircraft  propellers  and  helicopter  rotors  are 
responsible  for  two  types  of  acoustic  excitation.  One  type  is  due  to 
vortex  shedding  from  the  propeller  or  rotor  blades,  as  discussed  pre¬ 
viously.  The  other  type  is  caused  by  the  periodic  pressure  pulses 
emitted  from  the  propeller  or  rotor  blades.  This  second  type  of  ex¬ 
citation  is  most  severe  for  the  case  of  fixed  wing  aircraft  with  wing 
mounted  engines,  where  the  primary  aircraft  structure  is  outside  the 
cylinder  of  blade  rotation.  For  this  type  of  aircraft  (Groups  1  and  2), 
propeller  blade  passage  is  the  major  source  of  vibration.  For  helicopters 
(Group  10),  rotor  blade  passage  is  a  less  severe  source  relative  to  other 
vibration  sources,  but  still  may  be  significant. 

D-2.  Propeller  (or  Rotor)  Rotation  —  The  vibration  induced  by 
propeller  rotation  on  fixed  wing  aircraft,  or  rotor  rotation  on  helicopters, 
is  a  direct  function  of  the  magnitude  of  propeller  or  rotor  unbalance  and 
RPM.  For  fixed  wing  aircraft,  propeller  rotation  is  a  negligible  source 
o!  vibration  under  normal  conditions,  and  can  be  ignored  in  the  prediction 
procedure  being  studied  here.  For  helicopters,  however,  the  main 
rotor  rotation  may  be  a  significant  source  of  low  frequency  vibration. 

D-3.  Jet  Engine  Shaft  Rotation  and  Blade  Passage  —  Periodic 
vibration  due  to  jet  engine  shaft  unbalance  may  be  quite  significant  in  jet 
powered  aircraft,  particularly  those  with  fuselage  incorporated  engines 
(Group  5).  The  vibration  occurs  at  the  rotational  frequency  of  the  shaft 
(or  shafts)  and  all  harmonics  thereof.  Some  vibration  may  also  occur 
at  frequencies  corresponding  to  the  rate  at  which  turbine  and  compressor 
blades  pass  the  engine  stator  blades.  This  later  excitation,  however,  is 
usually  more  of  an  acoustical  comfort  problem  than  a  structural  vibration 
problem. 
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D-4. 


engine  type  aircraft  (Groups  1  and  10),  the  engine  ignition  contributes 
to  the  aircraft  vibration  environment  in  two  ways.  First,  the  actual 
ignition  explosions  produce  engine  vibration  which  transmits  in  part 
through  the  engine  mounts  into  the  airframe.  Second,  the  acoustical 
pressure  pulses  from  the  exhaust  induce  vibration  much  like  the  pro¬ 
peller  blade  passage  does.  The  first  source  is  usually  negligible 
under  normal  conditions,  but  the  second  source  can  be  significant. 

D-5,  D-6,  D-7,  D-8.  Reciprocating  Engine  Shaft  Rotation, 

Engine  Accessory  Equipment,  Auxiliary  Rotating  Equipment,  Gear  Box 
Noise  —  These  various  periodic  contributions  are  usually  negligible 
sources  of  vibration  in  fixed  wing  aircraft,  but  may  be  significant 
sources  in  helicopters.  In  particular,  gear  box  noise  is  often  a  major 
source  of  vibration  in  helicopters. 

D-9.  Rotor  Blade  Slap  When  helicopters  perform  sharp 
turn  maneuvers  or  approach  limiting  forward  speeds,  the  main  roto^i 
blades  moving  aft  (in  the  direction  of  air  flow)  may  experience  a  partial 
stall  condition.  This  repetitive  quasi- stall  of  the  blades  induces  a 
momentary  flow  separation  called  blade  slap.  This  condition  can  be 
a  significant  source  of  low  frequency  vibration,  as  well  as  acoustic 
noise.  Since  it  occurs  only  under  special  conditions,  however,  it  should 
be  handled  separately  and  should  not  be  included  in  a  general  prediction 
procedure. 


D-10.  Gunfire  — -  Gunfire  can  produce  considerable  vibration, 
particularly  in  the  structural  regions  near  the  gun  locations.  It  should 
be  considered  in  the  prediction  model  for  those  aircraft  which  are 
equipped  with  guns. 


D-11,  D-12.  Acoustic  Cavity  Resonances  and  Unstable  Aero¬ 
dynamic  Conditions  —  These  two  sources  of  periodic  vibration  will 
hopefully  not  be  present  since  they  imply  a  poor  design.  If  they  are 
present,  they  constitute  a  special  problem  which  must  be  dealt  with  by 
special  techniques.  The  inclusion  of  such  techniques  in  a  general  pre¬ 
diction  model  is  clearly  not  feasible. 

In  summary,  the  various  sources  of  random  and  deterministic 
vibration  may  be  classified  into  four  categories  as  follows. 

A.  Sources  which  are  considered  sufficiently  significant 
to  include  in  the  prediction  model 

B.  Sources  which  are  considered  significant,  but  which 
should  be  treated  separately  and  not  as  part  of  the 
prediction  model 

C.  Sources  which  are  not  considered  sufficiently  significant 
to  include  in  the  prediction  model 

D.  Sources  which  do  not  apply  to  the  aircraft  group  in 
question . 

The  classification  of  sources  for  each  aircraft  group  is  summarized  in 
Table  6.  Note  that  these  classifications  are  tentative  and  may  be  altered 
by  later  data  studies. 
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Table  6.  Classification  of  Sources  for  Each  Aircraft  Group 


—T 

I  H 


3.  2  ANALYTICAL  DESCRIPTION 

OF  STRUCTURAL  VIBRATION 

All  of  the  significant  random  sources  of  aircraft  vibration  sum¬ 
marized  in  Table  6  are  pressure  field  type  sources.  Hence,  as  a 
preliminary  step  to  establishing  relationships  between  these  sources 
and  the  resulting  structural  vibration,  a  review  of  the  theory  for  struc¬ 
tural  response  to  random  pressure  excitations  is  in  order. 

Consider  an  arbitrary  structure  which  is  subjected  to  a  random 
pressure  field,  p(x,  t)  .  The  pressure  field  may  be  described  by  a 
spatial  cross  spectral  density  function,  as  follows. 


G  (x,  x',  f)  =  E 
P 


lim  Y  P*  (x,f)  P^(x',f) 
T-*oo 


where 


f) 


r  p(^,  t)  e  dt 

-oo 


(11) 


X,  x'  =  two  different  vector  points  on  the  structure 

* 

P^(x,  f)  =  complex  conjugate  of  P^(x,  f) 


Note  that  the  spatial  cross  spectrum  can  also  be  defined  as 


Gp(x,  x',  f) 


f  R  (x,  x',t)  e  dr 

•^-oo  P 


(12) 
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where 


R  (x.x',  t) 
P 


lim  “T  /  p(x,  t)  p(x',  t  +  t)  dt 
T—oo  ^  -'-T 


From  [13]  ,  the  power  spectrum  for  the  acceleration  response  of  the 
structure  at  any  point  x  induced  by  this  pressure  field  is  given  by 


G  (x.f)  = 
a 


.4 


i  k  M.  fj  f, 


Z.2  ‘ 


(13) 


where 


=/7^i 
0  0  ‘ 


(x)  <}>^(x')  G  (x,  x',f)  dx  dx'  (modal  cross 
P  spectrum) 


M.  =  /  m(x)  4>.  (x)  dx  (modal  mass) 

0  ' 


(j>.(x)  =  mode  shape  for  ith  normal  mode 
Hj(f)  =  frequency  response  function  for  ith  normal  mode 
f^  =  undamped  natural  frequency  for  ith  normal  mode 
m(x)  =  surface  mass  density 
A  =  surface  area 


The  modal  cross  spectrum  in  Eq.  (13)  can  be  reduced  ti  a  more  con¬ 
venient  form,  as  follows. 


1  -  (f/f,) 


+  24,  f/f, 


The  power  spectrum  values  of  principal  concern  are  those  at  frequencies 
near  normal  mode  frequencies,  since  these  frequencies  will  be  the  peak 
values  in  the  spectrum.  Letting  f  =  L  ,  Eq.  (22)  becomes 


1 


(23) 


where  Qj  =  1  /(2t.)  • 

The  quantity  (k/i)  l<«>l  is  plotted  against  the  frequency  ratio 
f  /f,  for  various  values  of  damping  in  Figure  2.  It  is  clear  from  Figure  2 
that,  on  the  average,  the  contribution  of  the  i  ^  k  terms  in  Eq.  (13) 
falls  off  very  rapidly  with  normal  mode  frequency  ratio  for  small  damp¬ 
ing.  For  example,  if  0.  01  (Q  =  50) ,  the  contribution  of  the  kth 
mode  to  the  power  spectrum  at  the  ith  normal  mode  frequency  is  less 
than  10%  when  the  kth  mode  is  only  10%  from  the  ith  mode  frequency. 
Since  the  damping  ratio  for  most  aircraft  structures  is  relatively  small, 
it  appears  justified  as  a  first  order  of  approximation  to  ignore  the  i  ^  k 
terms  in  Eq.  (13).  With  this  approximation,  the  power  spectrum  for  the 
acceleration  response  of  a  structure  may  be  written  as 
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(24) 


G^f) 


M.  f. 

i  1 


Gj(f) 


By  further  assuming  that  the  generalized  masses  and  joint  acceptances 
are  constants  =  (w/g)  A  and  (f)  =  j  ,  Eq.  (24)  can  be 

further  reduced  to 


(25) 


2  2  r  i2 

On  the  average,  the  quantities  A  and  I  <t>j(x)  j  should  not  vary 
sharply  for'^similar  types  of  structures.  Hence,  as  a  first  order  of 
approximation,  it  appears  reasonable  to  assume  that  the  power  j 
spectrum  for  the  acceleration  response  of  a  structure  is  proportional 
to  the  average  power  spectrum  of  the  excitation  over  the  structure 
divided  by  the  square  of  the  surface  weight  density  for  the  structure. 
That  is, 


G  (x,f)r- 
a 


2 


(26) 
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3.  3  ANALYTICAL  DESCRIPTION 
OF  PRIMARY  SOURCES 


From  Table  6,  there  are  three  random  and  seven  deteri?iinistic 
sources  of  vibration  which  are  considered  sufficiently  significant  and 
otherwise  suitable  to  include  in  the  vibration  prediction  model.  These 
significant  sotiroes  are: 

R-1  Jet  acoustic  noise 

* 

R-2  Aerodynamic  boundary  layer  noise 

* 

R-6  Wake  turbulence  caused  by  auxiliary  lift  devices, 
drag  devices,  and  spoilers 

D-1  Propeller  (or  rotor)  blade  passage 

D-2  Helicopter  rotor  rotation 

D-3  Jet  engine  turbine  shaft  rotation 

D-4  Reciprocating  engine  exhaust 

D-6  Engine  accessory  equipment 

D-8  Gear  box  noise 

D-10  Gunfire 

Five  of  the  ten  sources  listed  above  (those  marked  with  an  asterisk) 
are  basically  pressure  field  type  sources;  that  is,  they  induce  vibration 
by  generating  fluctuating  pressure  fields  which  impinge  over  the  struc¬ 
tural  surface.  In  Section  3.  2,  it  was  established  that,  as  a  first  order  of 
approximation,  the  resulting  structural  vibration  may  be  considered  pro¬ 
portional  to  the  power  spectrum  of  the  generated  pressure  field.  Hence, 
for  the  purposes  of  a  vibration  prediction  model,  it  is  desirable  to  define 
the  power  spectrum  for  these  seven  sources  in  terms  of  engine  operating 
conditions  and/or  aircraft  flight  conditions,  where  possible. 


The  remaining  five  sources  (D-2,  D-3,  D~6,  D-8,  and  D-10)  are 
basically  mechanical  type  sources  which  induce  vibration  by  direct 
mechanical  transmission  of  energy  through  the  structure.  The  defini¬ 
tion  of  vibration  as  a  function  of  operational  parameters  is  more 
difficult  for  these  sources  since  the  mechanical  impedance  of  the  struc 
ture  is  the  primary  factor  which  controls  the  resulting  vibration  in¬ 
tensity.  The  use  of  detailed  impedance  information  is  not  appropriate 
for  the  type  of  general  vibration  model  of  interest  here.  Hence,  it 
will  be  necessary  to  use  more  general  statistic  techniques  to  describe 
these  mechanical  sources. 

Each  of  the  ten  significant  sources  of  aircraft  vibration  will  now 
be  discussed  in  terms  of  their  relationships  to  various  descriptive 
parameters  of  engine  operating  conditions  and  flight  conditions. 

3.  3.  1  Jet  Acoustic  Noise 

From  [14] ,  the  total  acoustic  power  radiated  by  a  jet  engine  at 
rest  is  related  to  various  engine  and  atmospheric  parameters  by 


P/-pAV®/c^  (27) 

where 

p  =  air  density 
A  =  jet  exit  area 

=  exhaust  gas  velocity 
c  =  speed  of  sound 
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POWER  SPECTRUM  LEVEL 


Furthermore,  the  power  spectrum  for  the  noise  is  approximately  as 
shown  in  Figure  3. 


Figure  3,  Spectrum  for  Jet  Acoustic  Noise 


In  Figure  3,  f  =  frequency  in  cps,  d  =  jet  exit  diameter,  and  c  =  local 

e 

speed  of  sound.  Note  that  for  typical  jet  engines  at  takeoff  conditions, 
the  spectral  peak  usually  occurs  at  a  frequency  between  100  and  300  cps, 
but  because  of  the  shallow  slope  of  the  spectrum,  significant  acoustical 
power  exists  over  the  entire  audio  frequency  range. 
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The  power  spectrum  for  the  acoustic  pressure  at  any  point  on  the 
aircraft  structure  due  to  jet  noise  is  a  function  of  the  directivity  pattern 
of  the  noise  as  well  as  the  structural  location  relative  to  the  jet  engine, 
the  aircraft  velocity,  wind  conditions,  and  the  presence  of  reflecting 
surfaces.  For  the  problem  at  hand,  it  will  be  necessary  to  ignore  vari¬ 
ations  in  the  directivity  pattern,  wind  conditions,  and  reflecting  sur¬ 
faces;  that  is,  these  factors  will  be  permitted  to  appear  as  variability 
in  the  resulting  prediction  model.  Structural  location  relative  to  the 
jet  engine  and  aircraft  velocity,  however,  cannot  be  ignored.  The  vari¬ 
ation  of  sound  pressure  level  with  structural  location  can  be  dealt  with 
by  deriving  a  different  set  of  A,  .(f)  weights  for  each  of  several  different 
structural  regions  or  zones.  The  variation  of  sound  pressure  level  with 
airspeed  could  be  dealt  with  in  a  similar  way,  but  it  appears  more 
appropriate  to  simply  restrict  attention  to  a  single  airspeed;  namely 
zero.  This  should  be  adequate  for  the  desired  prediction  model  since 
the  jet  noise  impinging  on  the  structure  is  greatest  when  the  engine  is 
at  maximum  thrust,  and  maximum  thrust  occurs  during  the  takeoff  roll 
when  the  airspeed  is  near  zero.  Furthermore,  as  the  airspeed  increases, 
the  directivity  pattern  of  the  jet  noise  shifts  aft,  rapidly  reducing  the 
contribution  of  jet  noise  to  the  vibration  environment. 

In  conclusion,  for  a  given  structural  zone  of  an  aircraft  at  zero 
airspeed,  the  power  spectrum  for  the  acoustic  noise  impinging  on  the 
structure  is  related,  as  a  first  order  of  approximation,  to  engine  and 
atmospheric  parameters  by  ' 


0  * 


G(n)  r-  pAV^'/c® 


\ 


where  is  dimensionless  frequency  (fd  /  V^c) ,  and  n  is  some 
number  probably  between  6  and  8  .  A  suitable  value  for  n  should  be 
determined  by  optimizing  the  correlation  coefficient  for  the  vibration 
model  during  the  regre’ssion  studies. 

3.  3.  2  Aerodynamic  Boundary  Layer  Noise 

From  [16] ,  for  Mach  numbers  less  than  2.  5,  the  mean  square 
value  of  the  pressure  generated  by  boundary  layer  turbulence  for  rela¬ 
tively  "clean"  structures  is  approximated  by 


2 

(T 

P 


jsi  (0.  006  q) 


2 


(29) 


where  q  is  the  free  stream  dynamic  pressure  (q  =  p  V  /  2)  .  Further¬ 
more,  the  power  spectrum  for  the  noise  is  approximately  as  shown  in 
Figure  4. 


Figure  4.  Spectrum  for  Aerodynamic  Boundary  Layer  Noise 
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In  Figure  4,  6*  is  the  boundary  layer  thickness,  and  is  given  by 
6*  fs  0.  0026  X  for  "clean"  structures  where  x  is  the  distance  from 
the  front  of  the  structure. 

The  power  spectrum  for  the  boundary  layer  pressures,  as  shown 
in  Figure  4,  is  relatively  uniform  for  dimensionless  frequencies  up  to 
f6*  /  V  fas  0.  2  .  Beyond  this  frequency,  the  spectrum  falls  off  gradually 
until  f6*  /  V  ^  2.  0  ,  where  the  spectrum  cuts  off  sharply.  For  a  50  ft 
long  airplane  traveling  at  1000  ft/sec  (about  700  mph),  this  means  that 
the  spectrum  starts  falling  off  at  about  3000  cps  for  a  point  halfway  back 
on  the  airplane,  and  at  about  1500  cps  for  a  point  at  the  rear  of  the 
airplane.  In  practice,  it  can  probably  be  assumed  that  the  power 
sp^trum  of  the  boundary  layer  pressures  is  white  noise  in  the  fre¬ 
quency  range  of  usual  interest  (0  to  2000  cps).  Hence,  as  a  first 
order  of  approximation,  the  power  spectrum  for  the  boundary  layer 
noise  pressure  on  the  structure  is  related  to  flight  conditions  by 


G(i)r-  q^ 


(30) 


3.  3.  3  Wake  Turbulence 

Wake  turbulence  caused  by  auxiliary  lift  devices,  drag  devices, 
and  spoilers  is  no  doubt  a  function  of  aircraft  flight  parameters  such 
as  airspeed  and  density.  However,  the  actual  geometry  of  the  devices 
involved  is  a  more  significant  influencing  factor.  Furthermore,  the 
problem  generally  arises  for  a  limited  range  of  flight  parameters, 
since  such  devices  are  used  only  for  special  conditions  such  as  landing 
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and  takeoff.  This  means  that  an  effort  to  correlate  wake  turbulence 
induced  vibration  and  aircraft  flight  conditions  may  not  be  very  profitable. 
Perhaps  a  better  approach  is  to  consider  the  wake  turbulence  induced 
vibration  as  a  separate  environment  which  will  be  predicted  using  a  gen¬ 
eral  statistical  evaluation  of  all  data  taken  in  each  of  various  structural 
zones  when  the  various  turbulence  inducing  devices  are  deployed.  For 
the  case  of  speed  break  deployment,  the  possibility  of  a  relationship 
between  vibration  and  dynamic  pressure  might  be  investigated.  Fortunately, 
such  studies  are  easy  to  accomplish  because  AFFDL  identifies  and  sep¬ 
arates  all  data  taken  with  and  without  the  deployment  of  such  devices, 

3.  3.  4  Propeller  (or  Rotor)  Blade  Passage 

Consider  first  the  case  of  fixed  wing  aircraft  with  wing  mounted 
propeller  engines  (Groups  1  and  2).  From  [12],  the  overall  sound  power 
level  generated  by  a  propeller  in  the  near  field  is  proportional  to 


,  3. 66M 

P'^IO  ^ 


3.  34+  2.  44M 

V  D  (z/D)  ^R^/528) 


(31) 


where 

P^  =  mechanical  power  to  propeller 

=  propeller  tip  Mach  number 

V  =  number  of  propeller  blades 

z  =  clearance  between  propeller  blade  tip 
and  fuselage 

D  =  propeller  diameter 


=  absolute  temperature  in  degrees  Rankine 
(degrees  F  +  460) 


The  power  spectrum  for  the  propeller  noise  (neglecting  vortex  noise) 
will  be  a  series  of  harmonic  components  with  frequencies  of 


iv  (RPM) 
60 


1—  If  3f 


(32) 


For  the  propeller  aircraft  of  concern  to  AFFDL,  the  fundaniei^tal  fre¬ 
quency,  f^  ,  is  approximately  50  cps  for  normal  cruise  conditions. 

The  amount  of  power  in  the  various  harmonics  is  a  function  of  . 

For  <  0.  8,  over  50%  of  the  power  occurs  at  the  fundamental  fre¬ 
quency  (f  ) .  As  M  “*1.0,  the  power  in  the  harmonics  increases 
relative  to  the  power  of  the  fundamental.  | 

As  for  jet  acoustic  noise,  the  power  spectrum  for  the  propeller 
acoustic  pressures  at  any  point  on  the  aircraft  structure  is  a  function 
of  the  directivity  pattern  of  the  blade  passage  noise  as  well  as  the  struc¬ 
tural  location  relative  to  the  propeller,  the  aircraft  velocity,  wind  con¬ 
ditions,  and  the  presence  of  reflecting  surfaces.  Hence,  it  will  be 
necessary  to  deal  with  the  propeller  noise  problem  in  the  same  manner 
as  suggested  for  the  jet  noise  problem  in  Section  3.  3.  1.  Specifically, 
a  prediction  model  for  propeller  induced  vibration  should  be  developed 
by  deriving  a  different  set  of  weights  for  each  of  several  different 

structural  zones  at  zero  airspeed  only.  If  this  is  done,  the  power  spec¬ 
trum  for  the  propeller  noise  impinging  on  the  structure  will  be  a  series 
of  harmonically  relate''  delta  functions  at  the  frequencies  given  in  Eq.  (32), 
and  with  a  total  mean  square  value  proportional  to  P  in  Eq.  (31) . 
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Now  consider  the  case  of  helicopters.  There  is  a  similar  problem 
of  vibration  induced  by  the  tail  rotor  blade  passage,  and  Eq.  f31)  applies 
to  this  problem.  For  the  case  of  the  main  rotor,  however,  the  problem 
is  somewhat  different.  Equation  (31)  applies  only  to  structure  outside 
the  cylinder  of  rotation  for  the  propeller  or  rotor.  Most  of  the  struc¬ 
ture  for  helicopters  is  inside  the  cylinder  of  rotation  for  the  main  rotor. 
A  model  for  this  vibration  is  not  known  to  be  available.  Hence,  it  will 
be  necessary  to  establish  a  prediction  model  for  the  helicopter  vibration 
induced  by  main  rotor  blade  passage  by  careful  inspection  of  past  data. 

3.  3.  5  Helicopter  Rotor  Rotation 

The  vibration  induced  by  rotor  rotation  on  helicopters  is  a  direct 
function  of  the  magnitude  of  the  rotor  unbalance  and  rotor  RPM,  as  well 
as  the  mechanical  impedance  of  the  supporting  structure.  The  frequency 
of  the  rotor  rotation  induced  vibration  is  generally  below  10  cps.  It  is 
unlikely  that  a  strong  correlation  would  be  obtained  with  flight  conditions 
and  engine  operating  conditions,  other  than  RPM.  Since  the  RPM  for 
helicopter  rotors  is  relatively  constant  for  most  flight  conditions  and 
rotor  unbalance  is  usually  unknown,  it  is  suggested  that  rotor  rotation 
induced  vibration  be  predicted  using  a  general  statistical  evaluation  of 
all  data  taken  in  each  of  various  structural  zones. 

3.  3.  6  Jet  Engine  Turbine  Shaft  Rotation 

The  problem  of  vibration  induced  by  the  jet  engine  turbine  shaft 
rotation  is  similar  to  the  problem  of  helicopter  rotor  vibration  discussed 
in  Section  3.  3.  5,  and  should  be  handled  in  a  similar  way.  Specifically, 
shaft  rotation  induced  vibration  should  be  predicted  using  a  general 
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statistical  evaluation  of  all  data  taken  in  each  of  various  structural 
zones.  The  fundamental  frequency  of  shaft  rotation  induced  vibration 
is  generally  between  50  cps  and  150  cps  for  most  jet  engine  operating 
conditions. 

3.  3.  7  Reciprocating  Engine  Exhaust 

From  [12] ,  the  overall  sound  power  level  generated  by  the  ex¬ 
haust  of  a  reciprocating  engine  is  proportionaLto  the  total  horsepower 
delivered  by  the  engine.  Furthermore,  the  spectrum  for  the  exhaust 
noise  is  basically  periodic  with  a  fundamental  frequency  equal  to  the 
ignition  firing  frequency  (from  100  cps  to  200  cps  for  most  engines 
and  operating  conditions).  The  power  content  of  the  harmonics  falls 
off  at  a  rate  of  about  3:1  per  octave.  The  vibration  induced  by  the  ex¬ 
haust  is  a  function  of  the  directivity  pattern,  structural  location,  etc.  , 
as  discussed  for  jet  noise  in  Section  3.  3. 1.  Hence,  it  will  be  necessary 
to  derive  a  difi.'er»?nt  set  of  weights  for  each  of  several  different 

structural  zones  for  various  different  flight  conditions. 

3.  3.  8  Engine  Accessory  Equipment 

The  vibration  induced  by  engine  accessory  equipment  tends  to  be 
localized  and  relatively  independent  of  flight  conditions.  It  may  occur 
at  frequencies  up  to  1000  cps.  The  most  important  factors  influencing 
the  magnitude  of  the  vibration  are  the  details  of  the  equipment  and  the 

I 

mounting  point  impedance  of  the  supporting  structure.  It  is  suggested 
that  accessory  equipment  induced  vibration  be  predicted  using  a  general  | 
statistical  evaluation  of  all  data  collected  in  each  of  several  strl^ctural 


zones. 


The  vibration  generated  by  helicopter  gear  box  noise  is  quite 
complex  and  may  occur  ever  a  wide  frequency  range.  Although  the 
magnitude  of  the  vibration  is  unquestionably  related  to  engine  operating 
conditions  and  perhaps  flight  conditions,  the  relationships  are  not  well 
defined.  It  is  suggested  that  gear  box  induced  vibration  be  predicted 
using  a  general  statistical  evaluation  of  all  data  collected  at  the  fre¬ 
quencies  corresponding  to  anticipated  gear  box  frequenciej. 

3.  3.  10  Gunfire 

The  vibration  induced  by  gunfire  is  clearly  not  related  to  engine 
operating  conditions  or  flight  conditions.  Hence,  the  contribution  of 
gunfire  to  the  vibration  environment  should  be  predicted  based  upon  a 
general  statistical  evaluation  of  data  collected  when  guns  are  fired. 


4.  REVIEW  OF  AVAILABLE  DATA 

The  flight  acoustic  and  vibration  data  available  to  AFFDL,  and 
the  acquisition  and  analysis  procedures  used  to  obtain  the  data  are  sum- 
marized  in  this  section. 

4.  1  SUMMARY  OF  AIRCRAFT 

As  of  1  January  1968,  AFFDL  has  documented  flight  vibration 
data  collected  on  a  total  of  13  different  aircraft.  Flight  acoustic  data 
were  also  collected  and  documented  for  two  of  those  aircraft.  The 
13  aircraft  are  classified  into  5  different  groups  depending  upon  their 
basic  mission  and  configuration.  In  Table  7,  the  specific  aircraft  in 
each  of  the  5  groups  are  presented  along  with  the  AFFDL  aircraft  code 
number  and  report  number  which  summarizes  the  results  of  the  acoustic 
and  vibration  surveys  for  each  of  the  aircraft. 

For  the  purposes  of  test  specifications  and  design  criteria,  AFFDL 
considers  an  aircraft  to  be  composed  of  a  possible  27  different  structural 
zones.  Each  of  these  structural  zones  and  the  code  number  used  to 
identify  the  zone  are  presented  in  Table  8.  The  data  measured  in  each 
zone  of  each  aircraft  are  further  classified  in  terms  of  the  aircraft  flight 
condition  at  the  time  of  measurement.  Each  flight  condition  and  the  code 
number  used  for  its  identification  are  summarized  in  Table  9. 

Using  the  code  numbers  in  Tables  7  through  9,  the  general  location 
and  flight  condition  for  any  vibration  measurement  can  be  readily  identified. 
For  example,  a  measurement  in  zone  03  of  aircraft  46  for  flight  condition 
07  would  be  a  measurement  in  the  aft  quarter  of  the  fuselage  of  an  F-102A 
during  normal  cruise.  Note  that  the  exact  location  for  each  vibration 
measurement  is  detailed  in  the  various  documentation  reports. 
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Table  8,  Summary  of  Structural  Zones 


Description  of  Structural  Zone 

Code  Number 

Forward  quarter  of  fuselage 

01 

Center  half  of  fuselage 

02 

Aft  quarter  of  fuselage 

03 

Vertical  and  horizontal  stabilizer,  including 
rudder  and  elevators 

04 

Outer  one -third  of  wing 

05 

Inner  two-thirds  of  wing 

06 

Engine 

07 

Rigidly  mounted  equipment  in  zone  01 

08 

Rigidly  mounted  equipment  in  zone  02 

09 

Rigidly  mounted  equipment  in  zone  03 

10 

Rigidly  mounted  equipment  in  zone  04 

11 

Rigidly  mounted  equipment  in  zone  05 

12 

Rigidly  mounted  equipment  in  zone  06 

1.3 

Rigidly  mounted  equipment  in  zone  07 

14 

Shock  mounted  equipment  in  zone  01 

15 

Shock  mounted  equipment  in  zone  02 

16 

Shock  mounted  equipment  in  zone  03 

17 

Shock  mounted  equipment  in  zone  04 

18 

Shock  mounted  equipment  in  zone  05 

19 

Shock  mounted  equipment  in  zone  06 

cO 

Shock  mounted  equipment  in  zone  07 

21 

Accessory  section 

22 

Helicopter  section 

23 

Rigidly  mounted  equipment  in  zone  22 

24 

Rigidly  mounted  equipment  in  zone  23 

25 

Shock  mounted  equipment  in  zone  22 

26 

Shock  mounted  equipment  in  zone  23 

27 
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Table  9.  Summary  of  Flight  Conditions 


Description  of  Flight  Conditions 

Code 

Number 

Tajci 

01 

Ground  Runup  (clean) 

02 

Takeoff  Roll 

03 

Takeoff 

04 

Takeoff  (assisted  i.  e.  A/B,  jato,  etc. ) 

05 

Climb  (normal) 

06 

Cruise  (normal) 

07 

Cruise  (speed  brakes  extended) 

08 

Cruise  (rocket  pod  extended,  missile  launcher  ex¬ 
tended,  gunfire  or  test  item  operating) 

09 

Cruise  (09  +  speed  brakes) 

10 

Vertical  Dive  (clean) 

11 

Vertical  Dive  (w/09) 

12 

Vertical  Dive  (w/speed  brake) 

13 

Vertical  Dive  (w/10) 

14 

Cruise  (flaps  extended) 

15 

Cruise  (gear  extended) 

16 

Cruise  (refueling  doors  open) 

17 

Inverted  Flight 

18 

Standard  Pursuit  Curve  (clean) 

19 

Standard  Pursuit  Curve  (w/gunfire) 

20 

Normal  Descent  (clean) 

21 

Normal  Descent  (speed  brakes  extended) 

22 

Normal  Approach 

23 

Normal  Approach  (w/flaps,  gear,  etc.) 

24 

1 
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Table  9  (continued) 


Description  of  Flight  Conditions 


Touchdown 
Landing  Roll 

Drag  Chute  or  Reverse  Thrust  (or  both) 

Bomb  Bay  Open,  Camera  Doors  Open  or 
Troop  Doors  Open 

Cruise  (turbulent  air) 

Cruise  (one  or  more  engines  out) 

Turn  (1.5g  or  less) 

Turn  (more  than  1.  5g) 

Auto -Rotation  —  Engine  at  Idle 
%  RPM  =  Rotor  RPM  (helicopter) 

Auto -Rotation  Stop  —  as  above  (helicopter) 

Hover  (helicopter) 

Rearward  Flight  (helicopter) 

Side  Flight  (helicopter) 

Emergency  or  High  Speed  Stop  (helicopter) 

Test  Stand  Runup  (w/mech.  or  man.  gov.  ) 

Test  Stand  Runup  (w/elec.  (auto)  gov.) 

Windup  During  Start  (before  eng.  fires) 

Windup  During  Start  (after  eng.  fires  and  before 

idle  speed  is  reached) 

Cruise  (unpressurized  —  abnormal  cond.  ) 
Ground  Runup  (w/one  or  more  eng.  out) 

Ground  Runup  (w/flaps  extended) 

Stall  (high  or  low  speed) 

Climb  (w/after  burner) 


Table  9  (continued) 


Description  of  Flight  Conditions 


Code 

Number 


Cruise  (A/B,  clean) 

Cruise  (A/B,  speed  brakes  extended) 

Cruise  (A/B,  flaps  extended) 

Ground  Runup  (A/B  or  water  injections) 

Ground  Runup  (A/B,  water  inject. ,  with  flaps  extended) 
Cruise  (flaps  and  gear  extended) 

Cruise  (flaps  and  speed  brakes) 

Cruise  (speed  brakes  and  gear  extended) 

Cruise  (speed  brakes,  flaps  and  gear  extended) 

Cruise  (A/B,  with  gear  extended) 

Cruise  (A/B,  flaps,  gear  extended) 

Cruise  (A/B,  flaps,  and  speed  brakes) 

Cruise  (A/B,  speed  brakes  and  gear  extended) 

Cruise  (A/B,  speed  brakes  and  flaps  extended) 

Descent  (with  gear  extended) 

Descent  (with  flaps  extended) 

Descent  (with  gear  and  flaps  extended) 

Descent  (with  gear  and  speed  brakes  extended) 

Descent  (with  flaps  and  speed  brakes  extended) 


56 


t 

I 

4.  2  DATA  ACQUISITION  EQUIPMENT 

AND  PROCEDURES 

Data  acquisition  procedures  and  techniques  often  have  a  profound 
impact  on  the  quality  of  the  resulting  data.  Hence,  it  is  important  to 
review  and  clarify  the  various  data  acquisition  instruments  and  pro- 
cedi  js  used  by  AFFDL  to  acquire  the  data  summarized  in  Table  7. 
Fortunately,  AFFDL  has  been  quite  consistent  in  the  instrumentation 
and  procedures  used  to  acquire  flight  acoustic  and  vibration  data.  The 
same  basic  procedures  were  used  to  acquire  data  in  all  the  aircraft 
listed  in  Table  7,  with  the  exception  of  aircraft  Code  Nos.  62  and  63 
(UH-IF  and  RF-4C).  At  a  point  in  history  just  prior  to  the  flight 
surveys  of  aircraft  Code  Nos.  62  and  63,  a  basic  change  was  made  in 
the  vibration  data  acquisition  system,  and  acoustic  data  acquisition 
was  introduced.  Further  changes  in  the  data  acquisition  system  are 
currently  being  contemplated.  For  the  case  of  currently  available  data, 
however,  the  vibration  data  acquisition  systems  used  for  all  flight 
vibration  surveys  may  be  divided  into  two  categories.  The  principal 
change  in  the  vibration  data  acquisition  system  was  from  velocity 
transducers  to  acceleration  transducers.  Hence,  the  two  different 
vibration  data  acquisition  systems  will  be  referred  to  as  the  velocity 
transducer  system  and  the  acceleration  transducer  system. 

4.  2.  1  Velocity  Transducer  System 

The  flight  vibration  surveys  of  all  aircraft  in  Table  7,  excluding 
aircraft  Code  Nos.  62  and  63,  were  performed  using  velocity  trans¬ 
ducers.  CEC  Model  4-102A-MD  velocity  pickups  were  used  in  aircraft 
Code  No.  49,  and  MB  type  124  velocity  pickups  were  used  for  the  other 
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aircraft.  The  pickups  were  mounted  at  various  locations  on  the  aircraft 
structure,  generally  in  clusters  of  three  so  that  vibration  along  the  three 
orthogonal  axes  would  be  sensed.  The  output  signals  from  the  pickups 
were  recorded  12  at  a  time  on  a  Davies  Model  501  14-channel  magnetic 
tape  recorder.  Since  there  were  usually  more  than  12  pickups  installed 
for  any  given  flight,  it  was  necessary  to  restrict  the  vibration  measure¬ 
ments  at  any  given  time  to  12  or  less  pickups.  By  using  a  selector 
switch,  data  from  all  pickups  were  recorded  in  a  grouped  sequential 
manner.  The  length  for  each  sample  record  was  fixed  at  5  seconds. 

The  exact  time  of  each  vibration  measurement  was  recorded  along  with 
pertinent  flight  parameters  including  engine  speed,  indicated  airspeed, 
altitude,  manifold  pressure  or  power  lever  angle,  percent  of  rated 
horsepower  or  percent  of  rated  thrust,  and  flight  condition  (takeoff, 
cruise,  etc. ). 

The  MB  Type  124  velocity  pickup  is  a  velocity  signal  generating 
device  with  a  nominal  sensitivity  of  96.  4  millivolts  per  inch  per  second, 
a  usable  frequency  range  of  5  to  2000  cps,  and  a  usable  temperature 
range  of  -65°  to  250° F.  A  typical  frequency  response  curve  is  shown 
in  Figure  5.  The  CEC  Model  4-102A-MP  has  somewhat  similar  char¬ 
acteristics. 

The  Davies  Model  501  recorder  is  an  FM  type  magnetic  type  tape 
recorder  with  a  carrier  frequency  of  10  kc,  a  usable  frequency  range  of 
3  to  2000  cps,  a  nominal  dynamic  range  of  45  dB,  and  a  tape  speed  of 
15  inches  per  second.  The  frequency  response  is  relatively  flat  over 
the  usable  frequency  range. 

Further  information  on  both  the  vibration  transducers  and  the  mag¬ 
netic  tape  recorder  are  available  in  any  of  the  documents  referenced  in 
Table  7,  excluding  the  reports  for  aircraft  Code  Nos.  62  and  63. 


58 


SENSITIVITY,  mv /inch/second 


T 


, 


1000 


FREQUENCY,  CPS 


Figure  5.  Frequency  Response  of  MB  Type  1Z4  Velocity  Pickup 


4.  2.  2  Acceleration  Transducer  System 

The  flight  vibration  surveys  on  aircraft  Code  Nos.  62  and  63 
were  performed  using  Columbia  Research  Laboratories  Model  902-H 
accelerometers  and/or  Endevco  Model  2215,  221  5C  and  2245B  accel¬ 
erometers.  The  accelerometers  were  used  in  conjunction  with  Endevco 
Model  2617  amplifiers.  As  for  the  velocity  pickups,  the  accelerometers 
were  mounted  at  various  locations  in  the  aircraft,  usually  in  clusters  of 
three,  to  sense  vibration  along  the  three  orthogonal  axes.  The  output 
signals  from  the  accelerometer  amplifiers  were  sequentially  recorded 
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The  Honeywell  Model  MTR-7200  recorder  is  an  FM  type  magnetic 
tape  recorder  with  a  carrier  frequency  of  54  kc,  a  usable  frequency 
range  of  dc  to  10,  000  cps,  a  nominal  dynamic  range  of  45  uB,  and  a 
tape  speed  of  60  inches  per  second.  The  frequency  response  is  relatively 
flat  over  the  usable  frequency  range. 

Further  information  on  both  the  acceleration  transducers  and  the 
magnetic  tape  recorder  are  presented  in  [26],  [28]. 

4.  2.  3  Microphone  Transducer  System 

The  flight  acoustic  surveys  on  aircraft  Code  Nos.  62  and  63  were 
performed  using  Gulton  Model  MA  299501  microphones,  in  conjunction 
with  Endevco  Model  2617  amplifiers.  The  microphones  were  mounted 
in  various  compartments  of  interest  to  sense  the  acoustic  noise  environ¬ 
ment  in  those  compartments.  The  output  signals  from  the  microphone 
amplifiers  were  recorded  on  a  Honeywell  .Model  MTR-7200  14-channel 
magnetic  tape  recorder  (the  same  instrument  used  to  record  the  output 
of  the  acceleration  transducers).  Record  lengths  were  5  seconds  or  longer 
and  the  time  of  each  measurement  along  with  pertinent  flight  conditions 
were  recorded. 

The  Culton  microphones  are  piezoelectric  crystal  devices  with  an 
upper  sound  pressure  level  limit  of  190  dB.  They  have  a  minimum  open 
circuit  sensitivity  of  1. 6  millivolts  for  120  dB  sound  pressure  level,  a 
usable  frequency  range  of  2  to  6000  cps,  and  a  vibration  sensitivity  of 
90  dB  sound  pressure  level  per  1  g  (rms)  vibration.  The  Model  2617 
amplifiers  are  ac  voltage  type  amplifiers.  The  frequency  response  for 
the  microphones  is  flat  to  'vithin  2  dB  over  the  usable  frequency  range. 

The  Honeywell  Model  MRT-7200  recorder  is  as  described  in 
Section  3.  2.  2.  Further  information  on  the  microphones  is  presented 
in  [26] ,  [28] . 
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4.  3  DATA  PROCESSING  PROCEDURES 

As  for  the  case  of  the  data  acquisition  equipment  and  procedures, 
the  data  processing  instruments  and  procedures  have  a  considerable 
influence  on  the  final  quality  of  the  data.  The  procedures  used  by  AFFDL 
to  reduce  vibration  data  have  been  consistent  for  that  data  collected 
using  velocity  transducers.  Specifically,  that  vibration  data  collected 
for  all  aircraft  flight  vibration  surveys,  excluding  the  surveys  of  air¬ 
craft  Code  Nos.  62  and  63,  were  performed  using  procedures  detailed 
in  WADC  Technical  Note  59-44  [29].  Different  instruments  and  some¬ 
what  different  data  reduction  procedures  were  introduced  at  a  time  prior 
to  the  flight  vibration  surveys  on  aircraft  Code  Nos.  62  and  63.  Further¬ 
more,  data  reduction  for  acoustic  data  were  introduced  at  that  time. 

4.  3.  1  Analysis  Procedures  for  Velocity  Data 


For  the  case  of  data  collected  using  velocity  transducers,  the 
reels  of  recorded  data  were  first  edited  in  the  laboratory  and  then  spliced 
into  endless  loops.  Each  loop  was  placed  on  a  Davies  Model  502  Tape 
Playback  System  and  continuously  recirculated  for  analysis.  A  frequency 
analysis  was  performed  on  the  data  using  a  6- channel  Davies  Model  510 
heterodyne  type  wave  analyzer.  The  bandwidths  for  the  analysis  were 
approximately  10  cps  in  the  frequency  range  below  500  cps,  and  30  cps 
in  the  frequency  range  above  500  cps.  The  averaging  time  constant  for 
the  analyzer  detector  circuit  was  approximately  1  second.  A  plot  of 
vibration  amplitude  versus  frequency  was  recorded  on  a  Brown  strip  chart 
recorder  for  each  channel  of  vibration  data.  The  resulting  frequency 
spectra  were  converted  into  discrete  data  points  for  handling  by  digital 
processing.  To  accomplish  this,  each  significant  peak  observed  in  each 


measured  frequency  spectrum  was  considered  to  be  a  data  point. 

Personal  judgment  was  used  to  define  a  significant  peak.  The  result¬ 
ing  data  points  were  punched  onto  IBM  cards.  Additional  cards  were 
punched  to  furnish  tie-in  data  which  described  the  pickup  locations,  the 
aircraft,  the  flight  conditions,  and  other  pertinent  information.  All 
information  was  then  handled  automatically  on  a  computer.  The  value 
of  each  spectral  peak  was  converted  from  a  velocity  to  a  double  ampli¬ 
tude  (peak-to-peak  displacement)  in  inches,  using  calibrations  which 
would  be  appropriate  if  the  spectral  peak  represented  a  sine  wave. 

The  resulting  data  were  then  presented  as  scatter  plots  of  double  ampli¬ 
tude  in  inches  versus  frequency  in  cps.  These  data  are  presented  in 
the  applicable  references  listed  in  Table  7,  [16]  -  [25],  [27].  Further 
descriptions  of  the  procedures  are  also  given  in  these  references,  as 
well  as  [29] . 

4.  3.  2  Analysis  Procedures  for  Acceleration  Data 

For  the  case  of  aircraft  Code  No.  62,  the  acceleration  data  re¬ 
cordings  were  edited  in  the  laboratory  and  then  spliced  into  endless  loops. 
Each  loop  was  placed  on  a  Honeywell  Model  3170  Tape  Reproducer 
System  and  continuously  recirculated  for  analysis.  The  data  were  ana¬ 
lyzed  using  a  6-channel  Honeywell  Model  9050  Automatic  Wave  Analyzer. 
The  analyses  were  performed  using  bandwidths  of  10  cps  in  the  frequency 
range  from  3  to  300  cps,  30  cps  in  the  frequency  range  from  300  to 
1000  cps,  and  100  cps  in  the  frequency  range  from  1000  to  10,000  cps. 
The  averaging  time  constant  for  the  analyzer  detector  circuit  was  1  sec¬ 
ond  in  all  cases.  The  resulting  data  were  reduced,  processed,  and 
plotted  exactly  as  was  done  for  the  velocity  data  discussed  in  Section 
4.  3.  1  (see  [28]  for  further  details). 
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For  the  case  of  the  acceleration  data  acquired  for  aircraft  Code 
No.  63,  the  data  were  edited,  spliced  into  endless  loops,  and  recirculated 
for  analysis.  The  data  were  analyzed  using  the  Honeywell  Model  9050 
Wave  Analyzer  and  a  Bruel  fcKjaer  Type  2111  Audio  Frequency  1/ 3-octave 
spectrometer. 

For  the  analysis  performed  using  the  wave  analyzer,  the  analysis 
bandwidths  were  10  cps  in  the  frequency  range  from  3  to  500  cps,  and 
40  cps  in  the  frequency  range  from  500  to  2400  cps.  The  averaging  time 
constant  for  the  wave  analyzer  detector  circuit  was  approximately  3  sec¬ 
onds.  The  resulting  frequency  spectra  were  converted  to  discrete  data 
points  as  before,  where  each  observed  peak  in  the  spectrum  was  con¬ 
sidered  a  data  point.  However,  in  this  case,  the  resulting  data  points 
were  converted  to  acceleration  power  spectral  density  levels  in  g^/cps, 
using  calibrations  which  would  be  appropriate  if  the  spectral  peaks 
represented  narrowband  Gaussian  random  noise.  No  effort  was  made 
to  distinguish  between  those  spectral  peaks  due  to  sine  waves  and  those 
due  to  narrowband  random  vibrations. 

For  the  analysis  performed  using  the  1 /3-octave  band  analyzer, 
the  data  were  analyzed  in  the  21  1/3  octaves  between  25  and  2500  cps. 

A  true  rms  detector  circuit  was  used.  The  resulting  frequency  spectra 
were  presented  as  rms  g's  versus  1 /3-octave  center  frequencies , (see 
[28]  for  further  details). 

4.  3.  3  Analysis  Procedures  for  Microphone  Data 

For  the  case  of  aircraft  Code  No.  62,  the  acoustic  data  recordings 
were  processed  and  presented  exactly  like  the  vibration  data,  as  dis¬ 
cussed  in  Section  4,  3.  2.  For  the  case  of  aircraft  Co.  No.  63,  the  acoustic 
data  were  processed  and  presented  exactly  like  that  vibration  data  analyzed 
using  the  1 /3-octave  analyzer,  as  discussed  in  Section  4.  3.  2  (see  [26] , 

[28]  for  further  details). 
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5.  PRELIMINARY  DATA  EVALUATION 


The  preliminary  evaluation  of  the  available  AFFDL  flight  vehicle 
vibration  data  immediately  led  to  an  important  conclusion.  Specifically, 
with  the  exception  of  the  RF-4C  survey  [26],  the  data  as  acquired, 
analyzed  and  presented  by  AFFDL  tends  to  accentuate  the  contribution 
of  periodic  vibration  sources  and  obscure  the  contributions  of  random 
vibration  sources.  In  other  words,  the  dc  ta  may  present  a  reasonable 
measure  of  periodic  vibrations,  but  they  definitely  provide  a  poor  mea¬ 
sure  of  random  vibrations.  The  reasons  for  this  evolve  principally 
from  the  data  analysis  procedures,  although  the  data  acquisition  pro¬ 
cedures  also  contribute. 

The  facts  which  led  to  the  above  conclusion  were  first  discovered 
during  initial  efforts  to  perform  regression  analysis  on  data  for  jet 
bombers  and  fighters  [19]  to  [25].  These  initial  regression  studies 
failed  to  indicate  a  correlation  between  structural  vibration  and  dynamic 
pressure  (q) .  It  is  interesting  to  note  that  this  result  was  observed 
during  a  previous  study  of  the  same  data  [30] ,  but  the  investigator  in 
that  case  chose  to  accept  the  result.  The  result  was  not  accepted  here 
since  past  studies  of  data  for  other  aircraft  [8]  ,  as  well  as  studies  of 
the  RF-4C  data  [31],  clearly  revealed  a  correlation  between  structural 
vibration  and  dynamic  pressure.  Furthermore,  the  theoretical  argu¬ 
ments  supporting  the  existence  of  this  correlation,  as  developed  in  Sec¬ 
tion  3,  are  quite  compelling.  Hence,  a  detailed  study  of  the  basic  quality 
of  the  AFFDL  vibration  data  was  pursued.  This  study,  in  turn,  led  to 
the  noted  conclusion.  Since  this  conclusion  is  profoundly  important  to 
the  determination  of  appropriate  data  evaluation  and  statistical  analysis 
prf>cedures  discussed  later,  it  is  desirable  to  clearly  present  all  of  the 
related  facts. 
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5.  1  SUMMARY  OF  DATA  LIMITATIONS 

There  are  a  number  of  factors  in  the  data  acquisition  and  analysis 
equipment  and  procedures  which  individually  tend  to  exaggerate  the  im¬ 
portance  of  periodic  contributions  relative  to  random  contributions  in 
the  vibration  environment  for  all  aircraft  surveyed  by  AFFDL,  except 
for  the  RF-4C  [26] .  They  may  be  summarized  as  follows. 

1.  Use  of  velocity  signal  generators 

2.  Frequency  response  function  for  velocity  signal 
generators 

3.  Use  of  common  sensitivity  for  all  recorder  inputs 

4.  Reduction  of  data  in  terms  of  the  amplitude  and 
frequency  for  spectral  peaks 

Each  of  these  contributing  factors  is  now  discussed. 

5.  1.  1  Use  of  Velocity  Transducers 

As  noted  in  Section  4.  2,  most  of  the  AFFDL  flight  vehicle  vibra¬ 
tion  data  have  been  acquired  using  velocity  type  transducers.  However, 
acceleration  type  transducers  would  have  been  more  appropriate  for  the 
following  reason.  There  is  a  general  tendency  for  aircraft  vibration  data 
to  have  velocity  spectra  which  decrease  with  increasing  frequency.  On 
the  other  hand,  the  acceleration  spectra  often  increase  with  increasing 
frequency  over  much  of  the  frequency  range  of  the  data.  This  fact  is 
clearly  illustrated  in  Figure  7  which  summarizes  data  for  Zone  02  of 
aircraft  Group  5,  Figure  7  covers  data  only  in  the  frequency  range  be¬ 
low  500  cps,  but  data  from  [21]  indicate  the  trend  continues  up  to  at  least 
2000  cps.  These  data  support  the  conclusion  that  velocity  transducers 


66 


AVERAGE  VELOCITY,  INCHES/SEC 


Figure  7.  Velocity  and  Acceleration  Spectra  of  Vibration  Data  for  Zone  2  of  Aircraft  Group  5 


produce  data  with  a  better  signal-to-noise  ratio  at  the  lower  frequencies 
and  a  poorer  signal-to-noise  ratio  at  the  higher  frequencies  than  would 
be  obtained  using  acceleration  transducers.  From  the  developments 
in  Section  3,  the  periodic  contributions  in  aircraft  vibration  data  occur 
principally  in  the  frequency  range  below  300  cps,  while  the  random 
contributions  are  prevalent  in  the  frequencies  above  this  range.  Hence, 
the  use  of  velocity  transducers  favors  the  periodic  contributions  over 
the  random  contributions.  The  use  of  accelerometers  will  generally 
produce  the  reverse  situation.  This  is  more  desirable  since  periodic 
components  can  be  retrieved  by  narrow  bandwidth  frequency  analysis, 
even  when  they  are  well  down  into  the  data  acquisition  system  noise. 

The  random  portions  of  the  data,  however,  cannot  be  retrieved  from 
the  noise.  Random  signals  can  be  extracted  from  random  noise  only  by 
cross -correlation  analysis  with  a  noise-free  replica  of  the  random  data. 


5.  1.  2  Frequency  Response  for 
Velocity  Transducers 

The  problem  discussed  in  Section  5. 1  is  compounded  by  the  fre¬ 
quency  response  characteristics  of  the  velocity  transducers.  Referring 
to  Figure  5,  it  is  clear  that  the  velocity  transducer  attenuates  the  signal 
being  measured  with  increasing  severity  as  the  frequency  increases 
above  100  cps.  For  example,  at  2000  cps,  the  signal  has  been  attenu¬ 
ated  by  more  than  6  dB.  However,  the  data  acquisition  system  noise  is 
not  attenuated.  Hence,  the  signal-to-noise  ratio  is  further  reduced  at 
the  higher  frequencies  where  the  random  contributions  in  the  data  are 
prevalent. 
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5.  1.  3  Use  of  Fixed  Sensitivity  into  all 
Channels  of  the  Recorder 


For  all  flight  vehicle  vibration  surveys  with  velocity  transducers, 
the  sensitivities  for  the  inputs  to  the  14- channel  magnetic  tape  recorder 
were  the  same  for  all  channels.  Hence,  the  recorder  sensitivities  were 
established  to  measure  the  maximum  signal  level  to  be  expected  for  any 
flight  condition  from  any  transducer  mounted  on  the  aircraft.  This 
means  that  flight  conditions  and  transducer  locations  of  relatively  low 
vibration  severity  resulted  in  transducer  signals  which  were  recorded 
well  down  from  the  maximum  input  level  to  the  recorder.  Since  the  re¬ 
corder  has  a  usable  dynamic  range  of  only  45  dB,  these  lower  level 
signals  were  often  recorded  near  or  even  below  the  noise  level  of  the 
recorder.  In  such  cases,  the  random  portions  of  the  signal  were  lost 
in  the  recorder  noise.  On  the  other  hand,  the  periodic  components  wiere 
probably  retrieved  by  the  narrow  bandwidth  frequency  analysis,  as  dis¬ 
cussed  in  Section  5.  1.  1. 

5.1.4  Data  Analysis  and  Presentation  Procedures 

The  principal  problem  is  introduced  by  the  procedure  of  reducing  the 
data  into  discrete  values  for  the  significant  observed  spectral  peaks.  When  a 
spectrum  analysis  is  performed  on  any  mixed  signal  consisting  of  random 
and  periodic  components,  the  periodic  components  tend  to  be  more  pro¬ 
nounced  in  the  spectrum.  This  is  due  to  the  fact  that  the  periodic  com¬ 
ponents  constitute  finite  amounts  of  energy  in  infinitesimally  narrow 
bandwidths.  Hence,  when  an  analysis  is  performed  using  a  narrow  band¬ 
pass  filter,  the  full  magnitude  of  a  periodic  component  will  be  passed 
by  the  filter,  no  matter  how  narrow  the  bandwidth  may  be.  For  random 
components,  the  mean  square  output  of  the  narrow  bandpass  filter  will 
be  proportional  to  the  bandwidth  of  the  filter.  Theoretically,  the 
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periodic  components  in  a  mixed  random  and  periodic  signal  will  have 
infinite  power  spectral  density  at  the  frequency  of  the  periodic  com¬ 
ponents.  In  practice,  the  periodic  components  will  not  appear  with 
infinite  density,  but  they  will  appear  with  a  very  high  density  relative 
to  the  random  components,  assuming  the  analysis  bandwidth  is  relatively 
narrow.  This  point  is  illustrated  in  Figure  8,  which  compares  the 
analysis  of  a  sine  wave  signal  to  that  of  a  narrowband  random  signal 
of  equal  rms  value  and  center  frequency.  The  narrowband  random  sig¬ 
nal  represents  the  response  of  a  5%  damped  structural  resonance  to 
random  excitation.  In  both  cases,  the  signals  are  superimposed  on  band - 
limited  white  noise  representing  independent  instrument  noise.  When 
these  data  are  analyzed  with  a  10  cps  wide  filter  (as  was  used  for  data 
analysis  by  AFFDL),  the  indicated  spectral  peak  for  the  sine  wave  signal 
is  2.  7  times  (8.  6  dB)  higher  than  the  background  noise,  while  the  spec¬ 
tral  peak  for  the  random  signal  is  only  1.  35  times  (2.  6  dB)  higher  than 
the  background  noise.  Since  the  data  reduction  procedure  requires  a 
subjective  opinion  as  to  what  constitutes  a  significant  spectral  peak,  it 
is  clear  that  peaks  due  to  sine  wave  signals  are  more  likely  to  be  noted 
than  peaks  due  to  random  signals. 

5.  2  LIMITATIONS  IMPOSED  ON 

DATA  EVALUATION 

At  first  glance,  the  discussions  in  Section  5.  1  might  imply  that 
vibration  data  collected  by  AFFDL  prior  to  the  flight  vibration  survey  of 
aircraft  Code  No.  63  (to  be  called  the  "previous  data")  are  not  usable 
for  the  statistical  evaluations  proposed  herein.  This  conclusion,  however, 
is  not  completely  warranted.  It  certainly  is  true  that  the  previous  data 
are  not  suitable  for  the  ultimate  goal  of  developing  a  general  prediction 
model  for  aircraft  vibration,  as  described  in  Section  2.  Nevertheless,  the 


previous  data  are  usable  for  estimating  the  periodic  portions  cf  the 
environment.  Furthermore,  they  are  useful  for  illustrating  the  various 
statistic  evaluation  techniques  to  be  developed. 

In  light  of  the  above  discussions,  the  previously  reduced  data  will 
be  used  to  perform  the  homogeneity  tests  and  spatial  variation  studies 
in  Sections  6  and  7,  where  it  is  understood  that  the  results  are,  at  best, 
applicable  only  to  the  periodic  portions  of  the  aircraft  vibration  environ¬ 
ment.  Results  for  the  random  portion  of  the  vibration  environment  must 
be  determined  from  future  data  collected  and  analyzed  using  current 
equipment  and  procedures  (similar  to  those  employed  for  aircraft  Code 
No.  63).  The  statistical  procedures  needed  for  the  development  of  a 
general  prediction  model,  covered  in  Section  8,  will  be  illustrated  using 
data  from  aircraft  Code  No.  63,  only.  The  final  development  of  a  I 
useful  vibration  prediction  model  will  require  considerable  data  from 
future  aircraft  vibration  surveys. 


6.  HOMOGENEITY  STUDIES 


For  simplicity,  it  would  be  desirable  to  have  a  single  vibration 
prediction  model  apply  to  all  structural  locations  of  all  aircraft.  On  the 
other  hand,  greater  accuracy  can  usually  be  achieved  by  developing  a 
separate  model  for  various  different  structural  locations  and  different 
aircraft.  Clearly,  a  compromise  is  required.  One  method  of  approach¬ 
ing  this  compromise  is  to  start  with  vibration  data  divided  among  the 
various  structural  locations  and  aircraft  for  which  data  have  been  ob¬ 
tained.  These  data  can  then  be  investigated  for  homogeneity  (equivalence) 
to  determine  which  structural  locations  and/or  aircraft  produce  similar 
vibration.  In  this  manner,  the  total  number  of  classifications  of  data 
hopefully  can  be  reduced  to  the  minimum  number  needed  to  describe 
significant  variations  in  the  vibration  environment.  For  the  case  of 
AFFDL,  flight  vibration  data  are  initially  classified  by  direction 
(vertical,  lateral,  and  longitudinal),  structural  zone  (see  Table  8),  and 
aircraft  model  and  group  (see  Table  7).  Hence,  the  homogeneity  studies 
will  start  with  the  data  divided  among  these  various  classifications. 

Determining  whether  or  not  two  sets  of  data  are  homogeneous  is  not 
as  straightforward  as  it  might  appear  at  first  glance.  Remembering  that 
the  available  data  represent  only  finite  samples  of  the  vibration  environ¬ 
ment,  it  is  not  likely  that  two  sets  of  data  would  be  identical  even  if  the 
data  represented  the  same  vibration  environment.  Some  type  of  statistical 
hypothesis  test  is  needed  to  determine  if  observed  differences  between 
two  sets  of  data  are  actually  significant,  or  simply  the  result  of  sampling 
errors.  Such  statistical  tests  are  described  and  illustrated  in  Section  6.  1. 
These  tests  are  applied  to  all  previous  AFFDL  vibration  data  in  Section  6.  2 
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Note  that  the  data  used  for  these  homogeneity  studies  do  not  properly 
describe  the  random  portion  of  the  aircraft  vibration  environment, 
as  discussed  in  Section  5.  Hence,  at  best,  the  conclusions  will  apply 
only  to  the  periodic  portion  of  aircraft  vibration. 


6.  1  DEVELOPMENT  OF  TESTING  PROCEDURES 

There  are  numerous  possible  ways  to  test  the  homogeneity  of 
sample  data,  but  two  specific  procedures  appear  appropriate  for  the 
problem  at  hand.  These  procedures  are  (1)  the  chi-square  test,  and 
(2)  the  multiple  rank  test.  Both  of  these  testing  procedures  are  now 
described  and  illustrated  using  AFFDL  vibration  data. 


6.  1.  1  Chi-Square  Test 


Assume  two  or  more  sets  of  sample  data  are  available  which 
represent  the  structural  vibration  for  two  or  more  different  situations 
of  interest  (different  orthogonal  directions,  different  zones,  different 
aircraft  within  a  group,  or  different  aircraft  groups).  It  is  now  desired 
to  test  the  hypothesis  that  the  two  or  more  sets  of  sample  data  are  from 
the  same  population.  Let  the  data  in  the  ith  frequency  interval  be 
divided  into  m.  amplitude  intervals,  as  illustrated  in  Table  10a.  Further¬ 
more,  for  any  given  frequency  interval,  let  n  be  the  number  of  ampli- 
tude  values  in  the  kth  amplitude  interval  for  the  jth  set  of  data,  as 
shown  in  Table  10b. 


Table  10.  Illustration  of  Data  Format  for  Chi-Square  Test 
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10b.  Observed  Number  of  Amplitude  Values  in  One  Frequency  Interval  for 
s  Sets  of  Data 


Now  let  be  the  true  probability  for  an  amplitude  in  the  kth  interval 

for  the  ith  set  of  data.  The  maximum  likelihood  estimator  for  P.,  is 

jk 

given  by 


P.,  =  n  /n. 
jk  jk  j 


The  s  sets  of  data  will  be  homogeneous  if  P, ,  =  P,,  =  P,,  =  .  .  .  =  P  ,  ; 

XiC  miC  .DK  sk 

k  =  1,  2,  3,  ....  m . 

For  simplicity,  consider  the  case  where  s  =  2  (two  sets  of  data 
are  available).  The  homogeneity  hypothesis  is  P,,  =  P,.  »  k  =  1,  2,  3, 

.  .  .  m.  Under  this  null  hypothesis,  the  estimate  of  a  common  value  for 
P^^  and  P2j^  •  denoted  by  P^,  is 


^k  "  ^k^*' 


where  the  notation  n  means  n.  summed  over  j,|  n  ,  = 

•  k  jk  I  *k  jk 

The  appropriate  test  statistic  is  then  '  ^ 


^2  _  ^  <°lk  -  "l  ^  <-2k  -  ”2 

'i^l  ''l".k'''  dl  "2".k'" 


P(1  -  P)  2  "ik^^.k  ■  "l*" 


where  P  =  n^/n 


The  X  statistic  defined  in  Eq.  (35)  has  a  x  distribution  with  (m  -  1) 
degrees-of-freedom  for  large  n^  and  . 

For  the  general  case  where  s  sets  of  sample  data  are  to  be  com¬ 
pared,  an  overall  test  statistic  is  given  by 

(36) 

The  X”  statistic  now  has  a  x”  distribution  with  s  (m  -  1)  degrees-of- 

freedorn  if  n  is  large  and  n.  n  ,  /n  >  2  for  all  i  and  k  . 

1  •k  —  •' 


Suppose  that  it  is  desired  to  test  whether  or  not  the  distribution 
of  vibration  levels  in  Zone  01  is  the  same  as  that  in  Zone  02  in  400-499 
cps  frequency  interval  for  Group  5  type  aircraft.  The  results  are 
presented  in  Table  11a.  In  order  to  satisfy  the  condition  that  n.  n  ,  /n  >  2, 

j  •  k  — 

it  is  necessary  to  pool  classes  10,  11,  and  12  together.  Then  the  results 
in  Table  11b  are  obtained. 

From  Table  11,  P  =  266/328  =  0.  81 1 .  Hence,  using  Eq.  (35), 


P(1  -  P) 


20 


266  P 


=  (6.  524)  (6.413)  =  41.88 
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Table  11.  Illustration  of  Chi-Square  Test  for  Homogeneity 


Table  11a.  Data  for  Group  5  Aircraft 
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Table  11b.  Pooled  Data  for  Group  5  Aircraft 


Amplitude  Intervals, 

10 

^  inches 

Sets 

6 

8 

10 

13 

20 

25 

32 

40 

(Zones) 

4 

5 

7 

9 

12 

19 

24 

31 

39 

50 

Total 

01 

9 

21 

53 

30 

41 

42 

17 

12 

21 

20 

266 

02 

8 

10 

25 

8 

8 

2 

1 

62 

tk 

17 

31 

78 

38 

49 

44 

18 

12 

21 

20 

328 

r; 

0 

II 

I 

I 

I 

I 

n 

0 

D 


78 


E 

E 

E 


The  5th  percentage  point  (95th  percentile)  of  the  distribution  with 
9  degrees  of  freedom  equals  16.  9.  Since  41.  88  >  16.  9  ,  the  hypothesis 
that  Zone  01  and  Zone  02  of  aircraft  Group  5  have  the  same  distribution 
in  the  400-499  cps  frequency  interval  would  be  rejected  at  the  5%  level 
of  significance. 

6.  1.  2  Multiple  Rank  Test 
2 

The  X  test  discussed  in  Section  6.1.  1  can  be  used  to  simultane¬ 
ously  test  more  than  two  sets  of  data  for  homogeneity.  The  result  of  the 
test,  however,  is  a  simple  "yes"  or  "no.  "  For  the  case  of  a  "no"  answer 
(rejection  of  the  homogeneity  hypothesis),  there  is  no  information  pro¬ 
vided  as  to  which  set  or  sets  of  data  are  nonhomogeneous.  This  problem 

2 

could  be  overcome  by  performing  individual  x  tests  on  combinations  of 

sets  taken  two  at  a  time.  Another  approach  which  is  considered  more 

suitable  is  to  use  a  nonparametric  rank  test.  A  rank  test  will  reveal  not 

only  an  overall  lack  of  homogeneity  among  s  sets  of  data,  but  also  which 

of  the  sets  are  nonhomogeneous  with  the  others.  On  the  other  hand,  it 

2 

should  be  noted  that  a  rank  test  is  less  "powerful"  than  a  x  test;  that 
is,  the  rank  test  involves  a  greater  risk  of  accepting  a  false  hypothesis. 

Consider  the  case  where  s  >  2  sets  of  sample  data  are  to  be 
compared.  Let  the  sample  size  for  the  jth  set  be  m  .  Now,  take  all 
of  the  sample  values  as  a  whole  and  arrange  them  in  order  of  increasing 
numerical  value  (rank),  as  illustrated  on  the  next  page.  When  two 
samples  have  the  same  numerical  value  (a  tie),  the  average  rank  is 
assigned  to  both  values. 
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2  3 


5  6 


8  9 


Value  Zg  x^  y, 


X  =  numerical  values  from  set  1 
y  =  numerical  values  from  set  2 
z  =  numerical  values  from  set  3 


Let  Tj  be  the  sum  of  the  ranks  for  the  sample  values  of  the  jth  set  of 
data.  Any  two  sets  of  data  can  be  compared  using  the  statistic 


d  =  T./n.  -  T  /n  ;  j  u 
j  j  u  u 


To  be  more  general,  assume  Sj^  sets  of  data  are  to  be  compared  to  s, 
different  sets  of  data.  The  desired  statistic  is  now 


From  [32]  ,  the  statistic  d  has  a  normal  distribution  with  a  mean 
value  of  zero  and  a  variance  of 


where  n  =^n.  +  ,  and  t  is  the  number  of  tied  scores  in  the  vth 

-Ui  j  ^  u  V 

rank.  For  example,  the  variance  of  d  given  by  Eq.  (37)  is 


(40) 


Hence,  a  null  hypothesis  of  homogeneity  can  be  tested  using  standard 
normal  distribution  tables  where  the  critical  region  is  given  by 


>  z 


a/2p 


(41) 


Here,  z  is 
o/2p 


o/Zp  percentage  point  of  the  standardized  normal  distri¬ 


bution  and  p  is  the  number  of  (multiple)  comparisons. 


To  illustrate  these  ideas,  assume  there  are  three  sets  of  data 
to  be  compared.  Further  assume  it  is  desired  to  compare  the  first 
set  of  data  with  the  second  and  third  combined,  the  first  and  the  third 
with  the  second,  and  the  first  with  the  third.  The  statistics  to  be 
calculated  are  then 


Example 


dj  =  Tj/n^  -  (T2+  T3)/(n2  +  n^) 


d^  =  (T^  +  T^)/(n^  +  n^)  -  T^/n^ 


^3  “  ■  ^3^”3 


The  Rank  Sum  Test  is  now  illustrated  using  the  same  data  used 
to  illustrate  the  test  in  Section  6.  1.  1.  Here,  s  =  2  and  p  =  1 .  Tl 


necessary  calculations  are  summarized  in  Table  12. 


Table  12.  Illustration  oC  Rank  Sum  Test 


From  Table  12,  it  follows  that 


d  =  47953/266  -  6003/62  =  180.  27  -  96.  82  =  83.  45 


2  r  (328)(329)  787,  192-328]  [ill 

'd  ^  L  12  '  (12)(327)  J  [  266  62  J 


=  174.  86 


and  hence 


d/o-j  =  83.45/13.  22  =  6.  31 
d 


Since  6.  31  >  1. 645 ,  the  hypothesis  that  the  vibration  in  Zones  01  and  02 
of  aircraft  Group  V  have  the  same  distribution  in  400  to  499  cps  frequency 
interval  is  rejected  at  the  5%  level  of  significance. 
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6.  2  APPLICATION  TO  AFFDL  VIBRATION  DATA 


As  previously  mentioned,  the  aircraft  vibration  data  collected 
by  AFFDL  is  initially  classified  in  terms  of  direction,  structural  zone, 
and  aircraft  model  and  group.  These  various  classifications  will  be 
used  as  the  basis  for  homogeneity  studies.  The  rank  sum  test  de¬ 
scribed  in  Section  6.  1.  2  will  be  used  for  the  desired  comparisons, 
as  required.  To  assure  consistency  in  the  data  used  for  the  com¬ 
parisons,  only  that  data  obtained  from  velocity  transducers  will  be 
employed.  This  eliminates  all  data  collected  on  aircraft  Code 
Nos.  62  and  63.  Furthermore,  only  data  in  the  first  six  structural 
zones  (Zone  Code  Nos.  01  through  06)  will  be  used.  These  six 
zones  cover  the  bssic  aircraft  structure,  which  is  of  principal  interest. 
Separate  tests  will  be  performed  on  the  data  in  each  of  20  contiguous 
frequency  increments  which  together  cover  the  frequency  range  from  5 
to  500  cps.  For  simplicity  in  later  presentations,  these  various  fre¬ 
quency  increments,  as  well  as  the  different  directions,  aircraft,  and 
structural  zones,  will  be  identified  by  a  code  .number.  These  code 
numbers  are  s\immarized  in  Table  13. 

6.  2.  1  Comparisons  of  the  Three  Orthogonal  Directions 

Consider  first  the  vibration  levels  measured  along  the  three 
orthogonal  axes.  The  question  of  interest  is  as  follows.  Are  there 
significant  differences  among  the  vibration  levels  measured  in  aircraft 
along  the  vertical,  lateral,  and  longitudinal  directions?  To  answer 
this  question,  rank  sum  tests  are  applied  to  the  data  segregated  among 
the  three  orthogonal  directions.  Separate  tests  are  performed  on  the 
data  for  each  frequency  increment  and  for  each  aircraft  group  .  The 
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data  for  the  different  structural  zones,  different  aircraft  within  a  group, 
and  different  flight  conditions  are  pooled  together  for  the  tests  to  maxi¬ 
mize  the  sample  size.  This  pooling  is  considered  appropriate  since 
vertical,  lateral,  and  longitudinal  measurements  were  made  in  all  zones 
of  all  aircraft  for  all  flight  conditions.  Note  that  only  data  from  tri- 
axial  measurement  locations  are  used  to  avoid  possible  bias  in  the  re¬ 
sults. 

The  detailed  results  of  the  homogeneity  tests  are  presented  in 
Section  1  of  Appendix  A,  and  summarized  in  Table  14.  The  results  are 
presented  in  terms  of  relative  severity  (A,  B,  C)  versus  frequency  for 
the  vibration  measured  in  the  three  directions  for  each  aircraft  group. 

A  common  level  of  severity  is  assigned  to  those  data  which  are  found 
to  be  homogeneous  at  the  1%  level  of  significance.  For  example,  if  the 
vibration  in  the  vertical  and  lateral  directions  are  found  to  be  equivalent 
at  the  1%  level  of  significance,  and  the  vibration  in  the  longitudinal  direc¬ 
tion  is  found  to  be  significantly  lower,  then  the  vertical  and  lateral  vibra¬ 
tion  are  identified  as  A  and  the  longitudinal  vibration  is  identified  as  B  . 

The  rank  sum  tests  establish  only  the  relative  severity  of  the  data 
and  not  their  absolute  value.  An  absolute  measure  of  the  data  severity 
is  given  by  the  mean  value  and  standard  deviation  versus  frequency  for 
the  vibration  measured  in  the  three  directions  for  each  aircraft  group. 
These  statistical  measures  are  calculated  in  Section  4  of  Appendix  A, 
and  summarized  in  Figure  9.  Note  that  the  results  in  Figure  9  do  not 
consider  possible  homogeneity.  The  sample  sizes  for  the  mean  value 
and  standard  deviation  calculations  vary  from  4  to  4800,  depending  upon 
the  direction,  frequency  increment,  and  aircraft  group.  Hence,  a  large 
indicated  difference  in  two  mean  values  may  not  be  significant  in  some 
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Table  14.  Rank  Sum  Test  Results  for  Direction  Comparisons 
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Figure  9a.  Mean  Value*  and  Standard  Deviations  for  the  Vibration 
in  the  Three  Orthogonal  Directions  for  Aircraft  Group  1 
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Figure  9b.  Mean  Values  and  Standard  Deviations  for  the  Vibration 
in  the  Three  Orthogonal  Directions  for  Aircraft  Group  2 
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Figure  9c.  Mean  Values  and  Standard  Deviations  for  the  Vibration 
in  the  Three  Orthogonal  Directions  for  Aircraft  Group  3 
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Figure  9d.  Mean  Values  and  Standard  Deviations  for  the  Vibration 
in  the  Three  Orthogonal  Directions  for  Aircraft  Group  5 
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cases,  while  a  small  indicated  difference  might  be  significant  in  other 
cases.  The  results  in  Figure  9  can  only  be  interpreted  based  upon  the 
conclusions  in  Table  14.  This  is  done  in  Figure  10. 

Figure  10  summarizes  the  mean  value  versus  frequency  for  the 
vibration  measured  in  the  three  directions  for  each  aircraft  group, 
where  homogeneous  data  are  pooled  and  identified  by  a  single  mean 
value.  Hence,  all  differences  indicated  in  Figure  10  are  significant 
differences  tin  the  statistical  sense. 

Inspection  of  the  results  in  Figure  10  immediately  reveals  an  im¬ 
portant  conclusion.  Specifically,  aircraft  vibration  levels  measured 
along  the  three  orthogonal  axes  appear  to  be  quit^  similar  on  the  average. 
From  Figure  10a,  for  aircraft  Group  1  (reciprocating  engine  transports), 
the  vibration  in  the  vertical  direction  is  somewhat  higher  than  in  the 
lateral  and  longitudinal  directions  for  frequencies  between  50  and  80  cps. 
This  is  the  frequency  range  for  blade  passage  excitation.  On  the  other 
hand,  the  longitudinal  vibration  appears  to  be  the  highest  between  100  and 
200  cps.  From  Figure  10b,  for  aircraft  Group  2  (turboprop  transports), 
the  vertical  vibration  is  consistently  higher  for  frequencies  below  20  cps, 
while  the  lateral  vibration  is  usually  the  most  severe  for  frequencies 
above  50  cps.  From  Figure  10c,  for  aircraft  Group  3  (jet  bombers), 
the  vibration  levels  are  surprisingly  similar  in  all  three  directions  over 
the  entire  frequency  range.  From  Figure  lOd,  for  aircraft  Group  5 
(century  series  fighters),  there  is  a  clear  tendency  for  the  most  severe 
vibration  to  be  in  the  vertical  direction  for  frequencies  below  200  cps, 
and  in  the  lateral  direction  for  frequencies  above  200  cps.  From  Fig¬ 
ure  lOe,  for  aircraft  Group  10  (helicopters),  the  vibration  is  most  intense 
in  the  lateral  direction  for  frequencies  from  14  to  40  cps,  and  in  the 
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igure  10c.  Summary  of  Average  Vibration  in  the  Three  Orthogonal  Directions  for  Aircraft  Group  3 


HELICOPTERS 


vertical  direction  for  most  frequencies  above  65  cps.  For  all  aircraft 
groups,  however,  the  differences  from  lowest  to  highest  average  levels 
at  any  frequency  are  usually  less  than  two  to  one  (6  dB). 

The  results  in  Figure  10  have  definite  implications  to  the  data 
acquisition  problem,  as  well  as  to  the  test  specification  design  problem. 
First,  AFFDL  customarily  collects  triaxial  measurements  at  most 
measurement  locations  during  a  flight  vibration  survey.  Since  the  vibra¬ 
tion  in  the  three  orthogonal  directions  is  so  similar  on  the  average, 
it  appears  that  a  more  appropriate  procedure  might  be  to  collect  single 
direction  measurements  at  three  times  as  many  locations.  The  measure¬ 
ment  locations  would  be  divided  equally  on  a  random  basis,  among  vertical, 
lateral,  and  longitudinal  measurements.  This  should  greatly  increase  the 
amount  of  significant  data  obtained  from  a  given  survey.  Second,  the 
vibration  test  specifications  in  MIL  STD  81  OB  do  not  differentiate 
among  the  three  orthogonal  directions  (the  test  levels  are  the  same  for 
vertical,  lateral,  and  longitudinal  vibration).  Since  the  average  differ¬ 
ences  in  the  vibration  along  the  three  directions  are  probably  small 
relative  to  other  uncertainties  in  the  test  specification,  it  appears  that 
this  is  an  acceptable  procedure  for  testing  aircraft  components. 

It  should  be  emphasized  that  the  above  conclusions  apply  to  the 
periodic  portion  of  the  vibration  environments  on  the  basic  aircraft  struc¬ 
ture  (Zones  01  through  06).  The  extension  of  these  conclusions  to  the 
random  portion  of  the  environment  will  require  a  similar  analysis  of 
future  data  which  properly  represent  the  random  contributions  in  the  en¬ 
vironment. 

One  final  point  should  be  noted.  For  the  analyses  in  this  section, 
the  vibration  data  were  divided  among  the  three  orthogonal  directions 
relative  to  the  aircraft  reference  frame.  This  is  considered  proper  since 


I 
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vibration  tests  are  performed  using  excitation  applied  sequentially  in 
these  three  directions.  Another  method  of  direction  classification, 
however,  would  be  in  terms  of  two  directions;  namely,  perpendicular 
to  the  structure  and  tangential  to  the  structure.  For  example,  a 
vertical  measurement  on  the  bottom  of  the  aircraft  fuselage  would  be 
perpendicular,  while  a  vertical  measurement  on  the  sidewall  of  the 
aircraft  fuselage  would  be  tangential.  If  the  AFFDL  data  were  classi¬ 
fied  in  terms  of  perpendicular  and  tangential  directions,  it  is  quite 
likely  that  different  results  would  be  obtained.  Specifically,  the  per¬ 
pendicular  measurements  would  probably  be  substantially  higher  than 
the  tangential  measurements. 


Comparisons  of  the  Various  Structural  Zones 


Now  consider  the  vibration  levels  measured  in  the  basic  structural 
zones  (Zones  01  through  06).  The  question  of  interest  is  as  follows.  Are 
there  significant  differences  among  the  vibration  levels  measured  in  the 
basic  structural  zones  of  aircraft  as  defined  by  AFFDL?  The  answer  to 
this  question  is  pursued  using  procedures  similar  to  those  employed  in 
Section  6.2.1.  Rank  sum  tests  are  applied  to  the  data  segregated  among 
the  various  structural  zones.  Separate  tests  are  performed  on  the  data 
for  each  frequency  increment  and  for  each  aircraft  group.  The  data  for 
different  orthogonal  directions,  different  aircraft  within  a  group,  and 
different  flight  conditions  are  pooled  together  for  the  tests  to  maximize  the 
sample  size. 

The  detailed  results  of  the  homogeneity  tests  are  presented  in 
Section  2  of  Appendix  A,  an:^  summarized  in  Table  15.  As  before,  the 
results  are  presented  in  terms  of  relative  severity  versus  frequency  for 
the  various  structural  zones  in  each  aircraft  group.  A  common  level  of 
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Table  15.  Rank  Sum  Test  Results  for  Zone  Comparisons 
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severity  is  assigned  to  those  data  which  are  found  to  be  homogeneous 
at  the  1%  level  of  significance.  The  mean  value  and  standard  deviation 
versus  frequency  for  the  vibration  measured  in  each  aircraft  group  are 
calculated  in  Section  4  of  Appendix  A,  and  summarized  in  Figure  11. 

The  sample  sizes  for  the  mean  value  and  standard  deviation  calculations 
vary  from  2  to  over  4000,  depending  upon  the  aircraft  and  frequency 
increment.  Combining  the  information  in  Table  15  and  Figure  11  leads 
to  the  final  presentation  in  Figure  12. 

Figure  12  summarizes  the  mean  value  versus  frequency  for  the 
vibration  measured  in  the  various  zones  of  each  aircraft  group,  where 
homogeneous  data  are  pooled  and  identified  by  a  single  mean  value. 

Hence,  all  differences  indicated  in  Figure  12  are  significant  differences 
in  the  statistical  sense. 

From  Figure  12a,  for  aircraft  Group  1  (reciprocating  engine 
transports),  the  vibration  levels  in  the  forward  quarter  of  the  fuselage 
tend  to  be  more  severe  than  in  other  zones  for  frequencies  above  50  cps. 
The  levels  are  also  high  in  the  center  half  of  the  fuselage  for  frequencies 
between  50  and  100  cps.  This  would  be  expected  since  the  propeller  blade 
passage  frequency,  which  is  a  principal  source  of  fuselage  vibration, 
falls  within  this  frequency  range.  From  Figure  12b,  for  aircraft  Group  2 
(turboprop  transports),  the  middle  half  of  the  fuselage  again  displays  high 
vibration  levels  at  the  frequencies  of  propeller  blade  passage  and  its 
harmonics.  The  outer  one-third  of  the  wing  also  displays  high  vibration 
levels  at  these  frequencies.  From  Figure  12c,  for  aircraft  Group  3  (jet 
bombers),  the  vibration  does  not  differ  widely  among  the  basic  structural 
zones,  except  for  the  inner  two-thirds  of  the  wing  where  the  levels  are 
substantially  lower  at  some  frequencies  above  50  cps.  The  results  for 
aircraft  Group  5  (Century  series  fighters)  in  Figure  12d  are  similar. 
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Figure  lia.  Mean  Values  and  Standard  Deviations  for  the  Vibration 
in  Various  Structural  Zones  of  Aircraft  Group  1 


104 


VIBRATION  DOUBLE  AMPLITUDE 


9  19  90 


FREQUENCY  ,  CPS 
Figure  11a  (Continued) 
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Figure  11c.  Mean  Values  and  Standard  Deviations  for  the  Vibration 
in  Various  Structural  Zo  ’rs  of  Aircraft  Group  3 
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Figure  11c  (Continued) 
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Figure  12a.  Summary  of  Average  Vibration  for  Various  Structural  Zones  of  Aircraft  Group  1 


Figure  12b.  Summary  of  Average  Vibration  for  Various  Structural  Zones  of  Aircraft  Group  2 
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Figure  12e.  Summary  of  Average  Vibration  for  Various  Structural  Zones  of  Aircraft  Group  10 


From  Figure  12e,  for  aircraft  Group  10  (helicopters),  the  differences 
among  the  average  vibration  in  the  three  zones  are  less  than  10  dB  at 
all  frequencies. 

The  principal  purpose  of  zoning  is  to  segregate  the  aircraft  struc¬ 
ture  into  categories  where  the  vibration  levels  in  each  category  are 
similar.  For  example,  assume  an  aircraft  structure  is  to  be  divided 
into  three  zones.  Further  assume  that  the  vibration  levels  over  the 
structure  vary  on  an  arbitrary  scale  between  1  and  10  units.  Then  an 
optimum  zoning  procedure  would  segregate  the  structure  such  that  the 
vibration  levels  vary  from  1  to  4  units  in  Zone  1,  4  to  7  units  in  Zone  2, 
and  7  to  10  units  in  Zone  3.  Conversely  ,  a  completely  inefficient  zoning 
procedure  would  segregate  the  structure  such  that  the  vibration  levels 
vary  from  1  to  10  in  all  three  zones. 

Inspection  of  the  results  in  Figure  12  suggests  that  the  zones  used 
by  AFFDL  to  segregate  the  basic  aircraft  structure  are  not  very  efficient, 
particularly  for  aircraft  Groups  3,  5,  and  10.  It  must  be  r<“*.nembered, 
however,  that  the  data  used  to  prepare  Figure  12  do  not  properly  represent 
the  random  portion  of  the  vibration  environment.  If  random  contributions 
were  accurately  reflected  in  the  data,  it  is  believed  that  more  striking 
differences  would  be  seen  among  the  various  zones,  at  least  for  aircraft 
Groups  3  and  5.  For  jet  bombers  and  fighters,  the  random  vibration  in¬ 
duced  by  jet  noise  tends  to  be  more  severe  towards  the  aft  end  of  the 
fuselage.  Furthermore,  the  boundary  layer  induced  vibration  varies  in 
spectral  composition  over  the  fuselage  length.  Hence,  significant  differ¬ 
ences  would  be  expected  between  the  vibration  in  the  forward  quarter  and 
the  aft  quarter  of  the  fuselage  for  these  aircraft.  Based  upon  these  con¬ 
siderations,  no  change  in  the  current  zoning  procedure  is  suggested  until 
further  studies  are  performed  on  future  data  which  include  a  proper 
representation  of  the  random  portion  of  the  environment. 
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6.  2.  3  Comparisons  of  the  Various 
Aircraft  Models  and  Groups 


Attention  is  now  directed  to  the  average  vibration  environment 
for  each  of  the  11  aircraft  models  divided  among  the  five  aircraft  groups, 
as  detailed  in  Table  13.  The  question  of  interest  is  as  follows.  Are 
there  significant  differences  among  the  vibration  environments  for  the 
aircraft  within  a  given  group,  and  for  the  five  different  groups? 

For  the  case  of  aircraft  comparisons  within  a  given  group,  rank 
sum  tests  are  applied  to  the  data  segregated  among  the  various  aircraft 
within  each  group  covering  two  or  more  aircraft  (Groups  2,  3,  and  5). 
Separate  tests  are  performed  on  the  data  for  each  frequency  increment. 
The  data  for  the  different  orthogonal  directions,  different  structural 
zones,  and  different  flight  conditions  are  pooled  together  for  the  tests 
to  maximize  the  sample  size. 

The  detailed  results  of  the  homogeneity  tests  are  presented  in 
Section  3  of  Appendix  A,  and  summarized  in  Table  16.  The  results  are 
presented  in  terms  of  relative  severity  versus  frequency  for  the  various 
aircraft  in  each  aircraft  group,  where  a  common  level  of  severity  is 
assigned  to  those  data  which  are  found  to  be  homogeneous  at  the  1% 
level  of  significance.  The  mean  value  and  standard  deviation  versus 
frequency  for  the  vibration  measured  in  each  aircraft  of  each  group  are 
calculated  in  Section  4  of  Appendix  A,  and  summarized  in  Figure  13. 

The  sample  sizes  for  the  mean  value  and  standard  deviation  calculations 
vary  from  3  to  12,  250,  depending  upon  the  aircraft  and  frequency  incre¬ 
ment.  Combining  the  information  in  Table  16  and  Figure  13  leads  to  the 
final  presentation  in  Figure  14. 

Figure  14  summarizes  the  mean  value  versus  frequency  for  the 
vibration  measured  in  each  aircraft  of  each  group,  where  homogeneous 
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VIBRATION  DOUBLE  AMPLITUDE  ,  INCHES 
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Figure  13b.  Mean  Values  and  Standard  Deviations  for  the  Vibration  J 

in  Various  Aircraft  of  Group  2 
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Figure  13c.  Mean  Values  and  Standard  Deviations  for  the 
Vibration  in  Various  Aircraft  of  Group  3 
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Figure  13d.  Mean  Values  and  Standard  Deviations  for  the 
Vibration  in  Various  Aircraft  of  Group  5 
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Figure  13e.  Mean  Values  and  Standard  Deviations  for  the 
Vibration  in  Various  Aircraft  of  Group  10 
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Figure  14a.  Summary  of  Averag 


AIRCRAFT  GROUP  3  '  -  VC  36  (JRB-52) 


data  are  pooled  and  identified  by  a  single  mean  value.  Hence,  all 
differences  indicated  in  Figure  14  are  significant  differences.  Note 
that  no  summaries  are  presented  for  aircraft  Groups  1  and  10  since 
data  from  only  one  aircraft  are  available  for  these  two  groups. 

Referring  to  Figure  14a,  which  summarizes  aircraft  Group  2 
(turboprop  transports),  it  is  seen  that  the  average  vibration  levels  in 
the  C-130  are  up  to  5  times  (14  dB)  higher  than  in  the  C-133  in  the  fre¬ 
quency  range  above  40  cps.  Below  40  cps,  the  vibration  environments 
are  similar  in  the  two  aircraft.  From  Figure  14b,  for  aircraft  Group  3 
(jet  bombers), the  vibration  environments  for  the  three  aircraft  are 
quite  similar  up  to  about  150  cps.  Above  this  frequency,  the  vibration 
levels  in  the  B-58  are  somewhat  higher  than  in  the  JRB-52  and  JRB-66. 
From  Figure  14c,  for  aircraft  Group  5  (Century  series  fighte  rs),  the 
vibration  environments  for  the  four  aircraft  differ  widely  in  some  fre¬ 
quency  increments,  with  the  F-106  levels  being  generally  the  most 
severe.  In  the  frequency  range  above  65  <:ps,  however,  the  average 
vibration  levels  in  the  four  aircraft  differ  by  less  than  three  to  one 
(10  dB). 

Now  consider  the  average  vibration  environments  for  each  of  the 
five  aircraft  groups.  Homogeneity  tests  are  not  required  here  since  the 
sample  size  of  the  data  within  each  group  is  quite  large  (over  500  in  most 
frequency  increments).  The  mean  value  and  standard  deviation  versus 
frequency  for  the  vibration  measured  in  each  aircraft  group  are  calcu¬ 
lated  in  Section  4  of  Appendix  A,  and  summarized  in  Figure  15.  A  com¬ 
parison  of  the  mean  value  versus  frequency  for  the  five  groups  is  pre¬ 
sented  in  Figure  16. 

The  results  in  Figure  16  must  be  interpreted  with  caution.  Noting 
that  the  maximum  difference  in  average  vibration  among  the  five  groups 
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131 


VIBRATION  DOUBLE  AMPLITUDE  ,  INCHES 


GROUP  5  -  X 

--  S  +  X 


Figure  15  (Continued) 

H 

0 

E 

E 

I 

E 


is  less  than  10  dB  at  all  frequencies,  the  first  inclination  might  be  to 
conclude  that  the  vibration  environments  for  propeller  transports,  jet 
bombers  and  fighters,  and  helicopters  are  not  strikingly  different. 

Figure  16  indeed  suggests  that  this  may  be  true  on  a  gross  average  basis, 
at  least  for  the  periodic  portions  of  the  environments  below  500  cps. 
However,  it  certainly  is  not  true  on  a  detailed  basis.  For  example, 
consider  the  vibration  levels  in  the  frequency  increment  from  50  to 
64  cps  for  turboprop  transports  and  jet  bombers.  From  Figure  15,  the 
average  vibration  in  turboprop  transports  is  less  than  2  times  (6  dB) 
higher  than  in  jet  bombers.  Now  restrict  attention  to  the  center  half  of 
the  fuselage.  From  Figure  12,  for  the  same  frequency  increment,  it  is 
seen  that  the  average  vibration  in  the  center  half  of  the  fuselage  for  turbo¬ 
prop  transports  is  30  times  (30  dB)  higher  than  in  the  same  zone  for  jet 
bombers.  If  the  comparison  were  further  reduced  to  an  individual  turbo¬ 
prop  transport  versus  an  individual  jet  bomber,  even  greater  differences 
might  be  observed.  In  purely  statistical  terms,  the  above  argument  is 
supported  by  the  relatively  large  standard  deviations  for  the  data  in  each 
group,  as  illustrated  in  Figure  15.  In  conclusion,  the  results  in  Figure  16 
serve  only  to' prove  that  nearly  all  pertinent  information  in  any  data  can 
be  obscured  if  sufficient  averaging  is  performed.  One  is  reminded  at 
this  point  of  the  story  about  the  man  who  drowned  in  a  lake  with  an 
average  depth  of  3  feet. 


7.  SPATIAL  DISTRIBUTION  STUDIES 


The  mean  values  and  standard  deviations  for  various  classifications 
of  past  AFFDL  vibration  data  have  been  calculated  and  evaluated  in  Sec¬ 
tion  6.  Such  information,  however,  is  not  sufficient  to  arrive  at  test 
specifications  and  design  criteria.  Specifications  and  criteria  must  be 
based  upon  an  extreme  vibration  level  as  opposed  to  an  average  level. 

To  establish  an  extreme  level,  it  is  necessary  to  determine  a  distribution 
function  for  the  vibration  levels  within  each  structural  zone  to  be  covered 
by  the  specification  or  criteria.  It  would  be  most  convenient  if  a  single 
typo  of  distribution  function  could  be  found  which  would  fit  the  spatial 
variation  of  vibration  levels  in  all  zones.  It  would  be  further  convenient 
if  the  distribution  function  in  question  were  of  a  classic  form  which  is  well 
tabulated. 

7.  1  CONSTRUCTION  OF  A  DISTRIBUTION 
FUNCTION  MODEL 

In  the  past,  various  models  have  been  used  to  describe  the  spatial 
distribution  of  flight  vehicle  vibration  levels.  In  some  cases,  a  classical 
function  is  assumed  and  fitted  to  the  data.  Examples  of  functions  which 
have  been  assumed  include  the  lognormal  distribution  function  [Z]  and 
the  Rayleigh  distribution  function  [8].  In  other  cases,  an  empirical  func¬ 
tion  is  developed  by  analysis  of  actual  data  [33] .  The  most  common  of 
these  techniques  is  to  assume  the  power  spectra  of  the  vibration  data  fit 
a  lognormal  distribution  function.  It  appears  logical  that  the  first 
approach  to  the  problem  should  be  to  establish  whether  or  not  the  log¬ 
normal  distribution  is  an  appropriate  fit  to  the  spatial  distribution  of  air¬ 
craft  vibration  data  collected  by  AFFDL. 
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7.  1.  1  Definition  of  the  Lognormal  Distribution 


The  random  variable  x  is  said  to  have  a  lognormal  distribution 
if  y  =  log  X  is  normally  distributed.  That  is,  if  y  is  normal,  then 

y 

X  =  e  has  a  lognormal  distribution. 

The  probability  density  function  of  x  is  given  by 


p(x)  =  ^  P(x)  =  £  P(y)  =  ^  p  (log  X) 


1  (log  x-ti)^ 

2  2 
(T 


=  0 


for  x  <  0 


(42) 


where  P(y)  denotes  the  normal  distribution  function. 

The  lognormal  density  function  has  positive  skewness,  and  is 
defined  by  two  parameters.  The  mean  of  x,  denoted  by  ,  is  de¬ 
fined  in  terms  of  the  mean  and  standard  deviation  of  y,  denoted  by 
p  and  (T  ,  respectively,  as  follows. 


j  -CD 


e  p(y)  dy 
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Correspondingly,  the  variance  of  x,  denoted  by  v  ,  is  defined  in  terms 

X 

of  the  mean  and  standard  deviation  for  y,  as  follows. 
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It  is  important  to  note  that  the  transformation  e^  transforms  the  mean 
as  well  as  the  variance. 

If  the  median  of  x  is  denoted  by  m  ,  then 

X 


V  1 

P(x  <  m  )  =  P(e'  <  m  )  =  P(y  <  lo(^  m  )  =  — 


Thus 


log  m^  =  |x  or 


(45) 


Using  the  invariance  property  of  the  maximum  likelihood  estimate, 

2 

the  maximum  likelihood  estimate  of  u  ,  <r  and  m  can  be  shown  to  be  | 

XX  X 


a2 

<r 
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m  = 
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(46) 


where  y  and  s  denote  the  sample  mean  and  variance  of  the  normal 
distribution  of  y  =  log  x.  The  above  estimates  are  not  unbiased. 

A  convenient  method  of  computing 

P(a  <  X  <  b) 


is  to  observe  the  following. 


1 


P(a  <  X  <  b)  =  P(log  a  <  log  x  <  log  b) 


loe  a  - 


cr  — 


where  and  tr  are  the  mean  and  standard  deviation  of  the  normally 
distributed  random  variable,  y  =  log  x,  and^(z)  denotes  the  standard¬ 
ized  normal  distribution.  For  example,  if  cr  =  1  and  p  =  0 ,  then 


^  ^  j  =  P  (log  X  <  P  +  (T^/Z) 


=  f’  7  f  0.6915 


7.  1.  2  Applications  to  Aircraft 
Vibration  Environments 


Consider  a  given  structural  zone  of  a  particular  aircraft  group. 
Assume  the  vibration  data  have  been  reduced  into  the  form  of  mean 
square  values  (or  average  power  spectra)  in  contiguous  narrow  frequency 
intervals.  Let  L  be  the  midpoint  of  the  ith  frequency  interval.  Let  the 
data  in  the  ith  frequency  interval  be  divided  into  m.  mean  square  value 
inter ^'^als,  where  is  the  midpoint  of  the  kth  mean  square  value  interval 
for  the  data  in  the  ith  frequency  interval.  Hence,  the  segregation  of  the 


( 


%  . 


data  is  as  illustrated  previously  in  Figure  10a,  excepthere  x  is  a  mean 
square  value  rather  than  a  peak  value. 

Now  let  p{x^)  denote  the  probability  density  for  the  number  of 
mean  square  values  in  the  ith  frequency  interval.  If  the  mean  square 
values  in  the  ith  frequency  interval  are  lognormally  distributed,  then 
p(xp  is  given  by  (42)  as 


1 

PiCx.)  - - =  e  for  X.  >  0 

x.<r.y2v  ^ 

1  1  * 


for  X,  <  0 
1 


where  p.  and  <r.  are  the  mean  value  and  standard  deviation,  respectively, 

for  the  random  variable  y^  =  log  x^ . 

Let  n.  be  the  sample  size  in  the  ith  frequency  interval,  and  n.,  be 
1  ik 

observed  number  of  values  in  the  kth  mean  square  value  interval  of  the 

2 

ith  frequency  interval.  The  maximum  likelihood  estimates  of  u.  and  <r, 

1 

are  given  by 


<r.  = 

i 


n.,  y.,  /n. 
ik  ik  1 

k=l 

in,  \  2 

i  , 

/  i  \ 

2  V  Vik 

■  S  -ik’^ik) 

k=l 

\k=l  1 

Hi  -  1 


where  =  log  Substituting  Eqs.  (49)  and  (50)  into  Eq,  (48)  gives 

an  estimate  for  the  probability  density  function,  p(x^) . 

Various  statistical  inferences  can  now  be  made  about  the  vibration 

environment  from  which  the  data  were  taken.  First,  the  1 00  a  percentage 

point  of  the  spatial  distribution  function  for  the  vibration  is  given  by  x  , 

a 

which  satisfies 


p(Xj)  dx^  =1-0- 


(51) 


It  follows  that 


-/3(x^)  =  P(x.  <  x^)  =  P(log  X.  <  log  x^) 


(52) 


where 'P(z)  is  the  standardized  normal  distribution  function.  Let  z  -le- 

a 

note  the  100  o  percentage  point  of  the  standardized  normal  distribution.  Then 
it  follows  that 


log  X  -  p. 
a  1 

- - A -  =  Z 

(T.  a 

1 


(53) 
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Hence,  the  100  a  percentage  point  of  p(x^)  is  given  by 


A  A 

z  O-.+p. 

ail 

X  =  e 
a 


(54) 


Because  of  the  functional  dependence  of  a  lognormal  variable  on 
a  normal  random  variable,  various  pertinent  problems  of  statistical 
inference  for  a  lognormal  probability  model  may  be  solved  ut^ing  the 
properties  of  a  normal  distribution.  For  example,  100(1  -  a)  % 
tolerance  limits  with  probability  y  for  the  lognormal  variable  x^  can 
be  constructed  as  follows. 


(55) 


where  K  is  the  same  constant  appropriate  for  the  corresponding  toler¬ 
ance  limits  for  the  normal  variable  y^ ;  that  is. 


Eq.  (55)  is  deduced  from  the  following  relation 
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7.  2  DEVELOPMENT  OF  TESTING  PROCEDURES 

The  degree  to  which  the  distribution  of  vibration  data  fits  an  assumed 

2 

distribution  function  model  can  be  tested  by  the  classical  y  goodness - 
of-fit  test,  as  follows.  Assume  the  mean  square  values  of  the  data  in  the 
ith  frequency  interval  are  divided  into  m^  mutually  exclusive  intervals, 
where  the  number  of  mean  square  values  in  each  of  the  intervals  is 
n^j^ ;  k  =  1,  2,  3,  .  .  . ,  m^  .  Let  the  corresponding  probability  of  a  mean 
square  value  being  in  each  of  the  intervals  be  P.  (G  ,  9_,  . . . ,  Q  )  ; 

XiC  L  U 

k  =  1,  2,  3,  .  .  .  ,  m^ ,  where  P.^  is  determined  from  an  empirical  distribu¬ 
tion  function  with  a  set  of  parameters.  If  the  true  values  of 
are  unknown,  the  statistic  for  a  test  is  given  by 


^  "^i  fn.,  -  n.  P.,{G, . G 

L  ik  1  ik  1 _ u^ 

~  A  A 

k=l  n.  P.,  (G,.  ....  G  ) 

1  ik  1  u 


2  2 

where  is  asymptotically  distributed  as  x  with  m^  -  u  -  1  degrees-of- 
freedom.  Note  that  u  number  of  degrees-of-freedom  are  lost  because 
u  parameters  are  estimated. 
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The  adequacy  of  the  fit  to  vibration  data  provided  by  the  log- 

2 

normal  distribution  defined  in  Eq.  (48)  can  be  tested  by  the  x  test 
given  in  Eq.  (57).  For  the  ith  frequency  interval,  the  statistic  used 
is 


A 

where  is  an  estimate  given  by 


A 

P(Xj)  dXj 


f 


(58) 


(59) 


The  term  p(x,)  is  an  estimate  of  p.(x.)  as  described  in  Section  7.  1.2, 

and  a  ,  b.,  are  the  lower  and  upper  limits  of  the  kth  mean  square 
Ik  Ik  2 

value  interval.  The  statistic  X  given  by  Eq.  (58)  has  an  asymptotic 
chi-square  distribution  with  m  -  3  degrees -of-freedom.  Evaluation 
of  Eq.  (58)  can  be  performed  conveniently  as  follows. 


c 

I 


^ik  '  ^  <»ik  ^  ’‘i  ^  •’ik’  =  P  ^k  ‘  n  ^  V 


where "^(z)  denotes  the  standardized  normal  distribution  function,  and 

A  A 

\x,  and  O’,  are  given  by  Eqs.  (49)  and  (50). 

^  ^  2 

The  X  test  described  above  can  be  implemented  to  perform  an 
overall  test  for  the  fit  of  c  sets  of  empirical  distributions.  For  ex¬ 
ample,  it  may  be  desired  to  test  the  overall  goodness-of-fit  of  the  log¬ 
normal  probability  model  for  all  frequency  intervals  in  a  particular  zone. 
The  desired  statistic  is  given  by 


(61) 


The  variable  has  a  x  distribution  with  degrees -of-freedom  given 
by 

c 

d.  f.  =  y  [  m.  -  2c  -  c 

rr  1 


(62) 


7.  3  APPLICATION  TO  AFFDL  VIBRATION  DATA 

The  procedures  developed  in  Section  7.  1  are  now  applied  to  test 
a  lognormality  hypothesis  for  the  spatial  distribution  of  aircraft  vibration 
using  AFFDL  data.  The  structural  zone  selected  for  the  test  is  Zone  02 
(middle  half  of  fuselage)  of  Aircraft  Group  5  (Century  series  fighters). 

As  for  the  homogeneity  tests  in  Section  6.  2,  data  for  aircraft  Code  No.  63 
is  omitted  to  assure  consistency  in  the  measurements.  This  particular 
zone  was  selected  because  (1)  it  provides  a  relatively  large  sample  of 
data,  and  (2)  it  represents  a  basic  structural  region  and  aircraft  group 
of  considerable  interest.  Since  only  spatial  variations  in  the  vibration 
levels  are  of  concern,  data  for  a  single  flight  condition  are  used  to  mini¬ 
mize  variations  due  to  different  excitation  characteristics.  The  flight 
condition  selected  for  the  test  is  Condition  07  (normal  cruise). 

The  flight  vibration  levels  during  normal  cruise  for  Zone  02  of 
Aircraft  Group  5  are  summarized  in  Table  17.  A  total  of  4273  data 
values  are  given.  Before  proceeding,  it  should  be  noted  that  a  critical 
assumption  must  be  made  concerning  these  data.  Specifically,  if  the 
data  in  Table  17  are  to  represent  the  spatial  distribution  of  vibration 
levels  in  Zone  02,  it  must  be  assumed  that  each  value  represents  a 
sample  of  the  vibration  at  a  randomly  selected  location  in  that  zone.  In 
reality,  this  assumption  is  not  totally  valid  for  two  reasons.  First,  the 
selection  of  transducer  locations  for  the  measurements  was  unquestionably 
influenced  by  such  nonrandom  practical  considerations  as  physical  accessi¬ 
bility  and  proximity  to  equipment  locations  of  interest.  Second,  some  of 
the  data  values  represent  repeated  measurements  at  the  same  point  on 
the  same  aircraft.  In  spite  of  these  possible  biasing  effects,  however, 
the  randomness  assumption  is  considered  sufficiently  acceptable  to  pro¬ 
duce  adequate  results. 
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The  data  in  Table  17  are  presented  in  terms  of  double  amplitude 
histograms  for  each  of  20  frequency  increments.  To  be  consistent 
with  the  development  in  Section  7.  1,  those  double  amplitude  intervals 
which  include  data  are  converted  to  mean  square  value  intervals  in 
Table  18.  This  step  is  not  actually  necessary  since  a  lognormal  distri¬ 
bution  for  double  amplitude  values  would  imply  a  lognormal  distribution 
for  moan  square  values  as  well.  Nevertheless,  the  step  is  made  to 
fully  illustrate  the  testing  procedure  as  it  should  be  applied  to  later  data 
which  include  an  accurate  representation  of  random  contributions. 

Such  data  would  normally  be  presented  in  terms  of  power  spectra 
(normalized  mean  square  values  versus  frequency).  Also  included  in 
Table  18  are  the  natural  logarithms  of  the  mean  square  value  interval 
limits  (a^,  b^^)  and  midpoints  (y^^) . 

The  mean  values  and  standard  deviations,  |x^  and  cr^ ,  for  the 
logarithms  of  the  data  in  each  frequency  increment  are  given  in  Table  19. 

These  values  are  calculated  using  Eqs.  (49)  and  (50).  With  this  informa- 

2 

tion,  the  values  of  are  calculated  for  each  frequency  increment  using 

Eq.  (58).  An  illustration  of  these  calculations  for  the  i  =  13th  frequency 

increment  (80-99  cps)  is  outlined  in  Table  20.  The  final  results  for  the 
2 

calculations  in  all  frequency  increments  are  summarized  in  Table  21. 
Note  that  the  calculations  were  performed  with  the  end  intervals  pooled 
as  required  to  make  the  expected  number  of  samples,  nP  ,  at  least  3. 

K 

Also  note  that  five  frequency  increments  (i  =  3,  4,  5,  6  and  9)  are  omitted 

from  the  final  results  because  the  sample  size  for  these  increments  is 

2 

not  considered  sufficient  for  the  x  test.  Actually,  the  sample  size  for 

many  of  the  frequency  increments  is  less  than  would  normally  be  desired 
2 

for  a  X  test  (at  least  200) .  Furthermore,  the  number  of  mean  square 
value  intervals,  m^^ ,  for  the  various  frequency  increments  is  not  always 
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Table  18.  Mean  Square  Value  Intervals  for  Lognormality  Test 


log  [b  X  10^°] 


log  [b  X  10^°] 


51.  0  X  10 
32.  0  X  lO"^ 

20.  2  X  10"^ 

n.  5  X  10"^ 
8.  14  X  lO' 
5.  1 2  X  10 
3.  21  X  10“ 

2.  02  X  10“ 

1. 28  X  lO'' 
80.  0  X  lO’^ 
51.  0  X  lO"^ 
32.  0  X  lo'^ 
20.  2  X  lO'^ 
12.  5  X  10‘^ 
8.  14  X  10‘^ 
5.  12  X  lO'^ 
3.  21  X  10'^ 
2.  02  X  lO"^ 

1.  28  X  10'^ 


17.  75 

17.  48 

80.  0  X  lO"®  8.  98 

17.  22 

17.  01 

51.0x10'®  8.  54 

16.  81 

16.  58 

32.0x10'®  8.08 

16.  35 

16.  12 

20.2x10'®  7.60 

15.  90 

15.  67 

12.  5  X  10'®  7.  12 

15.  44 

15.  21 

8.  14  X  lO'®  6.  70 

14.98 

14.  75 

5.  12x10  ®  6.  24 

14.  52 

14.  29 

3.  25  X  lO'®  5.  78 

14.  06 

13.  82 

2.  00  X  lO’®  5.  30 

13.  58 

13.  36 

1.  28  X  lO"®  4.  85 

13.  14 

12.  90 

78.1x10"^°  4.36 

12.  66 

12.  43 

50.  0  X  lO’^^  3.  80 

12.  20 

11.  97 

21.  1  X  10'^°  3.  05 

11.  74 

1 1 .  r)2 

12.  5  X  10’^°  2.  52 

11.  30 

11.  07 

8.  0  X  10'^°  2.  08 

10.84 

10.  61 

4.  5  X  10  1.  50 

10.  38 

10.  15 

3,  12x10^®  1.  14 

9.  92 

9.  69 

2.  00  X  10  0.  69 

9.  46 

9.  22 

1.12x10'^°  0.11 

'bj^  =  upper  limit  of  kth  interval;  a^  =  lower  limit  of  kth  interval 

;o::  log  X  lOlO]  +  log  [b^_^  X  10 10] 

2  ~  midpoint  of  kth  log  interval 
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Table  ZO.  Illustration  of  Calculations  for  Lognormality  Test 

Calculations  for  frequency  increment  i  =  13  (80-99  cps) 

Sample  size  n  =  845 


°<“k' 


< 

!  * 

'"k 

G(l^) 

2.  61 

00 

0. 9955 

2.  25 

2.  61 

0.  9878 

1.  91 

2.  25 

0.  9719 

1.  56 

1.  91 

0. 9406 

1.  20 

1.  56 

0.  8849 

0.  84 

1,  20 

0, 7995 

0  52 

0.  84 

0. 6985 

0.  17 

0.  52 

0.  5675 

-0.  18 

0.  17 

0.  4286 

-0.  55 

-0.  18 

0.  2912 

-0.  88 

-0.  55 

0.  1894 

-1.  26 

-0.  88 

0.  1038 

-00 

-1.  26 

0 

log  [a  X  10^®]  -  |i 


2 

*»  "k  -  “Pk 


1.  0000 
0.  9955 
0.  9878 
0.  9719 
0.  9406 
0.  8849 
0.  7995 


0.  5675 
0.  4286 
0.  2912 
0.  1894 
0.  1038 


0.  0045 
0. 0077 
0. 0159 
0. 0313 
0.  0557 
0.  0854 
0.  1010 
0. 1310 
0.  1389 
0.  1374 
0. 1018 
0.  0856 
0.  1038 


log  [b,  X  10^®] 


iiwy' 


Table  21.  Summary  of  Results  for  Lognormality  Test 


optimal  for  the  corresponding  sample  size.  Such  problems  are  un¬ 
avoidable,  however,  since  the  data  were  not  collected  and  reduced  with 
the  intent  of  performing  a  test. 

2 

Referring  to  Table  21,  the  value  of  ,  as  defined  in  Eq.  (61), 

is  1091. 8.  The  total  number  of  mean  square  value  intervals  is  m  =  146  , 

and  the  number  of  data  sets  is  c  =  15.  Hence,  from  Eq.  (62),  if  the 

2 

data  in  Table  17  has  a  lognormal  distribution,  then  X  is  distributed 

2  ^ 

approximately  like  y  with  101  degrees-of-freedom.  From  any  table 

2  2 
of  y  values,  it  is  seen  that  the  1%  point  of  y  with  101  degrees -of- 

2  2 

freedom  is  equal  to  136.  Hence,  since  X^  =  1091.  8  >  Xq  o1‘ 101  ~ 
a  lognormal  hypothesis  for  the  data  in  Table  17  is  rejected  at  the  1% 
level  of  significance. 

Technically,  the  foregoing  results  apply  only  to  one  structural 
zone  of  one  aircraft  group.  A  visual  inspection  of  histograms  (similar 
to  Table  17)  for  data  from  various  other  zones  and  aircraft  groups, 
however,  indicates  that  these  results  maybe  generalir.ed.  Specifically, 
based  upon  past  AFFDL  vibration  data,  the  spatial  distribution  for  air¬ 
craft  vibration  levels  is  not  lognormal  in  character.  This  fact  is  best 
illustrated  by  comparing  the  empirical  density  function  for  measured 
data  to  the  lognormal  density  function,  as  is  done  in  Figure  17.  The 
empirical  distribution  in  Figure  17  is  calculated  from  data  in  Table  20 
for  the  13th  frequency  increment  (80-99  cps)  of  Zone  02  for  Aircraft 
Group  5.  The  results  are  typical  of  those  obtained  in  other  frequency 
increments.  In  particular,  there  is  a  definite  tendency  for  the  distri¬ 
bution  of  log  values  to  be  skewed  so  that  the  higher  (but  not  extreme) 
mean  square  values  occur  with  somewhat  greater  probability  than  would 
be  predicted  using  the  lognormal  assumption. 


I  I  I  I  1 

-2<r  -lo'  M  l<r  Z  (T 

LOG  MEAN  SQUARE  VALUE 

A  =  7.4Ii  a  =  1.316,  n=845 

Figure  17.  Illustration  of  Spatial  Distribution  for  Vibration  Data 
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The  failure  of  the  data  to  pass  a  lognormality  test  does  not  rule 
out  the  possibility  that  a  lognormal  assumption  might  still  be  acceptable 
for  selecting  test  levels  and  design  criteria.  This  possibility  is  illus¬ 
trated  by  comparing  the  vibration  level  for  some  small  percentage  point 
of  the  actual  distribution  function  to  the  level  for  that  same  percentage 
point  of  a  fitted  lognormal  distribution  function.  As  before,  the  data 
for  frequency  increment  13  of  Zone  02  for  Aircraft  Group  5  are  used. 
Assume  a  prediction  at  the  2%  point  (a  98%  upper  limit)  is  desired. 

From  Table  17  or  20,  it  is  seen  that  the  2%  point  for  the  actual  data 
occurs  at  log  [a  x  10^^]  =  2.  61  ^  +  p  =  9.  46,  which  corresponds  to  a 
double  amplitude  of  0.  0032  inch.  From  Eq.  (54),  the  2%  point  for  a 
lognormal  distribution  fitted  to  the  mean  square  values  of  the  data  occurs 
at  log  [a  X  10^^]  =  2.  056  ?  +  p  =  8.  73,  which  corresponds  to  a  double 
amplitude  of  0.  0023  inch.  It  is  interesting  to  compare  these  results  with 
those  obtained  by  assuming  the  data  have  a  simple  normal  distribution. 
First  consider  the  case  where  the  mean  square  values  are  assumed  to  be 
normally  distributed.  From  Eqs.  (43)  and  (44),  the  mean  value  and 

A  -6  2 

standard  deviation  for  the  mean  square  values  are  p  =  0.  099  x  10  inch 
and  a  =  0.  212  X  10  ^  inch^ .  The  2%  point  is  then  0.  535  x  10  which 
corresponds  to  a  double  amplitude  of  0.  0021  inches.  Hence,  a  normal 
assumption  for  mean  square  values  provides  a  poorer  result  than  the  log¬ 
normal  assumption.  Now  consider  the  case  where  the  original  double 
amplitude  values  are  assumed  to  be  normally  distributed.  From  Table  17, 
the  mean  value  and  standard  deviations  for  the  double  amplitudes  are 
X  =  0.  00082  inch  and  s  =  0.  00081  .  The  2%  point  for  a  normal  distribution 
fitted  to  the  double  amplitudes  is  then  0.  0025  inch.  This  result  is  slightly 
better  than  achieved  using  the  lognormal  assumption. 
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The  results  of  the  spatial  distribution  study  to  this  point  are  incon¬ 
clusive.  For  the  problem  of  selecting  conservative  test  levels  and  de¬ 
sign  criteria,  it  appears  that  the  lognormal  distribution  for  mean  square 
values  may  provide  results  which  are  no  better  than  those  obtained  us¬ 
ing  a  straightforward  normal  distribution  assumption  for  peak  or  rms 
values.  Perhaps  a  more  accurate  model  for  the  spatial  distribution  of 
aircraft  vibration  should  be  sought  out  or  developed.  Further  studies  are 
required  to  resolve  these  issues.  It  is  not  considered  appropriate,  how¬ 
ever,  to  pursue  such  studies  at  this  time  because  of  the  limited  quality 
of  the  available  data.  Any  conclusions  that  might  result  from  further 
studies  of  the  currently  available  AFFDL  vibration  data  would  be  subject 
to  later  confirmation  and/or  modification  based  upon  similar  studies  of 
future  data  which  properly  represent  the  random  contributions  in  the 
vibration  environment.  Hence,  it  is  suggested  that  AFFDL  defer  studies 
of  this  problem  until  such  future  data  are  available  in  quantity. 
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8.  PREDICTION  MODEL  STUDIES 


Attention  is  now  directed  to  the  primary  purpose  of  this  effort: 
namely,  the  development  of  a  general  extrapolation  model  for  predictin 
aircraft  vibration  environments.  As  a  first  step,  it  is  appropriate  to 
review  the  information  needed  to  develop  such  a  model,  and  to  assess 
the  progress  made  thus  far  towards  acquiring  this  information. 

From  the  discussions  in  Section  2,  the  development  of  an  aircraft 
vibration  prediction  model  of  the  general  extrapolation  type  requires 
four  preliminary  steps,  as  follows. 

1.  Approximate  relationships  between  fundamental  source  param 
eters  and  the  resulting  vibration  environment  must  be  developed.  This 
step  is  required  to  define  the  structural,  engine,  and  flight  parameters 
which  should  be  measured  and  recorded  with  the  vibration  data,  and 

to  provide  an  initial  model  for  regression  studies. 

2.  A  large  collection  of  accurate  vibration  data  for  representative 
locations  in  various  aircraft  must  be  obtained  and  identified  with  the 
fundamental  parameters  defined  in  step  1.  This  step  provides  the  basic 
data  required  for  the  regression  studies. 

3.  Each  aircraft  type  of  interest  must  be  divided  into  appropriat«; 
structural  zones,  where  each  zone  will  be  the  basis  for  a  single  predic¬ 
tion  model.  This  step  is  necessary  to  reduce  the  complexity  and  re¬ 
dundancy  of  the  resulting  prediction  model. 
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4.  The  general  form  of  the  spatial  distribution  function  for  the 
vibration  within  each  zone  must  be  approximated.  This  step  is  needed 
to  provide  a  firm  statistical  basis  for  predicting  conservative  limits 
for  test  specifications  and  design  criteria. 

Step  1  is  dealt  with  in  Section  3.  Vibration  sources  are  identified 

and  relationships  are  developed  for  all  aircraft  groups.  However,  since 

the  design  of  propeller  type  aircraft  (Groups  1  and  2)  is  unlikely  in  the 

future,  these  aircraft  will  not  be  considered  further.  Attention  will  be 

restricted  to  jet  powered  aircraft,  as  represented  by  Groups  3,  5,  and 

10.  The  results  in  Section  3  indicate  that,  for  jet  powered  aircraft, 

the  critical  source  parameters  which  should  be  recorded  and  identified 

with  measured  vibration  data  are  p  (air  density),  A  (jet  exit  area), 

V  (exhaust  gas  velocity),  c  (speed  of  sound),  c  (local  speed  of  sound 
€  € 

in  jet  exhaust),  q  (dynamic  pressure),  w  (average  surface  weight  density), 
and  the  rpm  for  all  rotors,  gears,  and  other  rotating  equipment.  The 
suggestions  in  Section  3  for  analytical  relationships  between  these  vari¬ 
ous  source  parameters  and  the  resulting  vibration  environments  should 
provide  an  adequate  basis  for  initial  regression  studies.  It  should  be 
emphasized,  however,  that  the  relationships  might  be  modified  an  a  trial 
and  error  basis  as  required  to  obtain  good  correlation  coefficients  in  the 
regression  analyses. 

Concerning  step  2,  AFFDL  has  collected  considerable  quantities 
of  vibration  data  on  many  different  aircraft,  as  summarized  in  Section  4. 
The  data  are  well  documented  and  properly  identified  with  the  pertinent 
source  parameters  indicated  in  step  1.  In  most  cases,  however,  the 
quality  of  the  data  is  not  sufficient  for  use  in  the  development  of  a  vibra¬ 
tion  prediction  model.  The  reasons  for  this  are  detailed  in  Section  5. 
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Recent  revisions  in  AFFDL  data  acquisition  and  analysis  procedures 
should  provide  suitable  data  for  this  application  in  the  future,  but  at 
this  time,  such  data  are  not  available  in  quantity.  Data  of  reasonable 
quality  are  available  for  one  aircraft  only;  namely,  an  RF-4C  (Code 
No.  63).  Hence,  the  data  for  this  aircraft  will  be  used  to  illustrate  the 
desired  developments. 

Step  3  is  pursued  in  Section  6.  Proper  zoning  of  aircraft  is  in-  \ 

vestigated  based  upon  homogeneity  studies  of  past  AFFDL  vibration  data. 
Because  of  the  limited  quality  of  the  data,  firm  conclusions  cannot  be 
drawn.  Tentative  conclusions,  however,  are  as  follows.  Separate  pre¬ 
diction  models  are  not  required  for  the  vibration  along  the  three 
orthogonal  axes,  or  for  the  vibration  of  different  aircraft  within  a  group. 

On  the  other  hand,  separate  prediction  models  should  be  developed  for 
the  vibration  in  each  structural  zone  designated  by  AFFDL,  and  for  the 
vibration  of  each  aircraft  group.  To  facilitate  further  developments, 
it  will  be  assumed  that  the  above  conclusions  are  valid. 

Step  4  is  considered  in  Section  7.  A  lognormal  hypothesis  for  the 
spatial  distribution  of  mean  square  vibration  levels  over  a  structural  zone 
is  tested  using  past  AFFDL  vibration  data.  As  in  step  3,  firm  conclusions 
are  not  possible  because  of  the  data  quality.  The  tentative  conclusion, 
however,  is  that  the  lognormal  distribution  does  not  provide  an  acceptable 
description  for  the  spatial  distribution  of  aircraft  vibration.  In  fact,  it 
appears  that  a  simple  normal  distribution  for  peak  or  rms  values  might 
estimate  the  vibration  levels  at  small  percentage  points  with  as  much 
accuracy  as  a  lognormal  distribution  for  mean  square  values.  Neither 
distribution  function,  however,  produces  desirable  results.  To  facilitate 
further  developments,  the  normal  distribution  will  be  assumed  for  the 
spatial  distribution  of  peak  or  rms  vibration  levels  within  a  given  struc¬ 
tural  zone. 
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8.  1  MODEL  FORMULATION  FOR 

JET  BOMBERS  AND  FIGHTERS 

Separate  prediction  models  probably  will  be  required  for  jet 
bombers  (Group  3)  and  jet  fighters  (Group  5).  The  general  form  of  the 
models,  however,  will  be  similar  since  the  primary  vibration  sources 
for  these  two  aircraft  groups  are  the  same.  Hence,  they  will  be  con¬ 
sidered  together. 

From  Table  6  in  Section  3,  there  are  five  sources  of  vibration 
for  jet  bombers  and  fighters  which  are  considered  sufficiently  significant 
to  include  in  a  prediction  model.  In  review,  these  sources  are  as 
follows. 

R-1  —  jet  acoustic  noise 

R-2  —  aerodynamic  boundary  layer  noise 

R-6  —  wake  turbulence  caused  by  aerodynamic  devices 

(lift  devices,  drag  devices,  and  spoilers) 

D-3  —  jet  engine  shaft  rotation 

D-10  —  gunfire 

From  the  viewpoint  of  test  specifications  and  design  criteria,  the 
ultimate  requirement  is  to  predict  the  most  severe  vibration  levels  which 
occur  during  the  aircraft  service  life.  This  requirement  permits  the 
prediction  model  to  be  segmented  and  simplified.  For  example,  the 
maximum  vibration  induced  by  jet  noise  clearly  occurs  during  takeoff 
when  the  contribution  of  boundary  layer  noise  is  negligible.  Correspond¬ 
ing,  the  meiximum  vibration  induced  by  boundary  layer  noise  occurs  dur¬ 
ing  high  speed  flight  when  the  effects  of  jet  noise  are  small.  The  contri¬ 
butions  of  jet  noise  can  also  be  neglected  when  wake  turbulence  is  pre¬ 
dicted  since  the  engine  power  is  usually  low  when  the  various  aerodynamic 
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devices  are  deployed.  Guns  are  not  used  at  takeoff.  Hence,  the  worst 
case  for  gunfire  induced  vibration  is  probably  at  high  speeds  when  jet 
noise  contributions  are  small.  Based  upon  these  considerations,  it  is 
suggested  that  the  power  spectrum  for  jet  bomber  and  fighter  vibration 
environments  be  predicted  in  four  segments,  as  follows. 


1 .  For  Takeoff 

°<«to  = 

2.  For  High  Speed  Flight 

=  °<«r2  "  °<«d3 

3.  For  Deployment  of  Aerodynamic  Devices 

°'«ad  =  "  °<‘>d3  ^  °<«r6 

4.  For  Gunfire 

G(f)gj  =  *  °<«d3  ^  °<«dlO 


Note  that  a  separate  set  of  predictions  will  apply  to  each  structural  zone 
of  each  aircraft  group. 

Suggested  techniques  for  determining  each  of  the  contributing  power 
spectra  terms  in  Eq.  (63)  are  now  outlined. 


8.  1.  1  Jet  Noise  Contribution 

From  Eqs.  (26)  and  (28),  the  initial  model  for  the  vibration  in¬ 
duced  in  a  specific  zone  by  jet  noise  is 

(64) 

where 

A 

G  (fl)  j  =  average  power  spectral  density  for  the  acceleration 
response  versus  dimensionless  frequency 
A 

A(f2)  j  =  regression  coefficient  (estimated  weight)  versus 
dimensionless  frequency 

n  5  2 

X  ^  =  pAV”/c  w  (independent  variable) 

=  fdc  /V  c  (dimensionless  frequency) 
e  e 

p  =  atmospheric  density 
A  =  jet  exit  area 

=  exhaust  gas  velocity 

w  =  surface  weight  density  of  structure 
c  =  ambient  speed  of  sound 
c^  =  local  speed  of  sound  in  jet 

d  =  A/tt 

n  =  constant  between  6  and  8 
f  =  frequency  in  cps 

Equation  (64)  is  written  in  terms  of  dimensionless  frequency  to  account 
for  frequency  shifts  of  the  jet  noise  spectrum  due  to  variations  in  the 
parameters  which  define  the  dimensionless  frequency.  A  specific  pre¬ 
diction  maybe  readily  converted  to  power  spectral  density  versus  fre¬ 
quency  in  cps  using  the  ^^^ii^ition  for  . 


The  weight,  A(n)^^  ,  in  Eq.  (64)  is  estimated  by  a  conventional 
regression  analysis  of  "appropriate  data,  "  as  discussed  in  Section  2.  2 
(the  definition  of  "appropriate  data"  is  discussed  later).  For  the  problem 
at  hand,  only  one  independent  variable  is  involved.  Hence,  the  esti¬ 
mation  is  accomplished  using  Eq.  (10).  With  the  estimation  of  the  weight, 
Eq.  (64)  can  then  be  applied  to  predict  the  average  power  spectrum  for 
the  jet  noise  induced  vibration  in  a  similar  structural  zone  of  a  similar 
future  aircraft. 

The  final  goal,  of  course,  is  to  predict  some  conservative  upper 
limit  for  the  vibration.  The  regression  program  will  calculate  any 
desired  prediction  limit  for  the  data,  given  a  distribution  function  for 

the  data  about  the  regression  line,  as  detailed  in  Section  8.  3  to  follow. 

2 

Note  that,  because  of  the  w  term,  the  distribution  function  of  concern 
here  is  not  the  spatial  distribution  discussed  in  Section  7.  It  is  a  more 
complex  function  dependent  upon  the  spatial  distribution  of  both  the 
vibration  levels  and  the  surface  mass  density.  For  convenience,  the 
function  will  be  assumed  normal. 

Equation  (64)  provides  an  acceptable  model  for  regression  analysis, 
assuming  the  relationships  in  the  independent  variable  are  valid. 

If  there  is  some  uncertainty  about  any  terms  forming  ,  however,  the 
model  could  be  expanded  to  evaluate  these  terms  separately.  For  ex¬ 
ample,  there  is  uncertainty  concerning  the  proper  value  for  n  in  the 
v”  term.  Hence,  it  might  be  interesting  to  study  the  model  written  in 
the  following  form 


log  j  =  log  j  +  n  log  ^  V^J 


163 


Equation  (65)  constitutes  a  regression  equation  of  the  form  y  =  a^  +  a^ 

+  where  a^  =  log  j  ,  a^  =  1  and  a^  =  n.  Hence., the  re¬ 
gression  analysis  will  calculate  the  values  for  both  log  ^ 

which  provide  the  best  fit  tc  the  data.  The  procedure  could  be  extended 
to  solve  for  the  power  of  every  term  in  ,  but  such  a  procedure  is 
not  suggested.  As  the  number  of  terms  in  the  regression  model  is  in¬ 
creased,  the  variance  for  the  calculated  coefficients  is  also  increased. 

It  is  desirable  to  restrict  the  model  as  much  as  can  be  jusrified  by  known 
facts. 

The  next  problem  is  the  selection  of  "appropriate  data"  for  the  re¬ 
gression  analysis.  The  data  used  should  represent  as  closely  as  possible 
the  vibration  due  only  to  jet  acoustic  noise.  Such  data  may  he  obtained 
from  sample  records  of  the  vibration  during  takeoff,  preferably  at  the 
start  of  takeoff  after  the  engine  has  reached  full  power  but  before  sig¬ 
nificant  speed  has  been  achieved.  The  sample  records  should  be  re¬ 
duced  to  power  spectra  with  all  periodic  components  removed.  Since 
the  vibration  produced  by  jet  noise  is  random,  any  periodic  components 
that  might  be  present  must  represent  other  sources  (principly  engine 
shaft  rotation).  See  Section  8.  1.  6  for  further  discussion  of  data  reduction 
procedures. 

The  final  problem  is  the  determination  of  surface  weight  densities 
(w)  for  the  measurement  locations.  For  the  basic  structural  zones,  mea¬ 
surements  will  generally  be  made  on  frames,  stringers,  truss  sections, 
or  other  structures  with  clearly  defined  widths.  The  surface  weight  density 
for  such  structures  should  be  calculated  as  follows.  Compute  the  weight 
per  unit  length  of  the  structure  and  divide  by  the  width.  If  the  structure 
is  attached  to  a  skin  section,  add  in  the  surface  weight  density  of  the  skin. 


If  a  component  is  attached  to  the  structure  at  the  measurement  location, 
add  in  an  effective  surface  weight  density  of  the  component  given  by 


w 

e 


w 

e 


N  A 
a  s 


(66) 


where 


w^  =  total  weight  of  component 

=  number  of  attachment  points 

A  =  area  of  the  structure  in  the  general  region 
where  the  component  is  attached 


8-  1.  2  Aerodynamic  Noise  Contribution 

From  Eqs.  (26)  and  (30),  the  initial  model  for  the  vibration  in¬ 
duced  in  a  specific  zone  by  aerodynamic  noise  is 


where 


(67) 


r2 


A(f) 


r2 


average  power  spectral  density  for  the  acceleration 
response  versus  frequency 

regression  coefficient  (estimated  weight)  versus 
frequency 

q^/w^  =  V^/4w^  (independent  variable) 


p  =  atmospheric  density 
V  =  aircraft  velocity 

w  =  surface  weight  density  of  the  structure 
f  =  frequency  in  cps 
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The  procedures  for  calculating  the  coefficient  A(f)^2>  arriving  at 
a  conservative  prediction  limit  for  the  vibration,  and  for  determining 
proper  values  for  the  surface  weight  density  w,  are  the  same  as  dis¬ 
cussed  in  Section  8.  1.  1. 

Now  consider  the  problem  of  selecting  appropriate  data  for  the 
regression  analysis.  The  data  should  rej/resent  as  closely  as  possible 
the  vibration  due  only  to  boundary  layer  turbulence.  The  best  source 
of  data  would  be  from  sample  records  collected  during  clean  dives  with 
low  engine  power.  However,  sample  records  collected  during  any 
clean  cruise  condition  should  be  acceptable.  Care  must  be  exercised 
to  exclude  all  data  obtained  during  those  flight  conditions  when  lift,  drag, 
or  spoiler  devices  are  deployed,  or  during  climb,  acceleration,  or  other 
flight  conditions  where  the  engine  power  is  high  compared  to  the  resulting 
airspeed.  The  sample  records  should  be  reduced  to  power  spectra  with 
all  periodic  components  removed.  As  for  jet  noise,  the  vibration  in¬ 
duced  by  boundary  layer  noise  is  random.  Thus,  any  periodic  components 
which  might  be  in  the  data  must  represent  other  sources.  See  Section 
8.  1.  6  for  further  discussions  of  data  reduction  procedures. 

8.1.3  Wake  Turbulence  Contribution 

The  most  straightforward  approach  tc  establishing  the  conbribution 
of  wake  turbulence  is  to  compare  the  vibration  measured  at  similar  points 
under  similar  flight  conditions  with  and  without  various  lift,  drag,  and 
spoiler  devices  deployed.  By  reducing  the  data  in  terms  of  power  spectra 
with  all  periodicities  removed,  the  contribution  of  wake  turbulence  due 
to  the  deployment  of  a  given  device  would  be 


I 


t 


where 


G  (f)  =  measured  power  spectrum  with  wake  turbulence 

producing  device  deployed 

G  (f)  =  measured  power  spectrum  without  wake  turbulence 
producing  device  deployed 


The  above  procedure  should  be  repeated  using  data  collected  for 
different  values  of  q  to  determine  if  displays  a  q  dependence. 

If  it  does,  a  conservative  prediction  limit  for  G  (f)  .  should  be  estab- 
lished  by  a  regression  analysis  of  the  following  model. 


log 


’a  r a  r  * 


(69) 


If  it  is  found  that  B  =  2,  then  the  q  dependence  is  the  same  as  suggested 
in  Eq.  (b7).  Note  that  a  weight  term  (w)  is  not  suggested  in  Eq.  {b9). 
Since  wake  turbulence  is  predominately  a  low  frequency  phenomenon, 
the  structural  response  is  primarily  controlled  by  stiffness  rather  than 
mass. 

If  a  q  dependence  is  not  detected  in  the  calculated  values  for 
G^(f)^^  ,  then  a  conservative  prediction  limit  should  be  established  as 
follows.  For  the  measurements  in  each  structural  zone,  pool  together 
the  values  of  G  (f)  ,  determined  for  all  flight  conditions.  Calculate  the 
mean  value  and  standard  deviation  of  the  spectra  versus  frequency  for 
each  zone.  Fit  an  appropriate  distribution  function  (unknown  at  this  time) 
using  the  calculated  mean  and  standard  deviation  at  each  frequency. 
Determine  a  prediction  limit  by  selecting  the  spectral  density  level  at 
any  desired  percentage  point  of  the  fitted  distribution.  Further  dis¬ 
cussions  of  the  statistical  aspects  of  selecting  prediction  limits  are 
presented  in  Section  8.  5. 


b: 


It  should  be  mentioned  that  an  effort  was  made  to  evaluate  the 
contribution  of  wake  turbulence  using  currently  available  AFFDL.  data. 
Histograms  for  the  vibration  data  measured  in  various  zones  of  jet 
bombers  and  fighters  during  clean  flight  conditions  were  compared  to 
similar  data  measured  during  those  flight  conditions  when  various  wake 
turbulence  producing  devices  were  deployed.  In  many  cases,  the  com¬ 
parisons  did  not  reveal  any  significant  difference  in  the  average  vibra¬ 
tion  levels  for  the  two  situations.  This  failure  to  detect  the  contributions 
of  wake  turbulence  is  probably  due  to  the  limited  quality  of  the  available 
data.  On  the  other  hand,  wake  turbulence  may  not  be  as  significant  as 
originally  believed.  These  issues  must  be  resolved  by  applying  Eq.  (68) 
to  future  data  obtained  in  the  manner  specified. 

8.  1.  4  Jet  Engine  Shaft  Rotation  Contribution 

\'ibration  induced  by  jet  engine  shaft  rotation  is  one  contribu¬ 
tion  that  is  fully  represented  and  described  by  past  AFFDL  vibration 
data.  Jet  engine  shaft  vibration  is  periodic  with  principal  frequency  com¬ 
ponents  below  500  cps.  Hence,  past  Ai  FDL  data  reduction  techniques, 
as  summarized  in  Section  5,  tended  to  emphasize  the  vibration  from  this 
source.  This  fact  is  confirmed  by  an  inspection  of  the  shot  gun  data 
plots  in  [19]  -  [25]  .  It  is  seen  in  there  references  that  the  vast  majority 
of  data  points  are  clustered  at  the  frequencies  of  shaft  rotation  and  their 
harmonics. 

It  was  previously  concluded  in  Section  3.  3.  6  that  attempts  to  corre¬ 
late  engine  shaft  induced  vibration  with  engine  or  flight  parameters  other 
than  rpm  would  probably  not  be  productive.  This  conclusion  is  supported 
by  a  review  of  histograms  for  the  vibration  data  measured  in  jet  bombers 
and  fighters  during  various  different  flight  conditions.  The  review 
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indicates  that  variations  in  the  vibration  levels  at  shaft  frequencies 
and  their  harmonics  from  one  flight  condition  to  another  are  generally 
small  compared  to  the  variations  from  one  point  to  another  for  a  given 
flight  condition.  In  other  words,  the  variance  of  the  spatial  distribution 
within  a  zone  is  much  larger  than  the  variance  caused  by  changing  flight 
conditions.  This  point  is  illustrated  in  Table  22,  which  summarizes 
data  measured  in  Zone  02  (center  half  of  fuselage)  of  aircraft  Group  5 
(Century  series  fighters)  during  four  widely  varying  flight  conditions 
(takeoff,  climb,  normal  cruise,  cruise  with  afterburner).  The  data  are 
presented  for  the  frequency  interval  from  80  to  129  cps,  which  includes 
the  fundamental  of  shaft  rotation  (high  speed  shaft)  for  the  flight  con¬ 
ditions  in  question.  It  is  seen  in  Table  22  that  the  average  vibration 
levels  for  the  different  flight  conditions  vary  less  than  1,  5  to  1  while 
the  vibration  levels  for  any  one  flight  condition  vary  up  to  7 5* to  1. 

These  results  are  typical  of  the  data  for  all  zones  and  flight  conditions. 

If  the  past  AFFOLi  data  is  accepted  as  a  valid  representation  of 
the  periodic  vibration  environment  below  500  cps,  there  is  no  need  for 
further  data  collection  or  analysis  to  define  the  engine  shaft  induced 
vibration  in  jet  bombers  and  fighters.  The  desired  results  are  given 
directly  for  each  structural  zone  of  jet  bombers  by  Figure  11  (c),  and  for 
each  structural  zone  of  jet  fighters  by  Figure  11  (d).  The  results  are  pre¬ 
sented  as  mean  values  and  standard  deviations  for  displacement  amplitudes 
versus  frequency  for  each  zone.  A  conservative  limit  for  a  prediction  may 
be  obtained  from  these  results  by  assuming  a  specific  spatial  distribution 
function  for  the  data.  If  a  normal  distribution  for  peak  values  is  assumed, 
a  (1  -  a)  prediction  limit  for  the  ith  frequency  increment  is  given  by 


L. 

i; 


(l-a) 


z  B,  +  X. 
a  1  i 


(70) 
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Table  22.  Comparison  of  Vibration  Induced  by  Jet  Engine  Shaft  Rotation 
for  Various  Flight  Conditions 

Aircraft  Group  5  (Century  Series  Fighters) 

Structural  Zone  02  (Middle  Half  of  Fuselage) 

Frequency  Increments  13  and  14  (80  to  129  cps) 


Vibration 

Double 

Amplitude 

lO'^  in. 

Number  of  Data  Points  Measured  for 

Various  Flight  Conditions 

Condition  05 
(assisted 
takeoff) 

Condition  06 
(normal 
climb) 

Condition  07 
(normal 
cruise) 

Condition  48 
(cruise  with 
afterburner) 

8.  07-9. 

99 

1 

1 

6.  40-8. 

06 

1 

1 

5.  07-6. 

39 

1 

4.  02-5. 

06 

1 

2 

7 

3.  20-4. 

01 

9 

2 

14 

8 

2.  53-3. 

19 

28 

1 

13 

7 

2.  02-2. 

52 

24 

10 

30 

9 

1.  60-2. 

01 

39 

11 

39 

28 

1.  27-1. 

59 

34 

8 

37 

24 

1.  00-1. 

26 

27 

17 

80 

57 

0. 81-0. 

99 

30 

21 

81 

66 

0.  64-0. 

80 

31 

20 

159 

61 

0.  51-0. 

63 

50 

28 

179 

70 

0.  40-0. 

50 

63 

34 

243 

92 

0.  32-0. 

39 

71 

33 

286 

143 

0.  25-0. 

31 

73 

28 

223 

138 

0.  20-0. 

24 

35 

16 

84 

45 

0.  13-0. 

19 

4 

3 

16 

12 

Mean  Value 

1 0  in. 

0.  92 

0.  77 

0.  65 

0.  71 
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where 


z  =  100  a  percentage  point  of  standardized 
normal  distribution 

X.  =  mean  value  for  ith  frequency  interval  from 
Figure  1 1 

s.  =  standard  deviation  for  ith  frequency  interval  from 
Figure  11 

Further  discussions  of  the  statistical  aspects  of  selecting  prediction 

limits  are  presented  in  Section  8.  5. 

For  example,  assume  a  0.  975  (97.  5%)  prediction  limit  is  desired 

for  the  fundamental  shaft  vibration  in  Zone  02  (middle  half  of  fuselage) 

for  aircraft  Group  5  (Century  series  fighters).  The  fundamental  shaft 

frequency  for  the  high  speed  shaft  is  usually  in  the  frequency  increment 

from  80  to  99  cps.  From  Figure  11(d),  for  the  13th  frequency  increment, 
-  -3  -3 

X  =  0.  82  X  10  inch  and  8  =  0.  81x10  inch.  For  a  97.  5%  limit, 

z  =1.  96.  Hence,  L,  _  .  _o  =  0.  0024  inch  ps  1  g. 
a  13;  0.98  * 

8.  1.  5  Gunfire  Contribution 

The  vibration  induced  by  gunfire  is  basically  periodic  with  a  funda¬ 
mental  frequency  equal  to  the  firing  rate.  Hence,  the  contribution  of 
gunfire  is  easily  extracted  from  the  power  spectra  for  aircraft  vibration 
data  collected  during  the  operation  of  guns.  Unfortunately,  the  currently 
available  AFFDL  vibration  data  do  not  include  any  measurements  obtained 
during  gunfire.  Such  data  should  be  acquired  in  the  future.  The  contri¬ 
butions  of  gunfire  should  be  extracted  from  the  data  measured  at  all 
locations.  Mean  values  and  standard  deviations  should  then  be  computed 
for  the  data  in  each  structural  zone  at  the  firing  rate  frequency  and  its 
harmonics.  A  conservative  limit  for  future  predictions  maybe  obtained 
by  the  procedures  outlined  in  Section  8.  1. 4. 
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8.  1.  6  Suggested  Data  Reduction  Procedures 

The  prediction  models  suggested  in  this  section  for  jet  bomber 
and  fighter  vibration  environments  are  presented  in  terms  of  power 
spectra.  Theoretically,  the  data  used  to  develop  the  models  should 
then  be  reduced  in  terms  of  power  spectra.  In  practice,  however, 
this  is  not  always  necessary  or  even  desirable. 

First  consider  the  periodic  contributions  in  the  vibration  environ¬ 
ment  (jet  engine  shaft  rotation  and  gunfire).  The  power  spectrum  for 
a  periodic  vibration  is  theoretically  a  series  of  delta  functions.  Con¬ 
ventional  analog  and  digital  power  spectral  density  analysis  procedures, 
however,  will  display  each  periodic  component  (sinusoid)  as  a  peak  with 
finite  density  and  bandwidth.  These  peaks  must  then  be  converted  to  a 
mean  square  value  at  a  given  frequency  by  taking  the  product  of  the 
indicated  density  and  the  analysis  bandwidth.  The  data  may  be  further 
reduced  to  rms  or  peak  values,  as  desired.  From  the  data  reduction 
viewpoint,  the  only  requirement  is  that  the  analysis  bandwidth  be  suffi¬ 
ciently  narrow  to  distinguish  between  adjacent  periodic  components  and 
to  extract  the  periodic  components  from  the  background  random  contri¬ 
bution.  From  the  prediction  model  viewpoint,  the  data  pooling  operations 
(zoning)  will  produce  periodic  components  in  each  zone  at  numerous 
frequencies.  Hence,  the  predictions  must  be  made  in  terms  of  com¬ 
ponents  in  predetermined  frequency  increments  covering  the  frequency 
range  of  the  pooled  periodic  data.  The  selection  of  these  predetermined 
frequency  increments  need  not  be  influenced  in  any  way  by  the  bandwidth 
used  for  the  data  reduction.  The  increments  used  for  the  studies  in 
Sections  6  and  7  (see  Table  17)  are  acfc'eptable.  The  actual  calculations 
to  establish  a  conservative  prediction  level  should  be  performed  using 
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the  data  in  a  format  compatible  with  the  assumed  spatial  distribution 
function  for  the  data.  For  example,  if  a  lognormal  distribution  for 
mean  square  values  is  assumed,  the  data  should  be  evaluated  in  terms 
of  the  mean  square  values  for  the  periodic  components;  if  a  normal 
distribution  for  peak  values  is  assumed,  the  data  should  be  evaluated 
in  terms  of  the  peak  values  for  the  components,  etc. 

Now  consider  the  random  contributions  in  the  vibration  environ¬ 
ment  (jet  noise,  aerodynamic  noise,  wake  turbulence).  The  usual 
procedure  in  random  data  analysis  is  to  compute  power  spectra  using 
relatively  narrow  analysis  bandwidths  to  obtain  highly  resolved  spectral 
plots.  For  many  applications  (for  example,  transfer  function  analysis), 
highly  resolved  power  spectra  are  required.  Furthermore,  a  highly 
resolved  power  spectrum  provides  the  only  adequate  practical  method 
of  detecting  and  eliminating  periodic  components  in  the  data  spectrum. 
Hence,  it  appears  reasonable  that  the  data  should  be  analyzed  using 
relatively  narrow  analysis  bandwidths,  as  is  currently  done  by  AFFDL 
(see  Section  4).  This  does  not  mean,  however,  that  the  highly  resolved 
spectral  data  should  be  used  directly  in  the  regression  studies  to  estab¬ 
lish  a  vibration  prediction  model. 

It  is  wise  at  this  point  to  recall  the  ultimate  goal  of  the  regression 
analysis;  namely,  the  development  of  a  general  extrapolation  procedure 
for  predicting  aircraft  vibration  environments.  High  resolution  in  the 
final  predictions  is  not  believed  to  be  warranted  by  this  goal  since  con¬ 
siderable  frequency  blurring  will  already  have  been  introduced  by  the 
data  pooling  operations  required  to  develop  the  prediction  model.  It  is 
believed  that  predictions  with  relatively  coarse  resolution,  say  1/3 
octave  or  even  1/1  octave  bandwidths,  would  be  suitable.  This  predic¬ 
tion  resolution  fixes  the  resolution  bandwidth  for  the  power  spectra  data 
to  be  used  for  the  regression  analysis.  Hence,  for  this  application. 
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the  highly  resolved  power  spectra  should  be  converted  to  average 
power  spectral  densities  in  1/3  or  1/1  octave  bands.  This  may  be 
done  by  estimating  the  area  under  the  power  spectrum  between  the 
frequency  limits  of  the  1/3  or  1/1  octave  band  in  question,  excluding 
the  area  under  peaks  due  to  periodic  components.  This  procedure  is 
illustrated  for  the  1/1  octave  case  in  Figure  18. 


Figure  18.  Illustration  of  Random  Data  Analysis  Procedures 

One  problem  is  introduced  by  the  above  suggested  data  reduction 
procedures.  Since  the  data  and,  hence,  the  resulting  predictions  are  in 
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terms  of  averajge  power  spectra  in  1/3  or  1/1  octave  bands,  they  may 
well  underpredict  peaks  in  the  power  spectra  for  the  actual  random 
vibration  environment.  This  possibility  may  be  accounted  for  by  adding 
a  general  correction  term  to  the  resulting  predictions.  The  correction 
term  would  be  determined  empirically  by  noting  the  ratio  of  the  highest 
peak  to  average  density  for  each  1/3  or  1/1  octave  of  each  power  spec¬ 
trum  during  the  data  evaluation  illustrated  in  Figure  18.  Of  course, 
peaks  due  to  periodic  components  would  be  ignored.  The  mean  value 
and  standard  deviation  for  these  observed  peak  to  average  ratios  would 
then  be  computed.  The  correction  factor  would  be  determined  by 
selecting  the  ratio  corresponding  to  some  small  percentage  point  on  a 
normal  distribution  with  the  calculated  mean  and  standard  deviation. 

The  determination  of  a  correction  term  of  this  type  for  missile  vibration 
predictions  is  reported  in  [6]  .  The  results  suggest  a  correction  factor 
of  3  dB  for  1/3  octave  predictions  and  5  dB  for  1/1  octave  predictions. 

8.  Z  MODEL  FORMULATION  FOR  HELICOPTERS 

The  available  AFFDL  data  cover  both  reciprocating  and  jet  engine 
powered  helicopters  (one  of  each).  Since  the  design  of  future  reciprocat¬ 
ing  engine  powered  helicopters  is  unlikely,  jet  powered  helicopters  are 
of  primary  interest.  However,  the  principal  sources  of  vibration  for  the 
two  are  similar.  Hence,  data  from  both  types  may  be  used  as  a  basis 
for  predicting  the  vibration  environment  of  the  jet  powered  type. 

From  Table  6  in  Section  3,  there  are  seven  sources  of  vibration 
for  helicopters  which  are  considered  sufficiently  significant  to  include 
in  a  prediction  model.  In  review,  these  sources  are  as  follows. 


1.  D-1  —  rotor  blade  passage 

2.  D-2  —  rotor  rotation 

3.  D-3  —  jet  engine  shaft  rotation  (Jet  engine  powered 

only) 

4.  D-4  —  reciprocating  engine  exhaust  (reciprocating 

engine  powered  only) 

5.  D-6  —  engine  accessory  equipment 

6.  D-8  —  gear  box  noise 

7.  D-10  —  gunfire 

In  attempting  to  develop  a  practical  prediction  model  for  the  above 
sources,  two  problems  quickly  evolve.  First,  from  the  discussions  in 
Section  3.  3,  many  of  these  sources  are  poorly  defined  in  terms  of 
tractable  analytical  models.  Second,  all  of  the  sources,  except  gunfire, 
occur  simultaneously.  That  is,  it  is  difficult  to  isolate  the  contributions 
of  individual  sources  as  is  done  for  jet  bombers  and  fighters  in  Sec¬ 
tion  8.  1.  For  these  reasons,  an  attempt  to  develop  a  prediction  model 
for  helicopter  vibration  based  upon  individual  contributing  sources  is 
not  considered  worthwhile.  Instead,  the  problem  should  be  approached 
by  evaluating  the  vibration  environment  from  all  sources  measured  for 
various  pertinent  flight  conditions.  Currently  availi  DL  data 

should  be  acceptable  for  this  task,  at  least  in  the  freqii.  irv  range  below 
500  cps,  since  the  principal  sources  of  vibration  are  peiiodic. 

Histograms  for  the  vibration  levels  measured  by  AFFDL  in  the 
three  basic  structural  zones  of  helicopters  during  takeoff,  climb,  cruise, 
and  hover  are  presented  in  Section  5  of  Appendix  A.  An  inspection  of 
these  data  reveals  that  the  vibration  levels  for  the  four  different  flight 
conditions  are  surprisingly  similar.  This  fact  is  illustrated  in  Table  23, 
which  summarizes  the  mean  values  in  various  frequency  increments  for 
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Table  23.  Comparison  of  Average  Helicopter  Vibration  Levels  for  Various  Flight  Conditions 


400-499  51  0.087  30  1  0.078  160  0.064  I  82  0.061 


Zone  02  (middle  half  of  fuselage).  Note  that  results  are  presented  only 
for  those  frequency  increments  where  a  significant  number  of  data 
values  occur.  It  is  seen  from  Table  23  that  the  mean  values  for  the 
vibration  levels  measured  during  the  four  flight  conditions  differ  very 
little  in  most  frequency  increments.  Considering  the  sample  size  for 
the  mean  value  calculations,  many  of  the  indicated  differences  are 
not  even  statistically  significant.  These  results  are  typical  of  those 
for  the  other  zones  and  flight  conditions.  Hence,  if  the  past  AFFDL 
data  are  accepted  as  a  valid  representation  of  the  periodic  vibration 
environment  below  500  cps,  it  must  be  concluded  that  the  variation  of 
helicopter  vibration  levels  for  different  flight  conditions  is  negligible 
compared  to  the  variation  with  structural  location.  This  means  that 
the  vibration  environment  for  helicopters  (similar  to  those  surveyed 
by  AFFDL)  can  be  predicted  by  a  single  spectrum  for  each  structural 
zone,  independent  of  flight  condition.  The  desired  statistical  information 
is  given  directly  by  Figure  11(e).  By  fitting  an  appropriate  spatial  distri¬ 
bution  function  to  the  mean  values  and  standard  deviations  in  Figure  11(e), 
a  prediction  for  helicopter  vibration  amplitudes  versus  frequency  may 
be  determined  at  any  desired  percentage  point.  Such  prediction  curves 
for  the  three  basic  structural  zones  of  helicopters,  based  upon  the  2,  5% 
point  of  a  normal  distribution,  are  presented  in  Figure  19. 

Note  that  the  predictions  in  Figure  19  do  not  include  the  contribution 
of  gunfire.  No  measurements  obtained  during  gunfire  are  included  in  the 
currently  available  AFFDL  vibration  data.  Such  data  should  be  obtained 
in  the  future,  and  analyzed  to  generate  prediction  curves  similar  to  those 
presented  in  Figure  19.  These  new  curves  would  apply  to  those  flight 
phases  when  guns  are  used.  Further  note  that  the  predictions  in  Figure  19 
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probably  do  not  reflect  the  vibration  induced  by  "blade  slap"  during  high 
forward  speed  flight  or  sharp  maneuvers.  These  conditions  should  be 
analyzed  and  predicted  separately,  as  discussed  in  Section  3.  1. 

8.  3  REGRESSION  ANALYSIS  PROCEDURES 

From  Section  8.  1,  there  are  two  contributing  sources  of  vibration 
for  jet  bombers  and  fighters  which  must  be  evaluated  by  a  regression 
analysis  of  appropriate  flight  vibration  data.  These  sources  are  jet 
noise  and  aerodynamic  boundary  layer  noise.  The  suggested  regression 
models  for  these  two  sources  are  given  by  Eqs.  (64)  and  (67),  re¬ 
spectively.  In  both  cases,  the  model  is  of  the  same  general  form  and 
involves  only  one  independent  variable.  For  simplicity,  the  model  will 
be  written  in  the  form 


y  =  Ax 


(71) 


where  y  =  G(f)  and  x  =  X.  The  problem  now  is  to  establish  a  (1  -  a) 
upper  prediction  limit  for  y  given  x,  based  upon  a  collection  of  mea¬ 
sured  values  for  x  and  y . 

Let  y.  denote  the  ith  measurement  and  x.  denote  the  independent 
variable  associated  with  the  ith  measurement,  where  i  =  1,  2,  3,  ...,  in. 
The  required  calculations  are  as  follows  [34],  [35]. 


(72a) 
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(72g) 


t 


n-  1 ;  a 


the  a  percentage  point  of  a  student's  t  distribu¬ 
tion  with  n  -  1  degree-of-freedom 


s  is  defined  in  (72e) 
yx 


The  regression  coefficient  A  is  given  by  Eq.  (72a).  The  other  calcula 
tions  are  used  to  evaluate  the  model  and  select  prediction  limits,  as 
follows. 


8.  3.  1  Evaluation  of  Model  Linearity 

It  is  initially  assumed  that  the  relationship  between  x  and  y  is 
linear.  The  validity  of  the  assumption  should  be  checked.  This  may 
be  done  using  the  value  of  $  given  by  Eq.  (72b).  If  the  model  is  nonlinear, 
the  value  of  $  will  be  distributed  like  the  variable,  F  with  degrees-of- 
freedom  given  by  df^  =  1  and  df^  =  n  -  1[35].  Hence,  a  hypothesis  of  non¬ 
linearity  may  be  tested  by  comparing  the  value  of  $  with  the  tabulated  * 

value  for  F  at  any  desired  level  of  significance  o.  If  $  <  F 

1,  n- i; a  1,  n- 1;  o 

the  model  is  nonlinear.  If  §  >  F,  ,  ,  there  is  no  reason  to  question 

1 ,  n- 1;  O'  ’ 

that  the  model  is  linear. 


8.  3.  2  Evaluation  of  Model  Efficiency 

The  efficiency  of  the  model  may  be  evaluated  in  terms  of  the 

2 

correlation  coefficient  y  given  by  Eq.  (72c).  The  quantity  (1  -  y  ) 
is  a  measure  of  the  power  contributed  to  y  by  variables  other  than  x . 

If  the  model  provided  a  perfect  description  for  y  given  x  ,  y  would 
equal  unity.  On  the  other  hand,  if  there  were  no  relationship  between 
y  and  x,  y  would  equal  zero.  It  is  clearly  desirable  to  establish  if 
y  is  significantly  different  from  zero.  This  is  accomplished  using  the 
Fisher  "Z"  transformation  given  by  Eq.  (72d).  If  y  =  0  ,  the  value  of 
Z  will  be  normally  distributed  with  a  mean  of  zero  and  a  variance  of 
unity.  Hence,  a  hypothesis  of  significant  correlation  may  be  tested 
by  comparing  the  value  of  Z  with  the  tabulated  value  of  the  standardized 
normal  distribution  at  any  desired  percentage  point  a.  If  Z  <  z  ,  the 

Or 

correlation  coefficient  is  not  significantly  different  from  zero.  If 
Z  >  z^ ,  the  correlation  coefficient  is  significantly  different  from  zero. 
Note  that  the  foregoing  interpretations  of  y  are  valid  only  if  the  assump¬ 
tion  that  Aq  =  0  is  valid. 

The  efficiency  of  the  model  can  also  be  evaluated  in  terms  of  the 
(1  -  a)  confidence  interval  for  A,  as  given  by  Eq.  (72e);  if  the  lower 
confidence  limit  is  less  than  zero,  the  interpretation  is  the  same  as  for 
the  case  where  the  correlation  coefficient  is  not  significantly  different 
from  zero. 

8.  3.  3  Selection  of  Prediction  Limits 

The  final  prediction  limit  is  selected  at  any  desired  percentage 

point  using  Eqs.  (72f)  and  (72g).  The  term  Ly  provides  a  (1  -  or) 

y  K|  ct 

prediction  limit  for  each  value  of  y.  for  a  given  value,  x  ,  of  the  inde- 

X  iC 

pendent  variable.  Further  discussions  of  the  statistical  aspects  of 
selecting  prediction  limits  are  presented  in  Section  8.  5. 
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8.  4  ILLUSTRATION  OF  REGRESSION 
ANALYSIS  PROCEDURES 


The  procedures  outlined  in  Section  8.  3  are  now  illustrated  using 
data  collected  by  AFFDL  during  the  flight  vibration  survey  of  aircraft 
Code  No.  63  (RF-4C)  [26] .  Since  data  from  only  aircraft  is  available, 
the  illustration  is  limited  to  a  study  of  the  suggested  model  for  aero¬ 
dynamic  noise  induced  vibration,  as  given  by  Eq.  (67).  A  study  of  the 
suggested  model  for  jet  noise  induced  vibration  is  not  practical  because 
a  single  aircraft  provides  data  for  only  one  value  of  the  independent 

variable  X  ,  in  Eq.  (64). 
r  1 

Numerous  sample  records  of  the  vibration  at  various  locations  in 
Zone  02  of  the  RF-4C  were  provided  by  AFFDL  for  analysis.  The  sample 
records  were  reviewed  based  upon  the  structural  locations  and  flight  con¬ 
ditions  for  the  measurements.  Only  those  records  which  represented 
measurements  on  basic  structure  during  normal  cruise  flight  were  re¬ 
tained.  Since  the  independent  variable  in  Eq.  (67)  is  proportional  to 
2  4 

p  V  ,  every  effort  was  made  to  select  records  covering  a  wide  range 
of  airspeeds  and  altitudes.  However,  sample  records  for  flight  at  Mach 
numbers  around  0.  8  were  omitted  because  the  RF-4C  aircraft  displays 
unique  and  unrepresentative  vibration  characteristics  around  this  Mach 
number  [31] .  The  sample  records  which  passed  the  review  were  then 
carefully  edited  for  data  quality.  Only  those  sample  records  with  an 
acceptable  signal  to  noise  ratio  over  most  of  the  frequency  range  of 
interest  (10  to  3000  cps)  were  retained.  The  review  and  editing  pro¬ 
cedures  produced  a  total  of  20  acceptable  sample  records  covering  uhree 
transducer  locations  and  nine  combinations  of  airspeed  and  altitude. 

This  information  along  with  estimates  for  the  average  surface  weight 
density  at  each  measurement  location  are  detailed  in  Table  24. 

The  sample  records  for  the  measurements  summarized  in  Table  24 
were  reduced  to  average  power  spectra  in  1/1  octave  bands.  The  results 
are  presented  in  Table  25.  Note  that  data  at  frequencies  below  22  cps  are 
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surface  weight  density  of  the  noted  structural  location 
dynamic  pressure  for  the  noted  flight  condition 


Table  25.  Average  Power  Spectra  for  RF-4C  Vibration  Data 


Average  Power  Spectral  Densities 

in 

Xr2 

g^/cps  for  1/1  Octave  Bands  in  cps 

Code 

No. 

22- 

44- 

88- 

177- 

354-  ' 

707- 

1414- 

X  10® 

44 

88 

177 

354 

707 

1414 

2828 

1 

.0118 

.0016 

.0010 

.00071 

.00036 

. 00023 

.000036 

.  000032 

2 

1 

.  0596 

.0016 

.  0009 

.0011 

.0004 

.00018 

. 000045 

.000023 

3 

.  189 

.  0026 

.  018 

.032 

.0018 

.  0008 

.00007 

.  000028 

4 

.  449 

.  001 

.  0014 

.0008 

.013 

.013 

.  013 

.  0057 

5 

.  934 

.0064 

.071 

.  50 

.04 

.013 

.0016 

.00028 

6 

1.703 

.0018 

. 00079 

.  0008 

.04 

.040 

.050 

.020 

7 

2.  117 

.  0036 

. 00079 

.0005 

.025 

.  032 

.  057 

.  025 

8 

.  0118 

-- 

.00063 

.  0002 

.00016 

.  00025 

. 00008 

. 000023 

9 

.  0596 

.00073 

.0004 

.0003 

.0002 

.0004 

.00028 

.00025 

10 

.  189 

.0008 

.  0004 

.00018 

.0010 

.0018 

.0018 

.0004 

11 

.  934 

.  0013 

. 00071 

.00089 

.0113 

.018 

.  022 

.  0040 

12 

1.703 

.0014 

.0007 

.0005 

.013 

.023 

.036 

.0063 

13 

.  037 

.  0016 

.00089 

.00028 

.00023 

.  00007 

.000028 

. 000018 

14 

.  189 

.0032 

.0013 

.0004 

.00057 

.0029 

.  0025 

. 00028 

15 

.  596 

.  0018 

.0010 

.00045 

.004 

.011 

.014 

.  0008 

16 

1.417 

-- 

.  0016 

.0014 

.014 

.  013 

.029 

.  0028 

17 

2.  95 

.  115 

.  079 

.090 

.  102 

.  10 

.  14 

.023 

18 

5.  38 

.  58 

.29 

.  18 

.  18 

.  13 

.  16 

.051 

19 

.  538 

.0020 

. 00089 

.  0004 

.0036 

.  Oil 

.016 

.00036 

20 

.  43 

.0020 

.  0010 

.00036 

.0028 

.008 

.011 

.00025 

1 
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not  included  because  most  records  did  not  provide  an  adequate  signal 
to  noise  ratio  below  this  frequency.  Also  note  that  the  contributions  of 
periodic  components  were  not  removed  from  the  data,  as  suggested  in 
Soc  lion  8.  1.  2.  Periodic  components  due  to  jet  engine  shaft  rotation  were 
apparent  in  the  88  to  177  cps  octave.  They  were  not  removed  to  illustrate 
a  point,  which  is  discussed  later. 

The  d  l  a  in  Table  25  are  now  analyzed  using  the  procedures  out- 

I 

lined  in  Section  8.  3.  The  detailed  results  of  the  regression  analysis  are 
presented  in  Section  6  of  Appendix  A  and  summarized  in  Table  26.  It 
is  seen  from  Table  26  that  the  model  is  linear  and  provides  significant 
correlation  at  the  1%  level  of  significance  for  all  octave  bands  except 
the  88  to  177  cps  band.  This  is  the  frequency  band  which  included 
periodic  components  due  to  jet  engine  shaft  rotation.  Hence,  the  results 
in  this  band  clearly  illustrate  the  necessity  for  removing  periodic  contri¬ 
butions  from  the  data.  They  further  illustrate  the  ability  of  the  analysis 
procedures  to  detect  the  presence  of  contributing  vibration  sources  other 
than  aerodynamic  boundaiy  layer  noise. 

The  results  in  Table  26  confirm  that  the  prediction  model  for  aero¬ 
dynamic  noise  induced  vibration,  as  suggested  in  Eq.  (67),  is  acceptably 
linear  and  efficient.  The  indicated  correlation  coefficients  for  all  octaves 
(excluding  the  88  to  177  cps  octave)  are  from  0.  68  to  0.  90.  This  is  re¬ 
latively  high  considering  the  simplicity  of  the  model.  The  efficiency  of  the 
model  is  further  confirmed  by  the  fact  that  the  lower  confidence  limits 
for  the  coefficient  A  are  all  greater  than  zero.  Of  course,  when  these 
analysis  procedures  are  ultimately  applied  to  data  from  many  different 
aircraft,  the  results  may  not  be  as  good  as  achieved  here  using  data  from 
only  one  aircraft.  Nevertheless,  it  is  believed  that  they  will  be  meaning¬ 
ful  and  provide  an  acceptable  basis  for  establishing  conservative  vibration 
prediction  limits. 
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Table  26.  Summary  of  Regression  Analysis  Results 
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The  97.  5%  upper  prediction  limits  calculated  for  this  illustration 
are  presented  in  Figure  20.  Curves  of  this  type  represent  the  ultimate 
goal  of  the  analysis.  Note  that  the  predictions  are  in  terms  of  average 
power  spectral  densities  in  1  / 1  octave  bands.  If  it  is  desired  to  have 
predictions  which  will  envelop  anticipated  narrow  band  peaks  in  the  vibra¬ 
tion  power  spectra,  a  correction  must  be  applied  to  the  curves  in  Fig¬ 
ure  20,  as  discussed  in  Section  8.  1.  f'.  !<■  should  be  emphasized  that  the 

prediction  curves  in  Figure  20  are  presented  for  illustrative  purposes, 
only.  Because  of  the  limited  amount  of  data  used  to  arrive  at  the  curves, 
they  do  not  represent  a  valid  tool  for  predicting  the  vibration  environ¬ 
ment  of  future  aircraft. 

8.  5  STATISTICAL  CONSIDERATIONS  IN 

THE  SELECTION  OF  PREDICTION  LIMITS 

As  used  in  this  report,  a  (1  -  or)  prediction  limit  is  that  level  calcu¬ 
lated  from  past  vibration  measurements  v/hich  hopefully  will  exceed 
100(1  -  a)%  of  all  future  vibration  levels.  In  rigorous  terms,  such  a 
prediction  limit  should  be  based  upon  a  statistical  tolerance  limit  [36] . 

A  tolerance  limit  is  defined  as  that  level  calculated  from  past  observations 
which  will  exceed  at  least  100  (1  -  «)%  of  all  future  observations  with  a 
probability  of  P.  For  n  observations  of  a  normally  distributed  random 
variable  y ,  the  tolerance  limit  for  future  observations  is  given  by 


'P;a 


=  X  + 


s  K  _ 

X  n;  P;  a 


(73) 


where  x  and  s  are  the  sample  mean  and  standard  deviation. 


respectively. 


AVERAGE  OCTAVE  BAND  POWER  SPECTRAL  DENSITY  —  g/CPS 


and  K  is  the  tolerance  factor.  Tolerance  factors  for  various  values 
of  n  ,  P,  and  a  are  widely  tabulated  in  the  statistical  literature. 

The  prediction  limits  arrived  at  by  the  procedures  suggested  in 
this  report  are  not,  rigorously  speaking,  statistical  tolerance  limits. 

For  the  case  of  the  regression  analysis  procedures  described  in  Sec¬ 
tion  8.  3,  the  prediction  limit  given  by  Eq.  (72g)  is  similar  to  a  tolerance 
limit  [34].  For  the  other  cases,  however,  a  (1  -  n)  prediction  limit  is 
arrived  at  by  fitting  an  assumed  distribution  function  to  the  data,  and 
determining  the  a  percentage  point  of  the  fitted  distribution.  The  re¬ 
sulting  prediction  limit  would  be  equivalent  to  a  tolerance  limit  only  if 
the  fitted  distribution  were  in  fact  the  actual  distribution  function  for  the 
data.  This,  of  course,  is  not  the  case.  Nevertheless,  the  use  of  the 
percentage  point  procedure  for  selecting  prediction  limits  is  believed 
justified  for  the  following  two  reasons. 

1.  Tolerance  factors  are  well  tabulated  only  for  the  case 
of  normally  distributed  data.  The  data  of  interest  here 
are  not  normally  distributed,  as  discussed  in  Section  7. 

2.  The  potential  error  introduced  by  using  percentage 
points  of  fitted  distributions,  rather  than  tolerance  factors, 
to  establish  prediction  limits  is  negligible  compared  to 
the  potential  errors  associated  with  the  other  steps  needed 
to  arrive  at  vibration  predictions. 


9.  SELECTION  OF  VIBRATION  TEST  LEVELS 
AND  DURATIONS 


Having  predicted  an  aircraft  vibration  environment,  it  remains  to 
convert  the  predictions  into  practical  specifications  for  aircraft  com¬ 
ponent  vibration  tests.  Procedures  for  selecting  appropriate  test  levels 
and  durations  based  upon  environmental  predictions  are  investigated  in 
this  section. 

Historically,  there  have  been  two  basic  approaches  to  the  design 
of  vibration  test  specifications,  the  environmental  simulation  approach 
and  the  damage  simulation  approach.  The  environmental  simulation 
approach  involves  the  design  of  a  vibration  test  specification  which  will 

fe- 

conservatively  simulate  the  environmental  vibration  levels  and  exposure 
times  to  be  expected  in  actual  service.  This  approach  is  clearly  most 
suitable  for  missiles  and  spacecraft  applications  where  the  total  dura¬ 
tion  of  the  service  vibration  environment  is  relatively  shoit  and,  hence, 
can  be  reasonably  simulated  by  a  laboratory  test.  The  damage  simu¬ 
lation  approach  involves  the  design  of  a  vibration  test  specification  which 
conservatively  simulates  the  damaging  potential  of  the  environment, 
as  opposed  to  the  actual  environmental  levels.  The  use  of  this  approach 
permits  long  environmental  exposure  times  in  actual  service  to  be  re¬ 
placed  by  short  duration  vibration  tests.  This  particular  property  makes 
the  damage  simulation  approach  highly  desirable  for  aircraft  component 
vibration  test  specifications.  The  use  of  the  damage  simulation  approach, 
however,  does  pose  one  serious  problem.  Specifically,  it  is  necessary 
to  establish  a  model  for  the  mode  and  mechanism  of  failure  to  be  antici¬ 
pated  in  the  component  to  be  tested.  Only  on  the  basis  of  some  assumed 
mode  of  failure  can  proper  tradeoffs  between  vibration  level  and  vibration 
duration  be  established. 
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There  is  an  additional  problem  associated  with  the  selection  of 
vibration  test  durations  which  applies  to  test  specifications  based  upon 
either  the  environmental  simulation  or  damage  simulation  approach. 

This  problem  evolves  from  the  fact  that  laboratory  vibration  tests  are 
usually  performed  using  a  stationary  vibration  input,  while  the  vibration 
environment  of  flight  vehicles  in  actual  service  is  generally  nonsta¬ 
tionary.  Hence,  the  following  question  arises:  How  long  should  the 
duration  for  the  stationary  vibration  test  be  to  properly  simulate  a  non¬ 
stationary  vibration  environment?  For  the  case  of  missile  and  space¬ 
craft  component  testing  where  an  environmental  simulation  approach  is 
used,  the  general  procedure  is  to  establish  the  test  levels  based  on  the 
maximum  vibration  levels  which  are  anticipated  during  service,  and  to 
base  the  test  duration  upon  the  total  exposure  time  in  service.  This 
approach  often  produces  an  overly  conservative  test,  but  avoids  the 
problem  of  establishing  equivalences  between  nonstationary  and  sta¬ 
tionary  vibration.  For  the  case  of  test  specifications  based  upon  a 
damage  simulation  approach,  the  nonstationary- stationary  equivalence 
problem  can  be  dealt  with  more  directly.  In  fact,  it  can  be  handled  by 
exactly  the  same  procedures  used  to  make  tradeoffs  between  test  levels 
and  test  durations  to  compress  long  service  exposure  times  into  short 
test  durations.  This  fact  constitutes  a  significant  advantage  of  the 
damage  simulation  approach  over  the  environmental  simulation  approach. 

9.  1  FAILURE  CRITERIA 

The  possible  modes  of  failure  for  aircraft  components  are  numer¬ 
ous  and  generally  quite  complex  [37],  [38].  For  the  problem  at  hand, 
however,  it  is  necessary  to  hypothesize  that  failures  will  be  due  to  some 


specific  mechanism  which  can  be  described  by  a  tractable  mathematical 
model.  Two  different  types  of  failure  criteria  will  be  assumed  and  used 
independently.  The  first  is  a  fatigue  damag-'  criterion  and  the  second 
is  a  peak  value  criterion. 

9.  1.  1  Fatigue  Damage  Criterion 

The  exact  mechanics  of  structural  fatigue  are  not  well  defined. 
However,  a  number  of  models  for  describing  fatigue  damage  have  been 
proposed  over  the  years  which  produce  reasonable  fatigue  life  predic¬ 
tions.  One  of  the  earliest  and  simplest  of  these  models  is  Miner's 
hypothesis  [39],  which  maybe  stated  as  follows.  If  stress  cycles 
with  a  peak  stress  level  of  will  cause  a  fatigue  failure  of  a  given 
material,  then  nj^  <  stress  cycles  with  that  same  peak  stress  level 
will  expend  fractional  portion  of  the  material's  useful  fatigue 

life.  Hence,  if  streams  cycles  are  applied  at  various  different  peak  stress 
levels  (S. ,  S_,  . .  .  ,  S.  . . . ) ,  a  measure  of  fatigue  damage  is  given  by 

1m  1 

(74) 

i  i 


where  failure  is  assumed  to  occur  when  D  =  1  .  The  relationship  between 
the  peak  stress  level,  S^,  and  the  number  of  cycles  to  failure,  ,  is 
given  by  an  S  -  N  curve  for  the  material. 

S  -  N  curves  for  common  aircraft  materials  vary  widely  depending 
upon  the  alloy,  heat  treatment,  mean  stress,  stress  raisers,  and  other 
factors.  For  many  materials,  however,  the  S-N  curve  can  be  crudely 


approximated  by  two  straight  line  segments  when  presented  on  a  log-log 
plot,  as  illustrated  in  Figure  21.  The  stress  level  which  corresponds 
to  the  level  of  the  horizontal  segment  in  Figure  21  is  called  the  en¬ 
durance  limit.  The  S  -  N  curves  for  aluminum  alloys  do  not  generally 
break  as  abruptly  as  indicated  in  Figure  21;  that  is,  the  endurance  limit 
is  not  so  ideally  defined.  Nevertheless,  there  is  a  stress  limit  where 
stress  cycles  with  peak  values  below  that  limit  are  unlikely  to  produce  a 
fatigue  failure,  even  as  the  number  of  stress  cycles  becomes  very  large. 

For  the  moment,  let  the  break  in  the  S  -  N  curve  at  the  endurance 
limit  be  ignored.  The  S  -  N  curve  for  stress  levels  above  the  endurance 
limit  can  be  approximated  analytically  by 


(75) 


where  A  and  b  are  material  constants.  Substituting  Eq.  (75)  into 
Eq.  (74)  yields 


(76a) 


which  may  be  applied  directly  to  approximate  the  fatigue  damage  of  a 
material  due  to  sinusoidal  stress  loads  with  various  different  peak  levels. 
For  the  case  of  random  stress  data,  Eq.  (76a)  reduces  to 


where  p(S)  is  the  probability  density  function  for  peak  values  of  stress, 
and  n  is  some  equivalent  total  number  of  cycles. 


9.  1.  2  Peak  Value  Criterion 

Consider  now  those  failures  which  are  due  to  a  single  extreme 
value  of  strain,  stress,  velocity,  or  acceleration.  Included  would  be 
failures  due  to  a  stress  level  beyond  the  ultimate  strength  of  the  material, 
a  displacement  beyond  a  physical  space  limit  (a  collision  or  misalignment), 
or  an  acceleration  beyond  a  functioiial  performance  limit.  For  the  case 
of  periodic  vibration  environments,  such  peaks  are  easy  to  define  since 
they  occur  repeatedly  with  the  same  magnitude  for  each  cycle  of  vibration. 
For  random  vibration  environments,  however,  the  problem  is  more  com¬ 
plicated  since  random  data  generally  have  no  clearly  defined  peak  value. 
Hence,  the  possibility  that  a  random  vibration  will  or  will  not  exceed  any 
given  level  must  be  specified  by  a  probability  statement.  The  probability 
that  a  given  level  will  be  exceeded  is  a  function  of  the  power  spectrum 
and  probability  density  function  for  the  vibration,  and  the  total  observa¬ 
tion  time.  An  exact  closed  form  solution  for  the  probability  of  exceeding 
a  given  level  is  not  known  to  exist.  However,  for  the  case  of  vibration 
with  a  Gaussian  probability  density  function,  it  can  be  shown  [40]  that 
an  upper  bound  on  the  probability  of  exceeding  a  level  K  =  L/o-  >  3  is 
giv'en  by 


„  nx  -K^/2 
P  <  —  e 

—  IT 


(77) 
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where  K  is  the  ratio  of  the  peak  to  rms  vibration,  T  is  the  total 
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observation  time,  and  is  tt  times  the  expected  number  of  zero  cross¬ 
ings  per  unit  time.  In  terms  of  the  power  spectrum  G(f)  for  the  vibra¬ 
tion 


f 


Ztt 


^  G(f)  df 


G(f)  df 


1/2 


(78) 


There  are  other  approximations  for  the  probability  of  exceeding  a  given 
level  which  have  been  proposed  and  studied  by  various  people.  Included 
are  those  suggested  by  MacNeal,  Barnoski,  and  Bailie  [41]  ,  and  Mark 
[42]  ,  which  apply  to  the  response  of  simple  oscillators  to  random  ex¬ 
citation.  The  result  in  Eq.  (77),  however,  is  generally  simpler  in  form 
than  other  suggestions  and  is  considered  to  be  satisfactory  for  the  appli¬ 
cation  of  interest  herein. 


i 


i 


9.  2  APPLICATION  TO  COMPONENT  MODEL 

The  failure  criteria  discussed  in  Section  9.  1  must  now  be  applied 
to  a  specific  mechanical  model  which  approximates  the  dynamics  of  air¬ 
craft  components.  The  model  to  be  used  will  be  the  simple  mechanical 
oscillator  illustrated  in  Figure  22,  where  m  is  the  mass,  k  is  the  spring 
constant,  c  is  the  viscous  damping  coefficient,  and  a(t)  is  the  accelera¬ 
tion  of  the  foundation. 


199 


The  simple  model  in  Figure  22  is  clearly  a  gross  oversimplification 
for  the  complex  dynamic  characteristics  of  most  aircraft  components. 
However,  it  must  be  remembered  that  the  model  is  required  only  as  a 
basis  for  comparing  the  damaging  effects  of  two  different  environments 
on  the  same  component.  Furthermore,  if  it  is  assumed  that  sinusoidal 
tests  will  be  used  to  simulate  periodic  vibration  environments  and  random 
tests  will  be  used  to  simulate  random  vibration  environments,  the  en¬ 
vironmental  differences  of  interest  will  be  principally  dif'  irences  in  level 
and  duration,  and  not  differences  in  the  basic  characteristics  of  the 
vibration.  For  these  reasons,  it  is  believed  the  model  in  Figure  22  will 
provide  an  acceptable  basis  for  the  desired  comparisons. 

For  either  the  fatigue  damage  or  peak  value  criterion,  the  re¬ 
sponse  function  of  interest  for  the  model  in  Figure  22  is  the  displace¬ 
ment  of  the  mass  relative  to  the  foundation  (the  strain  in  the  spring). 

It  can  be  readily  shown  [43]  that  the  magnitude  of  the  frequency  response 
function  for  this  system  which  relates  an  acceleration  input  to  a  strain 
response  is 
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where  T  is  the  total  exposure  time.  Hence,  the  fatigue  damage  to  the 
model  is  approximated  by  Eq.  (76a)  as 


For  the  worst  case  where  resonant  vibration  occurs,  Eq.  (82)  becomes 


where  Q  =  1  /2  ^  . 

Now  consider  the  case  where  the  vibration  environment  is  a 
Gaussian  random  acceleration  with  a  uniform  power  spectrum  of  G ,  at 
least  over  those  frequencies  in  the  region  of  resonance.  Assuming  the 
system  Q  is  high  (greater  than  10),  it  can  be  sho  /n  [44]  that  the 
fatigue  damage  is  approximated  by 


where  r(  )  is  the  Gamma  function. 


(84) 


202 


9.  2.  2  Application  of  Peak  Value  Criterion 


For  the  case  where  the  vibration  environment  is  a  sinusoidal 
acceleration,  the  peak  value  criterion  is  of  no  use  since  any  given  level 
either  will  or  will  not  be  exceeded  on  the  first  cycle  of  vibration  (after 
transients  have  subsided).  In  other  words,  the  probability  of  a  failure 
is  not  dependent  upon  the  observation  time  beyond  one  vibration  periodic. 

For  the  case  where  the  vibration  environment  is  a  Gaussian  random 
acceleration,  tt  times  the  expected  number  of  zero  crossings  for  the  re¬ 
sponse,  as  defined  in  Eq.  (78),  becomes  =  2'n-f  .  Hence,  from  Eq.  (77), 

n 

the  probability  of  exceeding  any  given  level  K  =  L/cr  is  bounded  by 


P  < 


2f  T  e 
n 


-K^/2 


(85) 


9.  3  RELATIONSHIPS  BETWEEN  VIBRATION 
LEVEL  AND  EXPOSURE  TIME 

Assume  it  is  desired  to  substitute  a  short  duration  vibration  test 
for  a  long  duration  service  vibration  environment.  The  following  question 
arises:  How  much  should  the  test  level  be  increased  to  produce  the  same 
damaging  potential  as  the  service  environment?  This  question  will  now 
be  answered  using  the  results  in  Section  9.  2. 

9.  3.  1  Fatigue  Damage  Criterion 

For  the  case  of  sinusoidal  vibration  with  a  fatigue  damage  criterion, 
the  relationship  between  the  test  and  service  environmental  levels  and 
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durations  is  obtained  by  equating  damage  for  the  two  cases  using  Eq.  (82). 
Since  the  same  component  is  involved  in  both  cases,  it  follows  that 


_  b  _  b 
T  a  =  T  a 
t  t  s  s 


(86) 


For  random  vibration  with  a  fatigue  damage  criterion,  the  relationship 
is  given  by  Eq.  (84)  as 


^b/2  ^  ^b/2 

T  G  =  T  G 
t  t  s  8 


(87) 


Assuming  that  the  power  spectrum  for  the  test  and  the  8er\dce  environ¬ 
ments  are  similar  in  shape,  Eq.  (87)  reduces  to 


rr.  ^  n-  ^ 

T.  0-  =  T  0- 
t  t  8  8 


(88) 


where  <r  is  the  rms  value  of  the  random  vibration.  Noting  that  the  ampli¬ 
tude  of  a  sine  wave  is  proportional  to  its  rms  value,  it  is  seen  that  Eqs. 
(86)  and  (88)  provide  the  same  result,  which  may  be  written  as 

0  _ 

In 

It  remains  to  establish  an  appropriate  value  for  b. 
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The  value  of  b  for  aircraft  materials  varies  over  a  wide  range  de¬ 
pending  upon  the  alloy,  mean  stress,  etc.  For  standard  alloys  like 
7075-T6  subjected  to  a  stress  load  with  zero  mean,  the  value  of  b  usu¬ 
ally  falls  between  5  and  8  [45] .  A  plot  of  the  service  to  test  exposure 
time  ratio  versus  the  test  to  service  level  ratio  is  presented  for  these 
bounding  values  of  b  in  Figure  23. 

It  is  appropriate  at  this  time  to  reconsider  a  basic  assumption 

used  to  arrive  at  Eq.  (76).  Specifically,  Eq.  (76)  and,  correspondingly, 

the  relationship  in  Eq.  (89),  ignore  the  existence  of  an  endurance  limit 

for  the  material  being  tested.  In  practice,  it  is  quite  possible  that  the 

environmental  levels  in  service,  indicated  by  <r  ,  would  be  below  the 

s 

endurance  limit.  This  means  that  the  service  environment  would  produce 
far  less  fatigue  damage  than  would  be  predicted  by  Eq,  (76).  If  Eq.  (87) 
were  then  applied  to  arrive  at  a  test  level  with  reduced  exposure  time, 
the  test  level  (t^  might  be  above  the  endurance  limit.  This  could 
feasibly  cause  a  failure  which  would  not  occur  in  service.  However,  the 
opposite  is  not  true;  that  is,  the  use  of  Eq.  (89)  will  not  eliminate  a 
failure  in  test  which  would  actually  occur  in  service.  Hence,  the  fact 
that  the  endurance  limit  has  been  ignored  should  cause  Eq.  (89)  to  always 
produce  conservative  results. 

9.  3.  2  Peak  Value  Criterion 

For  the  case  of  random  vibration  with  a  peak  value  criterion,  the 
relationship  between  the  test  and  service  levels  and  durations  is  obtained 
by  equating  the  probability  of  exceeding  a  specific  level  L  =  K(r  using 
Eq.  (85).  Since  the  same  component  is  involved  in  both  cases,  it  follows 
that 
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which  reduces  to 


(91) 


The  problem  now  is  to  establish  an  appropriate  value  for  L. 

The  vibration  machines  used  for  random  vibration  tests  are 
usually  equipped  with  limiting  circuits  to  protect  the  vibration  machines 
from  damaging  loads.  These  limiting  circuits  restrict  the  peak  accelera 
tions  delivered  by  the  machines  to  no  more  than  3  times  the  rms  value 
of  the  vibration.  Hence,  in  test,  the  value  of  L  is  fixed  at  L  =  3(r^. 
Substituting  this  relationship  into  Eq.  (91)  yields 


(92) 


A  plot  of  the  service  to  test  exposure  time  ratio  versus  the  test  to 
service  level  ratio  using  the  peak  value  criterion  in  Eq.  (92)  is  presented 
in  Figure  23. 
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9.  3.  3  Comparison  of  Fatigue  Damage 
and  Peak  Value  Criteria 


It  is  clear  from  Figure  23  that  the  fatigue  damage  criterion  re¬ 
quires  a  greater  increase  in  level  to  account  for  reduced  exposure  time 
than  does  the  peak  value  criterion.  For  example,  if  a  1000  hour  service 
environment  were  to  be  simulated  by  a  1  hour  test,  the  fatigue  damage 
criterion  indicates  the  test  level  should  be  from  7.  5  to  12  dB  higher  than 
the  service  level,  while  the  peak  value  criterion  calls  for  a  test  level 
that  is  only  4  dB  higher.  Hence,  if  the  peak  criterion  were  used  to  select 
the  test  level,  the  test  would  bo  less  likely  to  produce  a  fatigue  failure 
than  the  service  environment.  On  the  other  hand,  if  the  fatigue  damage 
criterion  is  used  to  select  the  test  level,  the  test  would  be  more  likely 
to  cause  an  extreme  value  failure  than  the  service  environment.  To 
assure  that  the  test  will  be  conservative,  it  appears  desirable  to  use  the 
fatigue  damage  criterion  for  selecting  appropriate  test  levels  and  dura¬ 
tions.  See  Section  9.  7  for  further  comments  on  this  subject. 


9.  4  APPLICATION  TO  NONSTATIONARY 
ENVIRONMENTS 

Now  consider  the  problem  of  simulating  a  nons\tationary  vibration 
service  life  with  a  stationary  vibration  test.  Assume  the  vibration 
service  life,  whether  random  or  periodic,  can  be  described  by  a  proba¬ 
bility  distribution  function  for  rms  values,  as  illustrated  in  Figure  24. 
Note  that  the  abscissa  of  Figure  24  can  represent  the  rms  value  of  the 

M 

vibration  in  any  desired  narrow  frequency  increment.  Hence,  the  total 
vibration  profile  would  be  described  by  a  family  of  distribution  functions 
covering  the  frequency  rahge  of  interest. 
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Figure  24.  Illustration  of  Distribution  Function  for  Aircraft  Service 
Vibration  Environment 

Let  P(<r)  denote  the  distribution  function  for  the  vibration  profile.  The 
probability  density  function  for  the  vibration  profile  is  then  given  by 

(93) 

and  will  appear  as  illustrated  in  Figure  25. 
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RMS  VALUE  OF  SERVICE  VIBRATION 


Figure  25.  Illustration  of  Density  Function  for  Aircraft  Service 
Vibration  Environment 


In  Figures  24  and  25,  a-  denotes  the  maximum  rms  vibration  which 

m 

occurs  during  the  service  environment.  In  practice,  the  maximum  level 
will  correspond  to  the  highest  level  predicted  for  the  various  flight  con¬ 
ditions. 

I 

Let  T  denote  the  exposure  time  at  the  maximum  level  cr  (or 
some  other  arbitrary  level)  to  produce  the  same  fatigue  damage  as  the 
total  service  vibration  environment.  Then,  by  equating  damage  based 
upon  the  definition  in  Eq.  (88). 
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(94) 


T 

m 


b 

cr 

m 


where  cr.  is  the  rms  level  and  At  is  the  time  duration  at  the  ith  vibra- 
1 

tion  level.  Using  Eq.  (93), 


At.  =  T  AP((r.)  =  T  p(ir.)  A(r. 

1  S  1  S  I  1 


(95) 


where  T  is  the  total  environmental  exposure  time  in  service.  Sub¬ 
stituting  Eq.  (95)  into  Eq.  (94),  and  letting  the  increments  become  very 
small,  it  follows  that 


T 

J  O’  P(o’)  dor  (96) 

(T  0 

m 


Equation  (96)  gives  the  desired  exposure  time  T  for  vibration  with  the 

m 

maximum  service  level  a ^  that  will  produce  the  same  damage  as  the 

nonstationary  service  environment.  For  example,  assume  that  b  =  8 

and  that  the  vibration  occurs  during  service  with  equal  probability  at  any 

given  level  between  0  and  <r  .  The  density  function  for  this  case  would 

m 

be  uniform  with  a  density  of  l/cr  between  0  and  or  .  From  Eq.  (96) 

m  m 


(97) 
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Hence,  the  required  exposure  time  at  the  maximum  level  would  be  only 

1/9  the  total  service  exposure  time.  Note  that  once  T  has  been  de- 

m 

termined,  the  exposure  time  can  be  further  reduced  for  testing  purposes 

using  E(|.  (89),  where  T  and  cr  are  substituted  for  T  and  o'  . 

mm  s  s 

9.  5  DEVELOPMENT  OF  VIBRATION 
PROFILE  DATA 

The  task  of  developing  vibration  profile  information  for  aircraft 
of  concern  to  AFFDL  should  be  accomplished  in  bvo  steps.  The  first 
step  is  to  estimate  the  fractional  portion  of  the  total  aircraft  service  life 
spent  in  various  pertinent  flight  conditions.  The  second  step  is  to  predict 
a  conservative  limit  for  the  vibration  environment  during  each  of  the  per¬ 
tinent  flight  conditions.  By  ranking  the  various  flight  conditions  in  the 
order  of  increasing  vibration  level,  the  distribution  of  service  vibration 
levels  is  obtained  directly,  as  illustrated  in  Figure  26. 


PREDICTED  VIBRATION  LEVEL  ,  <7 

Figure  26.  Illustration  of  Flight  Vibration  Profile 
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Given  the  design  service  life  ,  the  total  time  required  to  produce 
equivalent  damage  with  a  stationary  vibration  at  the  maximum  service 
vibration  level  is  given  by  Eq.  (96). 

Consider  now  the  problem  of  estimating  the  fractional  portion  of 
time  spent  in  various  flight  conditions.  This  should  be  done  for  each 
new  aircraft  of  interest  based  upon  analytical  evaluations  of  design 
mission  requirements.  Guidance  may  be  obtained  from  empirical  flight 
data  for  previous  aircraft  of  similar  design  and  purpose.  The  flight 
data,  however,  must  be  in  a  form  which  is  compatible  with  the  require¬ 
ments  of  the  vibration  prediction  model.  For  example,  the  flight  vibra¬ 
tion  of  jet  bombers  and  fighters  is  predicted  by  Eq.  (67)  as  a  function 
of  dynamic  pressure,  q.  Hence,  the  empirical  flight  data  used  to 
establish  vibration  profiles  must  be  in  terms  of  the  fractional  portion 
of  time  spent  at  various  values  of  q.  It  should  be  noted  that  the  con¬ 
ventional  aircraft  flight  and  maneuver  loads  documents  issued  by  AFFDL 
do  not  generally  include  this  type  of  data.  These  loads  documents  do 
include  distributions  for  the  time  spent  below  arbitrary  altitudes  and  air¬ 
speeds,  but  such  data  cannot  be  extrapolated  to  distributions  for  the  time 
spent  below  arbitrary  dynamic  pressures  without  making  unreasonable 
assumptions.  It  is  strongly  recommended  that  AFFDL  acquire  and  in¬ 
clude  such  data  in  future  flight  loads  documents. 

Another  approach  to  the  estimation  of  times  spent  in  various  flight 
conditions  is  to  use  the  flight  condition  identifications  obtained  during 
AFFDL  flight  vibration  surveys.  If  it  is  assumed  that  vibration  measure¬ 
ments  are  made  at  random  during  each  survey  flight,  then  the  fractional 
portion  of  time  spent  in  any  given  flight  condition  may  be  estimated  from 
the  ratio  of  the  number  of  measurements  identified  with  that  flight  con¬ 
dition  to  the  total  number  of  measurements.  This  is  done  for  three 
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aircraft  (B-58,  F-IOIA,  and  H-37)  representing  aircraft  Groups  3,  5, 
and  10  in  Table  27.  If  the  flights  for  vibration  surveys  are  similar  to 
the  operation  missions  for  that  aircraft,  then  data  of  the  type  presented 
in  Table  27  may  be  used  to  establish  vibration  profiles. 

To  illustrate  the  ideas  in  this  section,  assume  a  specification  is 
required  for  the  vibration  testing  of  components  to  be  installed  in  a 
fighter  type  aircraft  similar  to  the  F-IOIA.  Further  assume  that  the 
vibration  levels  in  each  frequency  increment  for  the  structural  zone  of 
interest  are  predicted  using  the  models  suggested  in  Section  8.  1.  Let 
the  results  for  the  various  contributing  sources  in  octave  band  centered 
at  250  cps  (177  to  354  cps)  be  as  follows. 

-2  2 

10  g  /cps 

10  q  g  /cps 

in-8  2  2. 

10  q  g  /cps 

0 
0 

For  simplicity,  it  is  assumed  the  contributions  of  engine  shaft  rotation 
and  gunfire  (G,,  _)  are  negligible  in  this  frequency  increment.  Noting  that 
the  bandwidth  for  the  frequency  increment  is  177  cps,  the  rms  value  for 
the  vibration  during  each  general  flight  phase  is  given  by  Eq.  (63)  as 

Takeoff;  =  [(177)  10  J  =  1.  3  g 

r  -9  2l ^  -4 

Clean  Flight;  =  [  (177)  10  q  J  =  4.  2  x  10  q  g 

Flight  with  Deployment  -3 

of  Aerodynamic  Devices;  =  [(177)(10^+10  )qj  -1.4x10  qg 
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Flight  Condition 


Percent  of  Time  Spent 
in  Given  Flight  Condition 


Flight  Condition 


Description 


B-58  F-IOIA 


Auto-Rotation  (Helicopters) 

Hover  (Helicopters) 

Rearward  Flight  (Helicopters) 

Side  Flight  (Helicopters) 

High  Speed  Stop  (Helicopters) 

Grouiid  Runup  with  One  or  More 
Engines  Out 

Climb  with  Afterburner 

Cruise  with  Afterburner 

Cruise  with  Afterburner  and  Speed- 
brakes 

Cruise  with  Afterburner  and  Flaps 
Extended 

Ground  Runup  with  Afterburner 

Cruise  with  Flaps  and  Gear  Extended 

Cruise  with  Speedbrakes  and  Gear 
Extended 

Cruise  with  Speedbrakes,  Flaps,  and 
Gear  Extended 

Cruise  with  Afterburner,  Speedbrakes, 
and  Flaps  Extended 


2 


From  Table  27,  the  flight  conditions  which  fall  into  each  of  these  three 
flight  phases,  and  the  total  fractional  portion  of  time  within  each  phase 
are  as  follows. 


Takeoff  (02,  04)  =  8.  7% 

Clean  Flight  (06,  07,  21.  23,  48)  =  45.  2% 

Flight  with  Deployment  of 

Aerodynamic  Devices  (08-17,  22,  24,  28,  49,  50, 

55,  56,  61)  =  38.  8% 


The  remaining  flight  conditions  (01,  25-27,  44)  which  account  for  7.  3% 

of  the  service  life  will  be  ignored. 

A  probability  distribution  for  q  during  flight  with  and  without 

the  deployment  of  aerodynamic  devices  must  now  be  estimated.  For 

simplicity,  assume  q  is  uniformly  distributed  over  the  range  from  0 

to  1000  for  clean  flight,  and  from  0  to  300  for  flight  with  deployment  cf 

aerodynamic  devices.  The  maximum  vibration  levels  then  occur  during 

takeoff  where  or  =0",  =  1.  3  g.  The  time  duration  for  vibration  at  this 

m  to  ® 

level  required  to  produce  fatigue  damage  equivalent  to  the  service  en¬ 
vironment  is  given  by  Eq.  (96).  Assuming  b  =  8, 


T 

m 


(1.3) 


0. 

■ 


087(1.  3)®  +  0.  452(0.  21)®/9  +  0.  388(0.  42)®/9 


Hence,  for  this  illustration,  nearly  all  of  the  damage  occurs  during 
ground  runup  and  takeoff,  which  constitute  less  than  9%  of  the  total 
service  environment. 

To  carry  the  illustration  a  step  further,  assume  the  aircraft  is 
designed  for  a  service  life  of  2000  hours.  Based  upon  a  fatigue  damage 
criterion,  the  required  duration  for  a  vibration  test  at  the  predicted 
takeoff  vibration  level  would  be  174  hours  to  simulate  the  service  en¬ 
vironment.  However,  using  Eq.  (89)  or  Figure  23,  the  test  duration 
could  be  reduced  to  0.  68  hour  (41  minutes)  by  increasing  the  test  level 

to  2  times  the  predicted  takeoff  level.  Hence,  a  random  vibration 

2 

test  with  a  level  of  0.  015  g  /cps  in  the  frequency  range  from  177  to 
354  cps,  and  a  duration  of  41  minutes  would  be  appropriate  for  this 
illustration. 

The  above  illustration  considered  the  vibration  environment  in 
only  one  frequency  increment.  In  practice,  similar  calculations  would 
be  required  for  all  frequency  increments  covering  the  range  of  the  data. 
The  results  in  terms  of  the  required  test  duration  at  maximum  service 
levels  may  differ  from  one  frequency  increment  to  the  next.  This  is 
clearly  undesirable  from  the  testing  viewpoint.  Hence,  the  levels  in 
each  frequency  increment  should  be  modified  using  Eq.  (89)  to  produce 
the  same  required  test  duration  for  all  frequencies.  It  follows  that  the 
resulting  power  spectrum  for  the  test  may  not  be  similar  to  any  of  the 
power  spectra  for  various  phases  of  the  predicted  vibration  environment. 
The  power  spectrum  for  the  test  will,  in  effect,  represent  a  properly 
weighted  composite  for  the  varying  power  spectra  of  the  vibration  which 
occurs  during  the  service  life  of  the  aircraft. 
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9.  6  SELECTION  OF  TEST  LEVEL 

PERCENTAGE  POINTS 

The  levels  specified  for  a  vibration  test  should  never  be  increased 
to  exceed  the  maximum  predicted  service  environment  unless  a  corre¬ 
sponding  reduction  in  test  duration  is  made  in  accordance  with  Eq.  (89). 
Specifically,  no  "safety  factors"  should  be  added  to  the  test  levels,  as 
determined  by  the  procedures  suggested  herein.  Any  "safety  factor" 
needed  to  account  for  uncertainties  in  the  environmental  levels  will  have 
been  included  in  the  ^.rediction  procedure,  as  discussed  in  Section  8. 

The  selection  of  a  percentage  point,  a,  for  the  predictions  in  effect 
establishes  the  desired  degree  of  conservatism  for  the  predictions  and, 
hence,  the  resulting  vibration  test  specification.  The  smaller  the  value 
of  a,  the  greater  the  conservatism.  An  unresolved  issue  is  the  selec¬ 
tion  of  an  appropriate  value  for  a . 

In  the  past,  the  selection  of  percentage  points  for  vibration  pre¬ 
dictions  has  been  somewhat  arbitrary.  The  percentage  points  recom¬ 
mended  by  various  investigators  [2],  [5],  [8],  [33]  have  generally 
ranged  from  1%  to  5%,  providing  prediction  limits  which  hopefully  would 
exceed  from  95%  to  99%  of  the  anticipated  vibration  levels  in  a  given 
structural  zone.  Although  percentage  points  in  this  range  are  intuitively 
reasonable,  it  would  be  desirable  if  the  percentage  point  selection  could 
be  made  on  a  more  scientific  basis.  One  method  of  approaching  this 
problem  is  through  the  use  of  statistical  decision  theory  ("cost"  minimi¬ 
zation  procedures).  Such  a  method  is  suggested  in  [46],  and  outlined 
$ 

below. 

Consider  an  aircraft  component  which  must  be  vibration  tested  as 
part  of  its  qualification  for  service  use.  Assume  there  are  only  two 
possible  test  results;  pass  or  fail.  If  the  component  passes  the  test, 
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there  is  a  possibility  that  it  will  fail  in  service.  The  "cost"  associated 
with  this  possibility  is 


Ca  =  C„  Pj,  (98) 

where  is  the  "cost"  associated  with  a  service  failure,  and 
the  probability  of  a  service  failure.  If  the  component  fails  the 
must  be  redesigned  to  pass  the  test.  The  "cost"  associated  wi 
possibility  is 


teVt^it 
th  this 


C 


f 


(99) 


where  is  the  anticipated  "cost"  of  a  redesign  needed  to  pass  the  test. 
It  follows  that  the  total  "cost"  associated  with  the  test  is 


C,  =  P(L)  C  + 

t  cl 


[.-PC.; 


P(L)|  C^MCqPj,)  P(L) 


-  1  -  P(L) 


(100) 


where  P(L)  is  the  probability  of  passing  the  test  at  level  L.  Clearly, 
P(L)  is  a  function  of  the  percentage  point  used  to  establish  the  test  levels. 
The  desired  result  is  to  select  that  value  of  a  which  will  minimize  the 
total  "cost,  "  C^,  in  Eq.  (100).  From  [46],  a  lower  bound  on  the  optimum 
value  for  a  is  given  by  the  simple  expression 


0=1- 
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where  n  is  the  total  number  of  units  for  the  component  design  in 
question  to  be  manufactured  and  installed  in  service,  and  C^/C^  is 
assumed  to  be  less  than  unity. 

The  deviation  of  Eq.  (101)  involves  several  assumptions  which 
might  be  questionable  in  some  cases,  as  discussed  in  [46]  .  Further¬ 
more,  the  "cost"  terms  may  be  difficult  to  estimate  in  many  if  not 
all  cases.  Nevertheless,  it  is  believed  that  the  use  of  Eq.  (101)  with 
even  rudimentary  "cost"  estimates  provides  a  better  basis  for  select¬ 
ing  a  percentage  point  for  vibration  predictions  than  a  simple  intuitive 
guess.  Note  that  the  "cost"  terms  in  Eq.  (101)  may  be  expressed  in 
any  units  desired  (dollars,  relative  importance,  etc.  ).  Also  note  that 
the  "cost"  terms  appear  as  a  ratio.  This  greatly  simplifies  the  problem 
of  estimating  these  terms,  since  absolute  values  are  not  required.  For 
example,  it  might  be  unreasonable  to  estimate  the  "cost"  of  forcing  a 
pilot  to  eject  because  of  a  service  failure,  or  the  "cost"  of  delaying  a 
schedule  for  redesign  because  of  a  test  failure.  It  may  not  be  un¬ 
reasonable,  however,  to  estimate  the  relative  importance  of  these  two 
events. 

To  illustrate  Eq.  (101),  assume  a  component  is  to  be  tested  where 
the  ratio  of  importance  of  a  test  failure  to  service  failure  is  assessed  to 
be  0.  5.  Further  assume  that  16  units  are  to  be  manufactured  and  in¬ 
stalled  for  service  use.  The  prediction  limit  used  to  derive  the  vibra¬ 
tion  test  level  should  be 


100  (1  -  a)% 


100  [  1  -  0.  5] 


95.  7% 
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If  the  ratio  of  importance  had  been  0.  1,  the  test  level  wt'uld  be 


100  (1  -  o)%  =  100  J 


0.  1 


1/16 


99.  3% 


Hence,  as  the  anticipated  "cost"  of  a  service  failure  increases  relative 
to  the  anticipated  "cost"  of  a  test  failure,  the  required  test  level  in¬ 
creases.  If  the  ratio  of  importance  had  remained  0.  5,  but  the  number 
of  units  had  been  64,  the  test  level  would  be 


100  (1  -  a)%  =  100 


[1-0.5]"^^ 


98.  9% 


Hence,  as  the  number  of  units  to  be  manufactured  and  used  in  service 
increases,  the  required  test  level  increases. 

9.  7  CONCLUDING  COMMENTS  ON 

VIBRATION  TEST  SPECIFICATIONS 

The  developments  in  Section  9  are  presented  only  to  provide  some 
analytical  guidelines  which  might  be  employed  to  develop  improved  pro- 
ceduies  for  the  derivation  of  vibration  test  specifications.  Many  prac¬ 
tical  matters  of  considerable  importance  have  been  omitted  in  the  develop¬ 
ments  for  clarity  and  simplicity.  For  example,  the  endurance  limit  of 
structural  materials  and  the  nonlinear  response  properties  of  structural 
assemblies  are  ignored  in  the  derivation  of  Figure  23.  The  fact  that  many 
types  of  failures  may  not  be  related  to  fatigue  damage  is  also  ignored. 
Furthermore,  no  consideration  is  given  to  the  equipment  loading 
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(structural  impedance)  problem.  Finally,  the  suggested  guidelines 
inherently  assume  that  random  vibration  tests  will  be  used  to  simulate 
the  random  portions  of  the  environment,  while  sinusoidal  tests  will  be 
used  to  simulate  the  periodic  portions  of  the  environment.  In  practice, 
the  luxury  of  such  simplifying  assumptions  may  not  be  tolerable. 

Hence,  a  brief  discussion  of  these  matters  is  warranted. 

9.  7.  1  Endurance  Limit  Problem 

The  endurance  limit  problem  is  often  dealt  with  by  assuming  the 

service  environment  is  limited  to  that  duration  necessary  to  deliver 

6 

5x10  stress  cycles  to  the  test  item  of  interest.  The  theory  behind 
this  assumption  is  as  follows  [47],  [48].  If  a  structural  material  is 
going  to  experience  a  fatigue  failure,  it  will  generally  occur  before 
5  X  10^  stress  cycles  have  been  accumulated.  In  other  words,  the 
break  in  the  idealized  S  -  N  curve  shown  in  Figv  re  21  usually  occurs 
at  about  5  x  10^  cycles.  Hence,  if  a  fatigue  damage  criterion  is  to  be 
used  to  trade  off  test  levels  against  test  durations,  there  is  no  need 
to  consider  the  service  duration  to  be  longer  than  that  time  necessary 
to  accumulate  5  x  10^  stress  cycles  in  the  test  item. 

9.  7.  2  Structural  Nonlinearity  Problem 

There  is  little  question  the  c  the  nonlinear  properties  of  actual 
structures  often  tend  to  restrict  the  occurrence  of  extreme  displace¬ 
ment  values  in  the  structural  response  [49] .  This  means  that  the  trade¬ 
offs  between  test  levels  and  test  durations  suggested  in  Figure  23,  for 
either  a  fatigue  damage  or  peak  criterion,  may  deviate  widely  from  the 
actual  requirements.  One  suggested  method  for  overcoming  this  problem 
is  to  perform  all  tests  at  the  predicted  environmental  levels,  and  then 
limit  the  test  duration  solely  on  a  basis  of  the  5x10  stress  cycle  rule  [6] . 
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9.7.3  Nonfatigue  Failure  Problem 


Many  types  of  aircraft  component  failures  have  little  or  no  relation 
to  fatigue  damage  [37].  It  follows  that  vibration  tests  performed  at  ele¬ 
vated  levels  designed  to  simulate  the  service  fatigue  life  may  produce  un¬ 
representative  failures.  This  problem  can  be  greatly  reduced  by  separat¬ 
ing  the  testing  procedure  into  a  "structural  integrity"  test  and  a  "per¬ 
formance"  test.  The  "structural  integrity"  test  would  be  performed  using 
elevated  test  levels  to  simulate  the  service  fatigue  damage  with  relatively 
short  test  durations.  However,  functional  performance  of  the  test  com¬ 
ponent  would  not  be  expected  or  required  during  or  after  the  test  at  ele¬ 
vated  levels.  A  failure  would  be  assessed  solely  on  the  basis  of  a  struc¬ 
tural  failure  of  the  component.  Using  a  new  sample  component,  the  func¬ 
tional  performance  would  then  be  evaluated  by  a  second  test  performed  at 
the  predicted  envi’^onmental  levels  with  an  arbitrary  test  duration.  The 
theory  here  is  that  functional  performance  failures  may  be  extremely 
sensitive  to  elevated  test  levels  (in  an  unpredictable  manner),  but  prob¬ 
ably  not  very  sensitive  to  exposure  time. 

9.  7.  4  Equipment  Loading  Problem 

A  major  problem  in  the  specification  and  performance  of  aircraft 
component  vibration  tests  is  the  equipment  loading  problem.  For  the 
case  of  aircraft  components  which  are  relatively  heavy  compared  to 
their  supporting  structure,  there  may  be  a  significant  difference  be¬ 
tween  the  vibration  response  of  the  unloaded  structure  and  the  vibration 
of  the  structure  with  the  component  attached.  Vibration  test  specifica¬ 
tions  should  always  be  based  upon  predictions  for  the  structural  vibration 
with  the  component  of  interest  attached.  The  surface  weight  density 
term  (w)  in  the  random  vibration  prediction  models  suggested  in  Section  8 
will  partially  account  for  the  effect  of  component  loading.  The  situation 
is  further  improved  if  the  data  used  to  arrive  at  the  predictions  are 
measured  near  equipment  mounting  points.  Problems  may  still  arise. 
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however,  when  heavy  components  with  lightly  damped  resonances 
(sharply  peaked  impedance  functions)  are  tested.  The  pooling  operations 
used  to  arrive  at  the  predictions  tend  to  obscure  the  sharp  notches  in 
the  structural  response  spectra  induced  by  strong  resonant  reactions  of 
heavy  components  mounted  on  low  impedance  structures  [50].  A  com¬ 
mon  method  for  dealing  with  this  problem  is  to  permit  notching  of  the 
specified  test  levels  at  such  frequencies.  This  may  be  done  by  per¬ 
mitting  the  vibration  testing  machine  to  react  as  permitted  by  its  natural 
impedance  function,  by  analytical  corrections  based  upon  mounting  point 
impedance  measurements  for  the  supporting  structure,  or  by  specifying 
limits  on  the  response  of  the  component.  See  [51],  [SZ]  for  further  dis¬ 
cussion  of  this  problem. 

9.  7.  5  Environmental  Simulation  Problem 

As  noted  previously,  the  developments  in  Section  9  inherently 
assume  that  the  random  and  periodic  portions  of  the  vibration  environ¬ 
ment  will  be  fimulated  in  test  with  random  and  sinusoidal  excitations, 
respectively.  At  the  present  time,  this  procedure  is  not  always  em¬ 
ployed.  Specifically,  random  vibration  tests  are  often  omitted  in  the 
development  of  components  for  aircraft  applications.  Certs  inly  an  im¬ 
portant  reason  for  this  must  be  the  fact  that  MIL-STD-810B  does  not 
specify  random  vibration  tests  for  aircraft  components.  It  is  strongly 
recommended  that  any  future  Military  Standard  Specification  for  the 
environmental  testing  of  aircraft  components  include  provisions  for 
random  vibration  tests. 

A  second  simulation  problem  evolves  from  the  fact  that  periodic 
components  in  the  environment  (and  sometimes  the  random  components 
as  well)  are  usually  simulated  in  test  by  a  "frequency  sweep"  type  of 
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sinusoidal  excitation,  rather  than  by  a  series  of  steady  state  sine  waves 
with  fixed  frequencies  corresponding  to  the  anticipated  environmental 
frequencies.  Frequency  sweep  sinusoidal  testing  introduces  problems 
related  to  sweep  rate  and  total  sweep  duration.  Fortunately,  these 
problems  have  been  the  subject  of  considerable  past  study,  as  sum¬ 
marized  in  [47] ,  [48]. 
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10.  EXTENSIONS  TO  ACOUSTIC  ENVIRONMENTS 


The  prediction  of  acoustic  environments  might  be  approached  using 
the  same  general  statistical  procedures  suggested  for  the  prediction  of 
vibration  environments  in  Section  8.  In  at  least  two  respects,  the  develop¬ 
ment  of  appropriate  models  is  somewhat  easier  for  the  case  of  acoustic 
predictions.  First,  the  spatial  distribution  of  acoustic  levels  within  a 
given  structural  compartment  is  considerably  less  dispersed  than  the  distri¬ 
bution  function  for  vibration  levels  over  the  structure.  Second,  because  of 
high  modal  density  and  reverberation  effects,  the  power  spectra  for  random 
acoustic  levels  within  a  compartment  are  generally  smoother  in  shape  than 
the  power  spectra  for  the  structural  random  vibration  levels.  In  two  other 
respects,  however,  the  acoustic  prediction  problem  is  more  difficult. 

First,  the  acoustic  levels  in  a  given  aircraft  compartment  are  heavily  de¬ 
pendent  upon  the  details  of  the  compartment  soundproofing  materials  and 
their  installation.  Second,  a  major  portion  of  the  acoustic  environment  is 
contributed  by  the  airconditioning  system,  and  is  heavily  dependent  upon 
the  details  of  the  airconditioning  system  outlets. 

In  spite  of  the  above  not'^d  problems,  the  possibility  of  developing  an 
extrapolation  type  prediction  procedure  for  internal  acoustic  noise  in  aircraft 
should  be  investigated.  It  is  suggested  that  this  effort  be  pursued  using 
models  of  the  form  suggested  for  vibration  environments  in  Section  8,  with 
one  exception.  For  the  acoustic  prediction  problem,  an  additional  constant 
term  should  be  added  to  the  various  expressions  in  Eq.  (63)  to  account  for 
air  conditioning  noise.  For  example,  the  model  for  high  speed  flight  should 
be 

°'«hs  =  "  °<«r2  " 

where 

G(f)  =  Aq  regression  coefficient  to  account  for  the 
contribution  of  airconditioning  noise 

G(f)^2  “  definition  given  by  Eq.  (67) 

G(f)j,  -  definition  given  by  Eq.  (70)  ; 

d  j 

( 

1 


227 


11.  CONCLUSIONS 


The  studies  herein  define  a  procedure  for  extrapolating  data 
measured  on  previous  aircraft  to  derive  vibration  and  acoustic  design 
criteria  and  vibration  test  specifications  for  future  aircraft.  The 
suggested  procedure  requires  only  rudimentary  information  concerning 
the  future  aircraft  and  its  flight  conditions.  Hence,  it  may  be  used  to 
establish  preliminary  specifications  and/or  criteria  early  in  the  con¬ 
ceptual  design  phase.  The  basic  steps  in  the  suggested  procedure  are 
outlined  in  Figure  27.  Each  step  is  identified  with  numbers  corre¬ 
sponding  to  sections  of  this  report  which  relate  to  the  noted  step. 

The  practicality  of  each  suggested  analysis  technique  in  Figure  27 
has  been  verified  using  currently  available  AFFDL  vibration  data.  The 
current  data,  however,  are  not  entirely  adequate  for  all  steps  in  the  over 
all  procedure.  Further  data  acquisition  and  analysis  in  accordance  with 
the  techniques  suggested  herein  will  be  required  to  fully  implement  the 
suggested  procedure. 
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